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Geophysics. — The Earth's crust deformation in the East Indies. 
(Provisional Paper.) By F. A. VENING MEINESZ. 


(Communicated at the meeting of February 24, 1940.) 


§ 1, Introduction. 


For the further investigation of the results of the gravity expeditions 
with submarines of the Netherlands Navy, the writer has prepared new 
tables for regional isostatic reduction; he has applied the system published 
in the Bulletin Géodésique N°. 29, 1931, to five different degrees of 
regionality. This system distributes the isostatic compensation of each 
topographic element over a wider area than the area of the element itself 
and the distribution in a horizontal sense follows the rule that the amount 
of compensation is proportional to the vertical displacement brought about 
by the bending of the Earth’s crust under the weight of the topographic 
elements. For each element it adopts the bending-curve of the solution of 
HERTZ for the bending of an elastic plate floating on a fluid and sub- 
mitted to a concentrated surface load. Introducing the length / given by 
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where 


F = modulus of elasticity of the Earth’s crust, 


m == coefficient of Poisson, 
h = thickness of the Earth’s crust, 
6,— 0) = difference of the specific densities of the plastic substratum 


and the topography, 


the central part of the bending-curve of HERTZ was used up to a distance 
of 2.905 J from the concentrated load; the outward waves of small ampli- 
tude at greater distances from the center were neglected. The adopted 
curve is represented by fig. 1. So for each topographic element the isostatic 
compensation is assumed to be distributed in a horizontal sense up toa 
distance of t= 2.9051, having a maximum specific density in the center 
and this density falling to zero towards the border t according to the curve 
of fig. 1. The tables have been made for five different values of Lae. 
10, 20, 40, 60 and 80 km corresponding to values of the outer boundary t 
of 29.05 km, 58.10 km, 116.20 km, 174.30 km and 232.40 km. They have 
been published in the Bulletin Géodésique NO. 63, 1940 and details about 
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these laborious computations may be found in “Verhandelingen der Kon. 
Akad. v. Wetenschappen, Afd, Natuurkunde”’, Iste Sectie, DI. XVII, 
N°, 3, The Netherlands Geodetic Commission has defrayed the expenses 
of the computations. 
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Fig. 1. Curve of the distribution of the regional isostatic compensation. 


In these fundamental tables the compensation has been assumed to be 
distributed uniformly over a depth H below sea-level and the tables 
comprise all values of H from 0 to 60 km. From these tables other tables 
have been derived for practical use in which the compensation extends 
from a depth of 30km down to a depth of 30 + 4.45, s being the height of 
the topography above sea-level. This last assumption corresponds to that of 
the Airy-Heiskanen tables for local isostatic reduction based on the principle 
of a crust of a normal thickness of 30 km and a density of 2.67 floating on 
a plastic substratum of a density of 3.27, each element being in hydrostatic 
equilibrium independent of its surrounding and sinking till a local root at 
the lower boundary of the crust is formed reaching from a depth of 30 km 
down to 30 + 4.45 .s km. The new tables, therefore, adopt the same system 
of hydrostatic equilibrium of the crust but they assume the root at the 
lower boundary of the crust to be broader according to a bending of the 
crust, instead of its sinking locally under the load of the topography. 
Incase we have a submarine topography, we have to introduce a value of s 
with a negative sign, the topography being now below sea-level, and we 
have to put s== —0.615 d, if d is the depth, in order to take into account 
that the deficiency of mass has a specific density of 2.67—1.03; 1.03 has 
been assumed for the specific density of sea-water. The tables will be 
published in the course of this year. 

The Netherlands Geodetic Commission has used these tables for the 
isostatic reduction of all the gravity values obtained at sea. This paper 
contains a provisional publication of some conclusions derived from the 
results of these reductions for the East Indies. For each station six 
anomalies have been deduced, five according to the five degrees of 
regionality that have been mentioned and one for local compensation under 
the same conditions of vertical distribution. These last values as well as 
those for small regionality give an irregular picture of the distribution of 
the anomalies over the Archipelago; the anomalies for the three largest 
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degrees of regionality show simpler features. So probably a large degree 
of regionality fits better the actual conditions in the Earth's outer layers. 
The most regular picture seems to be given by the anomalies computed 
with the tables for 160 km, which corresponds to an outer radius ¢ of 
the compensation of 174.3 km. They are represented by the map of this 
paper which shows the isanomales for intervals of 50 milligals and as an 
example two gravity profiles are given in fig. 2a and b. The maps and 
profiles for 1—= 40 km and for /== 80 km also show acceptable distributions. 
A later publication by the Netherlands Geodetic Commission will contain 
the maps of the anomalies and a great many profiles for all the values 
of Land will thus allow a comparizon. The conclusions mentioned in this 
paper will be derived from the present map but nearly all of them would 
remain valid when the maps for /=40 km or = 80 km had been used. 


§ 2, The Map of the Regional Isostatic Anomalies ({—=60 km, 


In considering this map the writer thinks that there are several reasons 
why an attempt is justified to explain the anomalies by taking as a base 
the simple mechanical picture of a rigid Earth’s crust floating on a plastic 
substratum and subject to deforming forces; the present publication in 
fact is such an attempt. Of course he is perfectly aware that the real 
conditions in the Earth are more complicated. In the first place he will 
assume only one layer of constant density floating on a denser substratum 
of a specific density of 0.6 more than that of the crust. This is certainly 
not true as seismic and other evidence points to at least two and probably 
more rigid layers of increasing density downwards, but this simplification 
does not make a difference in principle for the main conclusions; 
they can be adapted without difficulty to the more complicated picture of 
more than one rigid layer. In the second place changes of state, e.g. from 
the crystalline to the amorphic, or differentiations and chemical reactions 
may occur. There is no doubt that these phenomena may affect the main 
occurrence; a further research may decide how far this will be the case. 
The object of this study is to see how far the above simplified picture will 
go in explaining the facts. 

As a first reason justifying this attempt the writer may mention the 
fact that the regional isostatic reduction, based on this assumption, un- 
doubtedly simplifies the anomaly map and so the way the topography is 
compensated is obviously in harmony with it. He must refer to the future 
publication of all the maps and profiles for proving this statement. For a 
second reason the writer may bring forward the success of the buckling 
hypothesis, which he advanced in 1930 for explaining the main feature of 
the anomaly map, the narrow belts of strong negative anomalies through 
the Archipelago, Although other attempts for an explanation have been 
made, it still remains the one that has found the widest support. It is based 
on the same mechanical picture and we shall now have to consider it a 
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little more in detail as it is in fact one of the principal parts of our attempt 
at an explanation of the anomaly field. 

As we shall point out in the next § an elastical plate floating on a fluid 
of greater density and submitted to a horizontal compressive force, will 
remain straight as long as this force is below a certain limit Do which we 
may call the buckling limit. When this limit is passed, waves will form and 
a slight increase of the force will then suffice for these waves to increase 
further in amplitude till one of them buckles. The writer assumes this to 
have occurred in the geosyncline in the East Indies under the effect of the 
tectonic forces. This corresponds to the great folding and overthrusting 
that has been found in the islands in the belt of anomalies, The buckling 
takes place in a downward sense and thus a great root of light crustal 
material is formed at the lower boundary of the Earth’s crust. Because of 
its displacing the denser substratum this represents a large defect of mass 
and so the strong negative anomalies may thus be explained. We may 
adopt the supposition that a downward wave buckles downwards and not 
an upward wave upwards, because of the fact that the energy needed for 
the first is proportional to the difference of density of the substratum and 
the crust, which may be estimated at 0.6, while for the latter it must be 
proportional to the density of the crust of 2.7. This theoretical consideration 
is confirmed by the beautiful experiments of Dr. PH. H. KUENEN 1), who 
compressed a plastic layer floating on water; gradually the layer forms 
waves and at a certain moment one of these waves buckles downwards. 
These experiments give a good support to our hypothesis. We may assume 
that in a further phase of the orogenetic development the upper layer of 
the crust will crumple upwards and gradually form a ridge that in its 
ultimate shape will compensate the defect of mass formed by the root. 
Thus the stage is reached of the complete mountain range with its com- 
pensation by a root of crustal material, about 44 times as large as the 
range itself, which we find in the Alps and other folded mountain ranges. 
The main point of this hypothesis, therefore, is that the root is formed 
before the range itself and that this stage is now reached in the East 
Indies. The formation of a surface ridge is already going on in the East 
Indies; in most parts the belt of anomalies coincides with a submarine 
ridge or with islands as is shown e.g. by the profiles of fig. 2a and 2b 
but this ridge is still far too small to compensate the negative anomalies 
caused by the root. The formation of this surface ridge may also be 
brought about by the rising of the buckling zone because of the read- 
justment of the isostatic equilibrium or by other causes; we shall come 
back to this point. The buckling hypothesis may also explain another 
feature, i.e. the well-known fact that an orogenetic cycle begins by the 
formation of a geosyncline; this may be interpreted as the downward 


1) Pu. H. KUENEN, The negative isostatic anomalies in the East Indies. Leids, Geol. 
Meded. VIII, 2 (1936). 
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wave which afterwards, when the compression continues, buckles inwards. 
These few words may suffice for giving an idea of the main points of 
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Fig. 2a. Gravity profile over Seroea and Tanimber Isl, Drawn line: regional 

isostatic anomalies (1 — 60 km, tf = 174.30 km), Dotted line: the same anomalies 

after subtracting the estimated effect of the root, Gravity Stations indicated by 
their number. 
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Fig. 2b. Gravity profile over Str. Bali southwards, Drawn line: regional 
isostatic anomalies ([— 60 km, t — 174.30 km), Gravity Stations indicated by 
their number, 


the buckling hypothesis. The writer should now wish to look at the distri- 
bution of the phenomenon through the Archipelago. In the first place we 
may notice on the map that the anomalies point to a compressive force 
mainly working in one direction as indicated by the arrow, and not to 
stresses working in all directions. Where the anomaly belt cuts this 
direction under an angle of at least 45° the anomalies are large, but where 
the angle is small the anomalies are likewise small, as e.g. in the profiles 
west of Sumatra, in the profile from the Banda Sea towards New Guinea 
and also in the E.W. profile over Strait Surigao, north of Mindanao; this 
last profile outside the area of the map does not show negative anomalies 
at all. As the total of the negative anomalies in a profile over the belt is 
a measure for the size of the root that has formed at that place and, there- 
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fore, also for the amount of the shortening of the crust, this shortening 
must have been smaller for these last profiles than for the profiles else- 
where. This finds an explanation by the assumption of the compression 
working in the indicated direction because in that case the component of 
the compressive force in the sense of these profiles is small, The writer 
considers it a favourable point for the mechanical picture advocated in 
this paper that these differences in intensity of the negative anomalies 
along the belt may thus be explained by a simple supposition. 

The foregoing considerations about the anomaly-belt have already been 
expressed by the writer in previous papers; in what follows he wishes to 
add a new view-point partly based on the results of the new isostatic 
reductions. These results have led viz. to a slight alteration of the isano- 
males near the island of Soemba. The anomalies west of the island of 
Timor bring about a pointed course for the isanomale of —50 milligal 
near the S.E, Cape of Soemba and this has brought the writer to draw the 
isanomales in that neighbourhood as they are given on the map. This leads 
to a new vision on the irregularity of the belt of anomalies in this area. 
While the old anomaly map showed an interruption of the belt, the new 
one only points to a shift of the axis. This is better acceptable from the 
view-point of the buckling-hypothesis; an interruption would be difficult 
to understand because the question rises how the crust could be shortened 
to both sides of the interruption without the block of Soemba taking part 
in this shortening. A shift of the axis along a line which is about in the 
direction of the compression could, however, be better explained. In this 
connection it is worth while to notice that this shift occurs exactly where 
the Australian continent begins to face the belt; for showing this the map 
has been provided with the contour line of 1000 m depth. This coincidence 
suggests a causal relation and it is easy to think of a satisfactory explana- 
tion. Considering the great amount of sedimentation at the continental 
edge we could understand that a continent can be surrounded by a belt of 
weaker crustal material and in that case the buckling axis could not continue 
from its course south of Java straight towards the east but it should have 
to shift towards that zone of weakness; in its further course to the east it 
continues to follow the continental edge. We thus get a satisfactory 
explanation of the curious twist of the island ridge and of the anomaly belt 
in this neighbourhood. Here again the mechanical view-point provides us 
with a simple hypothesis. The fact that the island of Soemba does not show 
great folding or overthrusting is in good corroboration. 

For further considerations about the buckling zone the writer wishes to 
refer to a future and larger publication, which will allow to go more into 
detail. Here other evidence regarding the mechanical picture advocated in 
this paper asks our attention. When we look at the anomaly map a perio- 
dicity in the occurrence of belts of negative and of positive anomalies 
strikes us. We notice it when drawing a line from the Celebes Sea towards 
the S.S.E.; we see here a regular succession of positive and negative 
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anomalies, more or less in the shape of belts, and the distance of successive 
features is of the order of 150 to 200 km. In a N.S. profile through Java 
we find the same: a negative belt in the northern half of Java, coinciding 
with the oil-syncline, a positive belt over the southern half of the island 
extending over a strip of the Indian Ocean, the belt of large negative 
anomalies south of it and lastly a belt of smaller positive anomalies which 
does not show in the isanomales but which may be derived from the 
figures. West of Sumatra we find a similar periodicity. And in the West 
Indies which show a far-going analogy to the East Indies and where 
likewise a belt of strong negative anomalies has been found, a similar 
periodicity may be noticed. So it looks as if it is a regularly occurring 
feature and the question arises what this periodicity may mean. 

Here again the writer thinks that the mechanical view-point may provide 
us with a better explanation than other lines of interpretation would be 
able to give; the obvious view-point is to interpret these alternating belts 
as brought about by a wave-formation in the rigid crust. The waves must 
give alternating deviations from the isostatic equilibrium of the crust, the 
upward wave causing positive anomalies and the downward one negative 
anomalies. We may find an important corroboration of this idea in the 
fact that the wave-length of 300—400 km is in keeping with that derived 
from the results of the new regional isostatic reduction; in the next § the 
writer will give a theoretical treatment of the phenomenon which will show 
the agreement. Another corroboration worth while mentioning is the fact 
that the profiles of the ocean bottom in many of the gravity profiles west 
of Sumatra and south of Java show a regular wavelike topography; the 
amplitude of the wave, i.e. the height of the topography, corresponds to 
the amplitude of the wave in the gravity anomalies when assuming no 
isostatic compensation and this is of course as it ought to be according to 
our supposition because the wave is entirely a deviation from the isostatic 
equilibrium. We do not find this same relation to the topography as clearly 
elsewhere in the Archipelago, but this may be explained by the complicated 
deformations of the surface layers and by the other surface effects of 
erosion, sedimentation and volcanological activity. The writer must again 
refer to the future publication which will contain all the gravity profiles 
for the corroboration of these statements; an instance may be seen in 
fig. 2b, a gravity profile going S.S.W. from Strait Bali. He may also refer 
to that publication for a detailed comparizon of this interpretation with the 
old attempt made in 1934 in “Gravity Expeditions at Sea’, Vol. II 1), 
where he attributed the positive anomalies in the deep basins, e.g. the 
Banda Sea and the Celebes Sea, to more or less local convection currents 
in the substratum and the belts of positive anomalies south of Java and 
west of Sumatra to an upwards drag of the crust because of the regional 
readjustment of the isostatic equilibrium of the crust over the belt of strong 


1) Publ. Neth. Geod. comm., Delft 1934, 


285 


negative anomalies, Both explanations meet with difficulties when we go 
deeper into the matter and so the writer considers the present attempt not 
only as simpler but also as better, 

The assumption of waves in the crust in the present period implies the 
assumption of a still working compressive force, This is in harmony with 
the fact that in a very recent period, in the pleistocene, folding has occurred 
in the oil-geosynclines of the Archipelago and also in E. Celebes and 
Boeton. It is also in keeping with the strong seismicity of the Archipelago. 
So we assume the crust to be subject to a compressive force in a horizontal 
sense, but we know that it is also subject to an upward force in the belt 
of strong negative anomalies because of the tendency to restore the isostatic 
equilibrium 1). We have, therefore, to take up the problem of the deforma- 
tions of the crust by this combination of forces. We shall take this up in 
the next § and we shall have a look at the consequences of these results 
for the recent geology in the ensuing §. 


§ 3. Formulas for the deformation of a floating elastical plate under 
the effect of a vertical force and of horizontal compression. 


We assume the plate to have constant thickness and infinite dimensions. 
We further assume the problem to be two-dimensional and we shall here 
only treat of the simplest case of only one concentrated vertical force in 
upward direction, We put the origin of the coordinates in the point where 
the force acts and we take the horizontal axis as the X axis and the 
vertical one as the Y axis (fig. 3). We assume the deformations of the 
plate to be small enough that we may neglect the second and higher powers 
of dy/dx. In these assumptions we have neglected the curvature of the 
Earth, but it is easy to see that its only effect is an upward shift of the 
plate of D/6,9R, where D is the horizontal compressive force, 6, the 
specific density of the substratum and R the Earth’s radius. Using the 
same letters and indications as introduced on the first page of this paper 
and introducing the quantities 1, and Dp given by 
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we find the following differential equation for the displacement y of the 
plate in vertical sense 

1D) 2d’y 
Dy, | dx? 
We see at once that for D=Dy the solution shows a simple sine-curve 
and a further investigation by means of no longer neglecting the second 
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1) It is likewise subject to the vertical forces exerted on the crust by the topographic 
features but these forces are smaller and we shall neglect them here. 


286 


and higher powers of dy/dx shows that the amplitudes are quickly 
increasing when D increases slightly. So we see that Do is the buckling- 
limit of the Earth’s crust, ie. when D reaches this limit the deformations 
are growing so large that one of the waves will give way and a root will 
form. In the further investigations we shall suppose D< Dp and for 
simplification we shall introduce the angle 6 given by 


Pe 
So B may vary from 45° for D=0 to 0° for D—Dp. 

The quantity J, has the dimensions of a length in the same way as the 
quantity 7 introduced on the first page by formula 1. When, however, 
the crustal wave is formed below the sea-surface we have not to introduce 
6, in formula 2 but 6, — 1.03, i.e. diminished by the specific density of 
sea-water. If we indicate the length we thus get byl», we have the following 
relation between the three values 
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and introducing 0; — 6) —0.6 and 6, —3.27 we find 
lp== 06561 b==0,719 (eee a eee eee Ok) 


So the value of 1=60km which has been adopted for the anomalies of 
the map, corresponds to J, = 39.4 km and I, = 43.1 km. 

After introducing formula 5 in the equation 4, we find the general 
solution of 4 in the following shape 


y=Ae-%* cos(bx +) + A’e%* cos(bx+ 9’). . . (7A) 


in which A, A’, » and q’ are integration constants and 
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When D does not reach the buckling-limit, ie. when a > 0, this formula 
represents two waves of the same length of which one is increasing with x 
and the other diminishing, These waves are caused by other forces working 
on the crust besides the compression D and obviously the constants A 
and A’ must be chosen in such a way that the waves are dying out to 
both sides of the outward force which causes them, If there is no outer 
force working on the crust, A as well as A’ are zero and the crust does 
not undergo any deformation save of course the thickening by elastical 
compression; it remains plane. In our case of one vertical force working in 
the origin of the coordinates, we have the following set of formulas for y 
and its derivatives for positive values of x, representing a wave dying out 
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for greater values of x, and for negative x we find the same curve towards 
the side of the negative x. 
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The length L of the waves is 
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and as f varies from 45° for D—0 to 0° for D—D, the wave-length 
varies from 27|/21, = 8.90 1, to 2ml, —6.281,. Taking e.g. the case we 
shall treat of presently that D—0.6 Dy and assuming a submarine wave 
of the crust and so, according to 160 km, Il, 43.1 km, we find 
L = 303 km. This is in good keeping with the wavelength of the anomaly- 
map as has been stated in the previous §. 

For the determining of the integration-constant A we may use the value 
of the vertical force Py working on the crust in the origin. If Ho is the 
total shearing-force in a cross-section perpendicular to the crust for x 
infinitely small but positive, we have, if the index 9 indicates values for 
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Po =—26,g1, A sin (y — f) 
or 
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We may derive the second integration-constant y from the value of dy/dx 
for x =0. When the crust had not been broken in that point, we should 
have had to introduce dy/dx —0. When, however, we identify this point 
with the buckling zone, other values are also possible, and so we have in 
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general to introduce a certain given value for (dy/dx) . We may thus 


derive pm from 
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Introducing A and p in 8A we find the curve for y. | 
In fig. 3 three cases are given, all for the same value of Po, i.e. 
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and for each case three assumptions have been made for D. The drawn 
curve corresponds to no compression, ie. D0, the point-dot curve to 
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Fig. 3. Theoretical curves of crustal deformation under the effect of an upward 
force Po and a horizontal compressive force D. 


D=0.6 Do and the dotted curve to D—0.882 Dy. On purpose only 
positive values of (dy/dx) have been chosen, because as a consequence of 
the downward buckling in a previous period, positive values seem more 
likely than negative values. They appear also to follow from some of the 
profiles, e.g. those south of Java of which fig. 2b is an instance, The 
vertical scale of the curves has been chosen in such a way that for the 
following acceptable values of Po, 6, and 1, the vertical scale is ten times 
exaggerated with regard to the horizontal scale, i.e. both scales are the 
same as those of the topography in fig. 2a and 2b. For Py a value has been 
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adopted of 7.2 X 1012 gram per cm of the third dimension; this corresponds 
to a root of a cross-section of 1200 km2 and a deficiency of specific density 
of 0.6 and this again is in keeping with the negative anomalies south of 
Java. For 0, and 1, the values of 2.24 and 43.1 km have been adopted 
corresponding to a wave-formation below sea-level. Introducing this in 
formula 3 we find for the buckling-force Dy 8.3 X 1013 gram per cm of 
the third dimension, i.e. about 114 times as much as Pp, It is of course 
possible to introduce different values of (dy/dx), to both sides and also 
to derive equations for more than one vertical force working on the 
crust, but the writer must refer to a future publication for these more 
complicated solutions. 

We may derive interesting conclusions from fig. 3. In the first place 
we see that an increase of the compressive force D may bring about rising 
as well as sinking for the central zone, depending on the slope of the crust 
to both sides. Fig. 3b gives a case where this zone does not move at all. 
In all cases the waves to both sides increase in amplitude when D increases. 
When D=0, i.e. when there is no lateral compression, these waves are 
insignificant and they become larger and larger when D increases. When 
D is near the buckling-limit they become nearly as large as the one in the 
central zone. 

So we see that the effect of a change in the compression, even incase 
the force Py is kept constant, shows a great diversity, depending on 
circumstances, and this diversity is still more striking when we include 
cases where the slope of the crust to both sides of the central zone differs. 
This rather surprising diversity provides us with a means to adapt this 
solution to many different cases of behaviour of the crust and to bring 
them all back to one common origin. To point out these possibilities is 
one of the aims of the writer in this paper. In the next §' he will make use 
of these possibilities for a few cases in the East Indies but an exhaustive 
treatment is out of the question here; the writer must defer this to a 
future occasion. He wishes to lay stress on the fact that the basic assump-~- 
tion opening these possibilities is a simple and acceptable one, the supposi- 
tion of a compressive force working in the crust; the presence of the 
vertical force Py) is no assumption but an observed fact directly derived 
from the gravity anomalies. 

With a view to an hypothesis he wishes to make in the next § the 
writer has marked the points of maximum curvature in the downward 
waves of the curves of fig. 3. It is easy to derive their location from the 
equations 8 by equalling 8D to zero which gives 


et me eee (2) 


For D=0 these points of maximum curvature are farthest away from 
the deepest points of these downward waves; for increasing D they get 
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nearer to these points and they coincide with them for the pure sine-curve 
which is obtained for D equal to the buckling limit. 


§ 4. Geological interpretation of the equations of § 3. 

According to the summary given by J. H. F. UMBGROVE in “Geological 
History of the East Indies’ 1) the recent geological history of the islands 
of Timor, Tanimber, Kei and Ceram, all situated in the belt of strong 
negative anomalies, is the following. After the last period of strong 
folding and overthrusting in the upper miocene, a period of denudation of 
the elevated islands and subsidence followed so that in the pliocene the 
sea could invade them. Then trough faults and ‘‘graben” have been formed 
in the axis of the islands and about parallel to them, and finally towards 
the end of the pleistocene elevation occurred again and faulting. At the 
same time probably of this last elevation the troughs to both sides of these 
islands have sunk away. 

In the light of the results of the previous § we may perhaps try to 
explain these occurrences by changes in the compressive force D. A 
decrease of this force may have brought about the subsidence in the first 
part of the pliocene and this is in keeping with the fact that in no part 
of the Archipelago folding occurred in that period. The writer does not 
think that the compressive force can have disappeared entirely because 
we may surmise that in that case the belt of negative anomalies would have 
succeeded in restoring its isostatic equilibrium. As since that time there 
has not occurred enough folding in this zone to give an easy explanation 
for these great anomalies, we must probably keep to the idea that they 
have at least for the greatest part originated during the great folding period 
in the upper miocene or earlier, and this probably implies that in no 
period since then the compression has disappeared entirely. 

Towards the end of the pliocene or the beginning of the pleistocene 
folding begins to occur again in the Archipelago, e.g. in the oil-geo- 
synclines, in Boeton and in E. Celebes and this appears to point to an 
increase of the compressive force. There seems reason to bring this in 
connection with the formation of the “graben” on Timor and the other 
islands and to attribute this likewise to an increase of the compressive 
force; such a reaction of the old buckling zone does not seem difficult to 
explain. The great vertical movements towards the end of the pleistocene 
indicate a further increase of the compression, which according to the 
results of the previous §, must gradually have brought about the waves 
that in the present period cause the alternating belts of positive and 
negative anomalies that have been found. The rising of the islands of 
Timor, Ceram etc. and the sinking away of the troughs to both sides could 
thus be attributed to this wave-formation; the deformation there would, 


1) Bulletin of the American Association of Petroleum Geologists, Vol. 22, 1 Jan 
1938, 
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therefore, have to belong to the type of fig. 3a or between fig. 3a and 3b 

This explanation of the vertical movements makes it possible to under- 
stand that in other parts of the central zone slighter or even no rising has 
occurred, This may explain the great differences of elevation in this zone, 
the alternation of islands and lower parts between. South of Java we 
have a particularly low part of the zone and so we should have to classify 
this part under the types of fig. 3b or perhaps even 3c. We may find two 
remarkable confirmations of this view-point. In the first place the topo- 
graphical profiles, as e.g. that of fig. 2b, show a downward slope towards 
the central part and so a value of (dy/dx) leading to one of these types 
may be expected. As the horizontal and vertical scales of fig. 26 and 
fig. 3 are the same, we may directly compare the slopes in both figures. 
In the second place we find zones of positive anomalies at rather large 
distances of the central zone to both sides, of which that over the southern 
part of Java is one, and this is in better agreement with fig. 3c than with 
fig. 3a. 

On the other hand the profile of fig. 2a over the Tanimber Islands 
comes near to the type of fig. 3a and this is in keeping with this part of 
the central zone having risen to a higher position than that south of Java. 
For making this clear the writer has subtracted from the regional anomalies 
of part of the profile, as given by the drawn curve in fig. 2a, the negative 
anomalies caused by the root of crustal material and the resulting anomalies 
are given by the dotted curve of fig. 2a. These anomalies may thus be 
assumed to be caused by the wave-formation of the crust and the curve 
in fact does fairly well agree with the type of the curves of fig. 3a. The 
writer considers this dotted curve as a first and still rather uncertain trial 
to free the anomaly-curve from the effect of the root. As we do not know 
its exact dimensions and shape this effect may be different and so the 
remaining part of the anomalies would also be different; this would mainly 
affect the part over the central zone. In consequence of this it might be 
possible that a nearer consideration would put the type of the dotted 
curve between that of fig. 3a and fig. 3b. It could not, however, quite be 
of the type of fig. 36 or still further away from fig. 3a. 

In any case the inner Banda arc, represented in this profile by the island 
of Seroea, would coincide with a downward wave and thus it would cor- 
respond to the downward wave that in the profiles over Java and Sumatra, 
e.g. in the profile of fig. 2b, represents the oil-geosyncline. In view of the 
fact that both features occur in the same topographic axis, the axis of the 
greater and smaller Sunda Islands, this seems acceptable. The question 
arises whether we could thus understand the presence of strong volcanic 
activity in this axis. The writer is inclined to answer this question in the 
affirmative; a downward wave of the crust must bring about tension in 
the lower half of the crust or at least a diminishing of the pressure and 
this might explain the rising of magma. If this explanation is correct the 
location of the volcanoes ought more or less to coincide with the location 
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of the maximum curvature of the crustal wave and this has led the writer 
to indicate these spots in the curves of fig. 3. We see that they are not 
found in the axis of the downward wave but shifted towards the side of 
the central zone and more so for smaller values of the compressive force, 
i.e. for waves that die out quicker with the distance to the central zone. 
For the profiles over Java this seems to be more the case than for the 
profile over the Tanimber Islands and the location of the volcanoes to the 
south of the oil-geosyncline appears in good keeping with this surmise. 
Besides this view-point the writer still thinks that the supposition he has 
formerly 1) made regarding the occurring of volcanoes in belts on the 
inside of curved parts of the orogenetic belt, has a bearing on the subject, 
especially to explain that no volcanoes occur on the outside of a curve. 
This problem of the distribution of the volcanic activity in the Archipelago 
is, however, a complicated one, requiring much expert knowledge, and so 
the writer does not think he can do justice to it, certainly not in the few 
lines he can devote to it here. 

An important point has still shortly to be mentioned. In his paper on 
“The negative isostatic anomalies in the East Indies’ KUENEN mentions 
a problem brought forward by Mac GILLAvVRY as an objection to the 
buckling hypothesis, i.e. the question how it is possible that the mean 
elevation of the eastern half of the Archipelago is below sea-level and 
probably lower than it has been in former geological periods, although the 
sial crust has been compressed; it looks as if this ought on the contrary to 
have brought about a rising of the mean elevation. There is no doubt 
that if the profile of fig. 2a over the Tanimber Islands should have to be 
identified with the deformation of the crust as given by fig. 3a, we must 
assume that the undeformed position of the surface should have to lie at 
a depth of some 5000 m or more. And so it appears as if there has been 
a considerable sinking of this whole area in recent periods. It is a great 
problem to find a cause for such a phenomenon which the writer is not 
able to discuss here, but he may shortly mention the possibility that this 
sinking could be explained by assuming a great convection-current system 
in the substratum with its descending current under this part of the 
Archipelago and rising currents under Asia and Australia or under one of 
these continents. Because of the drag exerted on the crust by this current 
such a current could at the same time explain the compressive force in the 
Archipelago. 

The writer should not want to conclude this paper about the mechanical 
conception for explaining the gravity field and several geological features 
in the Archipelago, without mentioning also a serious difficulty for 
adopting it. The question is that the forces and stresses which must be 
assumed to work in the crust are considerable and larger than is generally 
supposed to be possible. If we adopt a value for D of 0.6 Dp and if we 


1) FF. A. VENING MEINESZ, Gravity Expeditions at Sea, Vol. II, pp, 112 and 125, 
Publication of the Neth. Geod. Comm. Delft, 1934. 
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take for Do the value derived in the previous §, we find the value of D to 
be 5 X 1013 gram per cm of the third dimension. Supposing the crust to 
have a thickness of 40 km we thus obtain a mean value for the compressive 
stress of 12500 kg/cm? and we get still larger figures for those parts of 
the crust where the bending causes additional compression. The writer 
cannot go deep here in the question at what figures we have to estimate 
these additional bending stresses and so he has to reserve that for a future 
occasion, but he wishes shortly to mention that they would become 
enormous if we kept to our simplified conception of the crust being one 
elastical plate. There is no doubt we have to come to a view-point nearer 
to the actual conditions in the crust for the attack on this question. In two 
ways we can understand that the bending stresses would be less. In the 
first place we may assume the crust to be composed of several layers which 
can more or less slide over each other and in the second place we can 
presume that the deformation for these large stresses, at least for parts of 
the crust, has not been of the elastical type but of the plastical or the 
block-faulting type. For taking up these problems we have to write the 
equations 1 and 2 in a more general shape by substituting 
Se ea, 
ee 12 (m2 — 1) 


in which f is in general the ratio of the moment of force M working in the 


da 
crust to the value re of the curvature. We thus get 
aC 


— p-8Y Ge fa 
Mee ee (132A) '=|/ ga ica) 


The equations 6 for the values of 1, and /, remain the same and also all 
the other equations and formulas. 

The writer cannot go further in these problems here and so he wishes 
to conclude his paper with the general remark that this attempt for giving 
a mechanical explanation seems successful for many questions which are 
otherwise difficult to explain or to bring under one common view-point, 
but that the stresses and strains implied are large. He thinks, however, 
that it will not be easy to escape such a conclusion. 

The writer wishes to express his thankful acknowledgements to his 
colleagues UMBGROVE, KUENEN, ESCHER and RUTTEN for their helpful 
discussions; many view-points mentioned here may in fact have originated 
in these discussions, especially in the long talks he had about these matters 
with the former two. He wants further to acknowledge the assistance of 
Mr. J. vAN Dijk of the geological Institute of Utrecht for the making of 
the map and drawings. Lastly he wishes to recall that the whole gravity 
material is due to the great cooperation of the Netherlands Navy and to 
the helpful assistance of the Captains MANTE and VAN DER Kun, of the 
Officers and of the men of the submarines on board of which the observa- 
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Physics. — La pression du toit sur le charbon pres du front dans les 
exploitations par tailles chassantes. Par F. K. TH. VAN ITERSON. (II). 


DEUXIEME CHAPITRE. 


Répartition des pressions aprés écrasement du charbon. 


(Communicated at the meeting of January 27, 1940.) 


§ 4. La pression du toit sur les piliers de charbon. Il sera tenu compte 
des propriétés différentes de charbon et schiste, les deux considérés comme 
matiéres élastiques. 


Dans le paragraphe précédent nous avons considéré le rocher, dans le 
sein duquel des galeries sont creusées, comme un milieu isotrope, homogéne 
et élastique. Nous le supposons chargé en dessous de la limite d’écoulement 
et acceptons que la loi de HOOKE est applicable. 

Mais les propriétés élastiques du charbon et du schiste sont tellement 
différentes que la répartition des pressions autour d'une galerie creusée 
dans la veine de charbon sur toute son épaisseur, ne peut pas étre traitée 
comme nous le faisions pour les galeries en terrain homogéne. Pas plus 
les pressions dans le toit et le mur de schiste au-dessus et au-dessous d'un 
pilier de charbon, ne se répartissent de la méme maniére qu’autour des 
piliers laissés en terrains homogénes. 

Pour prouver cette thése, il faut maintenant indiquer en quoi ces matiéres 
se distinguent. Plus tard nous parlons de la nature cassante du schiste 
partout dans la mine ot il est mis 4 nu et de sa plasticité quand il est caché 
a notre vue et soumis dans tous les sens a une pression élevée mais inégale 
selon les axes principaux. Nous ferons mention de la fragilité du charbon, 
gui aprés rupture se comporte comme une matiére incohérente, mais ici 
nous donnons des résultats d’essais sur la résistance 4 la pression, le 


coéfficient d’élasticité et de contraction latérale des deux matiéres en 
question. 


LEON Morin, Quelques effets de pression de terrains. Bulletin de la Soc. de l'industrie 
Minérale 1912, p. 29, 
L’écrasement du charbon se produit 4 40 kg par cm? 
D Gls Senigas nancies 2 BO 4 4 
i: Hi an moyens a 170 ,, 
5, ¥, +3 durs a 220 
HEISE—HERBST, Bergbaukunde 1932, 
Résistance a la compression: 
charbon 35— 290 kg par cm? 
schiste 250— 850 ,, 
greés 290—1840_,, 


” ” 


9 ” ” 


” ” 


” 
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MULLER, Gliickauf 1930: 
schiste — 560 1 815 kg par cm? 


grées — 7704 800 oie 
PHILIPS, Annales des Mines de Belgique. Les roches houilléres 1938, i, Bless 
4 normal a la couche — paralléle a la couche. 
charbon dur anglais 208— 243 kg par cm? 
schiste Als 400— 800 ae 
= 260— 420 
grés zB 700—1400 ,, » » 
24 = 450—1300 ,, 


Ne ferait-on pas mieux de renoncer A |’élaboration d'une théorie apte 
a réconcilier la mécanique appliquée avec l'art de exploitation de mines 
quand on voit une telle divergence dans les propriétés des matériaux? Les 
difficultés s’entassent quand on se rend compte que la stratéfication peut 
étre irréguliére, les couches varient d'argile plastique a cuerelles dures. 
Le plissement des couches a changé le pendage, provoqué des dérange- 
ments, des failles, ouvert des clivages dans le charbon, transformé le rocher 
en diaclases et en sus de toutes ces irrégularités géologiques il y a les 
ruptures, le foissonnement, les dislocations occasionées par les travaux en 
tailles dont l'influence se fait souvent sentir jusqu’a cinquante métres des 
havées et davantage. 

Pour démontrer la complexité du probléme nous donnons quelques 
chiffres suggérant une connexion entre les dimensions des cubes d’essais 
et la résistance a la compression du charbon, 


Résistance a la compression 


P eae Dimensions des Reo cence en 
Ase! | ES oe Et cin cubes d’essai en cm kg par cm? 
a Sea ] 
Rice-Ensiau Anthracite de | 6.4—25.4 175 (6 fois) 
Bull. Bur. of Mines | Pennsylvanie 25, 4—-30)5 | 141 (3° fois) 
1929, p. 303 | me 75 57.4 (1 fois) 
~ 140 21.6 (1 fois) 
Lawall-Holland Virginia 7.6 100—300 
Trans, A. I. M.E. 
1932 
MULLER en Gliick- | Haute Silésie A ep) 210 
auf 1930 | 415) 165 
Bureau of Mines Charbon dur 125—400 
selon MULLER en Charbon tendre 30—150 
Gliickauf 1930 


Tous ces essais sont exécutés sur des cubes entre surfaces frottantes, 
et par conséquent sont erronés. FOppL1) et mieux encore Ros et 


1) FOpPPL, Technische Mechanik, III, p. 71. 
20* 
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EICHINGER !) ont démontré qu'il faut diminuer ces chiffres au moins d'un 
tiers pour avoir la résistance a l’écrasement de la couche. 

Chaque mineur connait au charbon en site deux espéces de plans de 
clivage, les plus manifestes qui s entrecroisent réguliérement sont dus 
a la poussée géologique et traversent la veine sous une inclinaison déter- 
minée et les plans de clivage paralléles au front de taille qui sont dus a 
Vexploitation de la couche. On reconnait que les premiers sont d'origine 
tectonique parcequ’on les trouve quelquefois remplis de calcite ou autres 
minérais, et les derniers se produisent sous nos yeux avec un crépitement 
bien connu, Il nous parait probable que la divergence entre les résultats 
d’essai comme la dépendance de la résistance a la compression des dimen- 
sions du bloc, sont causées par la fissuration des morceaux de charbon. 

Quand on ajoute a ces variations de la substance du terrain, les disparités 
de la pression dues aux poussées tectoniques et plus encore aux affaisse- 
ments provoqués par l’abatage, on désespére pouvoir exprimer en formules 
la distribution de la pression autour de la taille. Mais les régularités avec 
lesquelles les effets de pression se manifestent dans les exploitations nous 
incitent a continuer l'étude. 

Que sait-on des modules d’élasticité E et de la contraction latérale m? 
On trouve des chiffres chez LowEns, Gliickauf 1933, p. 34-38, MULLER, 
Gliickauf 1930, p. i601, Lawall-Holland Trans. A.I.M.M.E. 1932 et 
surtout chez Puitips, Annales des Mines de Belgique 1938, p. 564. Les 
résultats des essais sont variables mais pas tellement comme la résistance 
a la compression. Nous donnons ici quelques chiffres qu'on peut accepter 
comme moyennes, 

Pour le grés et le schiste E et m sont des fonctions de la pression. On 
peut accepter 

FE = 600.000 kg/cm? aux faibles pressions 

E = 100.000 kg/cm? aux pressions d’écrasement 

m diminuant en méme temps de 16 a 3,4. 

Aussi 4 quelque distance des exploitations ou le rocher est soumis a une 
pression dans tous les sens m ~ 3,5. 

En contraste avec le rocher pour le charbon dans la couche ces modules 
sont des constantes, on peut accepter 


E = 50.000 kg/cem2 m=2. 


Seulement pour les charges trés faibles on trouve pour m un peu plus 
au laboratoire, mais en chargeant davantage m tombe rapidement a 2 et 
reste constant. 

La constatation que m==2, ce qui veut dire que le charbon ne varie 
gue fort peu en volume quand on le comprime, qu’il se comporte comme le 
caoutchouc dans ce respect, tandis que le rocher au contraire ne se dilate 


1) Ros et EICHINGER, Versuche zur Klarung der Frage der Bruchgefahr. Diskussions- 
bericht No. 28. Nichtmetallische Stoffe, Ziirich 1928, p. 14. 


presque pas latéralement quand 
surprenantes pour l’abatage 1), 
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il est comprimé, méne a des conséquences 


ee Mates Module d’Elasticité Module de Decon 
ae FE en kg par cm? m 
LOWENS, Gliickauf Grés 1 450000 = 400000 3,4 
IGE), ja, Shee | Schiste sableux 4 275000 = 329000 
Schiste argileux t 547000 = 637000 
MULLER, Gliickauf Grés 1 380000 = 510000 
1930, p. 1601 Schiste 600000 
Schiste argileux 744000 
GREENWALD, Journ. | Grés de Pittsburgh ~ 200000 
appl. Physics 1937 Schiste de idem ~~ 200000 
S 2 
PHILIPS eres 100000 — 600000 | $_ 1° 2 72 ka/em 
( 7,64 a 430 
Ann. Mines de : 
; ( 50) & 710) * 
Belgique 1938 | Schiste 133000 — 400000 (34 a 840 
| DE) er 7h0) % 
Charben 50000 — 58000 | $ “> * 
| (2,03 a 430, 
Lawall-Holland | Charbon de 1 26000 — 43000 O59 
Trans. A.I.M.M.E. | Virginia = 2000 = 50000 = 3,28 
1932 | 
| 
MULLER, Glickauf | Charbon de 
1930 2) | Silésie 24000 — 55000 


1 veut dire perpendiculaire a la couche. 
= veut dire pression paralléle a la couche. 


Référons-nous 4 la figure 1. La nous avions supposé le charbon comme 
indéformable, nous y reviendrons, mais pour le moment nous constatons 


1) D. W. PHILIPS and A. HACK, The physical properties of coal measure rocks. 
Seventeenth Annual Report of the Safety in Mines Research Board 1939, p. 66. 

2) Pour le module d’élasticité de la roche houllieuse on trouve encore quelques 
chiffres dans: 

O, FLEISCHER, Bergmannische Betrachtungen und Versuche tiber Gebirgsbewegungen 
im oberschlesischen Steinkohlenbergbau zur Klarung von Gebirgsdruckfragen. Diss, T. H. 
Breslau 1933, 

O, MULLER, Untersuchungen an Karbongesteinen zur Klarung von Gebirgsdruckfragen. 
Diss. T. H. Breslau 1930, MULLER donne ici comme moyen £ pour le charbon 25000 kg 
par cm?, 

K. STOCKE, H. HERZMANN und H. UDLUFT, Gebirgsdruck und Plattenstatistik. Elasti- 
zitatsverzuche an karbonischen Gesteinen Oberschlesiens, Z. Berg-, Hiitt.- u. Sal. Wes. 
82, 307—354 (1934). & pour le schiste en moyen 450.000 kg par cm?. F pour le grés en 
moyen 200.000 kg par cm?. 
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que les bases de notre calcul viennent a nous manquer. Nous avions seule- 
ment calculé avec l'élasticité du rocher, qui est minime. Pour la couche de 
charbon au contraire la détente est grande en vertu du petit module 
d’élasticité et provoque une grande expansion vers le vide et en meme 
temps la veine se contracte beaucoup en épaisseur. Le déplacement du 
front fait glisser le charbon par les surfaces du rocher et le détache. Mais 
que signifie m= 2? 

Mettez un cube ou un cylindre de charbon dans un moule indéformable 
qui l’englobe parfaitement et mettez la pression p sur la surface libre, alors 
cette pression régne dans tous les sens dans le charbon et contre les parois 
du moule, Otez la pression et les parois sont dégagées. On est tenté a 
croire qu'il en est de méme dans la mine. Quand on met a nu la surface 
du charbon on dégage toit et mur de la pression sur une certaine étendue. 

Comme le disait Lord KELVIN: ,,On ne connait un phénoméne que quand 
il est possible de l'exprimer en formule’. Nous tacherons de traiter le 
probléme selon les lois de |’élasticité. Pour ne pas trop compliquer cette 
tentative nous négligeons le frottement du charbon le long du toit et du 
mur. L'observation de tous les jours nous apprend que ces surfaces se 
présentent comme trés lisses en général, le charbon y glisse sur une dizaine 
de cm, souvent les surfaces ont l’air d’étre graissées. Parfois méme une 
couche mince limoneuse lubrifie ces surfaces. Nous avons déja parlé des 
plans de clivage paralléles a la direction de la pression dus a l'abatage et 
pour simplifier la représentation mentale nous considérons la couche de 
charbon comme étant fendue en d'innombrables colonnes juxtaposées. 
Toujours nous traitons le probléme a deux dimensions y et z. 


$= -£ _g:+8 


Zz 


Fig. 11. Notations pour calculer l'augmentation de la pression vers les bords 
sur un pilier, sur une bande de charbon inexploitée. 


Par raisonnement nous croyons que la pression sur le pilier de charbon, 
laissé entre deux vides (fig. 11), augmente vers les bords. Le cas extréme 
étant celui de la bande indéformable, ou bien la bande extreémement mince 
ou trés large, traité dans le chapitre I (p. 99, 1938, figure 10) oa les 
pressions tendent vers l'infini quand on se rapproche des bords. 


I] parait impossible de trouver la solution exacte pour la répartition de 


BSS) 


la pression sur une bande élastique comprimée entre deux hémi-espaces 

de matiére a autre module d’élasticité:; mais on peut chercher une solution 

serrant entre limites la répartition réelle et alors approcher ces limites. 
Nous essayons dans ce but p en fonction de y 


2 4 6 
P= Po Bo thee mis 


Calculons la déformation de la surface de charbon et celle du mur et nous 
pourrons choisir po, k, 1, m, etc. de telle maniére que les surfaces s'épousent 
en autant d’endroits qu’on a d'inconnus, le résultat sera que les deux 
courbes se confondent. Déja avec trois inconnus Po la pression au milieu et 
deux coéfficients k et / on obtient un résultat satisfaisant a notre but. 

Dans l'annexe 1 nous donnons quelques détails du calcul pour trois 
cas: b—a, b= 2a et b= 4a c.-a-d. largeur du pilier de charbon une, deux 
et quatre fois sa hauteur et nous discutons a présent les résultats repré- 
sentés dans les figures 12, 13 et 14. 


Po P 


a + 


~— 


Fig. 12. Répartition de la pression du toit et du mur sur une bande de charbon 
de largeur 2b — 2a, E rocher = 360.000 kg par cm?. FE charbon — 60000 kg 
par cm?, 
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Fig. 13. Répartition pour une bande de charbon de largeur 26 — 4a. 
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Fig. 14. Répartition pour une bande de Charbon de largeur 2b = 8a, 
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Dans chaque figure la courbe inférieure donne la déformation de la 
surface du charbon 4a échelle exagérée pour les coéfficients k et | adoptés. 
Ces déformations, ces raccourcissements sont en méme temps une mesure, 
représentent aussi la distribution de la pression. La courbe supérieure 
représente l'impression du mur (le sol) sous cette distribution de pression. 


On remarque: 


1. Au bord méme la tangente de la courbe inférieure est peu inclinée 
et celle de la courbe supérieure prend une inclinaison de 90°. Notre solution 
ne donne donc pas la répartition mathématiquement exacte des pressions 
au ras du front de charbon. Mais l'influence locale des coins s’amortit a 
peu de distance. 


2. La répartition exacte de p selon la largeur du pilier de charbon est 
située entre nos deux courbes, qu'on peut serrer en ajoutant plus de termes 
a la formule. 


3. Pour les piliers relativement étroits que nous avons examinés, avec 
b une, deux ou quatre fois a, l’'augmentation de la pression vers les bords 
reste limitée, mais elle s’accroit avec la largeur du pilier, ou ce qui revient 
a la méme courbe pour la distribution de pression, l’augmentation vers les 
bords s’accroit avec amincissement de la couche, jusqu’a devenir infinement 
grande pour les piliers trés larges ou trés minces, sauf au ras méme du 
charbon ot la solution ne vaut plus sur une étendue trés restreinte. 


4. Nous résumons nos calculs dans le tableau suivant. 


| Augmentation de la pres- 


ee sion et de la compression 
Formule approximative 


Largeur 


2b 


comparée 


a 


l'épaisseur 2a du pilier de 


pour calculer l'augmen- 


du pilier aux bords com- 


| parée a celles du milieu. 
| Le rapport réel est entre 
| 
| 
| 


ors tation de la pression 
ar vers les bords. ces deux chiffres 
(py — Po): P 

ye me 

nee 0,06 po & ze i) 0,13 et 0,12 
cy? y" 

pie=wa 05115 po (i a - 0,252 eethO.23 
Te tee 

b=4a 0,22 po (Fi a 5) 0,49 et 0,44 

b= ce 0 0 


5. Le calcul démontre que les pressions augmentent vers le ras du front 
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de charbon, mais pour les piliers de largeur restreinte dans une mesure 
moderée. Il n'est pas dit qu’aux faibles profondeurs la résistance a la 
compression soit surpassée, 

Aux Etats Unis, ott le charbon est plus dur et moins fissuré que chez 
nous l’exploitation par chambres et piliers est trés en vogue. Le charbon 
n étant pas écrasé par le resserrement entre toit et mur doit étre scié par 
des machines haveuses et détaché a coups de tir. 

Notre calcul donne une idée de ce que serait la distribution de la 
pression si l’on compare le pilier de charbon avec un tampon de caoutchouc 
séparant toit et mur. Nous avons considéré les surfaces de glissement entre 
le charbon qui se dilate beaucoup et le toit et le mur qui n’expendent pas 
latéralement, comme parfaitement lisses et la seule remarque qu’on doit 
faire est qu'ainsi nous avons négligé le frottement. Plus loin nous verrons 
qu'un peu de frottement change complétement la répartition de la pression, 
dans ce sens, que la charge au ras du front diminue et augmente vers 
l'intérieur du pilier. 


§ 5. La pression du toit sur la couche de charbon, cette matiére con- 
sidérée comme parfaitement plastique. 


Quand on exploite la veine de charbon, quand on met a nu le front et 
le soulage de la pression, le charbon se détend vers le vide. Le volume de 
charbon resserré entre toit et mur qui s’approchent, ne change pas. Le front 
se déplace. Nous verrons que le frottement qui s’oppose a l’avancement 
de la masse, excite une pression latérale, qui s’accroit vers l'intérieur 
jusqu’a ce que le charbon soumis a la compression dans tous les sens, et a 
raison de m— 2, de volume constant, résiste a la compression comme de 
la matiére absolument dure. 

Pour exprimer le probléme posé en formules et en nombres nous I'abor- 
dons par la théorie de la plasticité, Nous considérons le charbon comme de 
la matiére parfaitement plastique, se comportant de méme que l'argile, 
le plomb, le fer chauffé au rouge, et comprimé entre des surfaces paralléles 
ou les tensions se disposent en vertu des lois de I’elasticité, Les calculs sont 
donnés dans l’annexe 3. Ils nous permettent d’évaluer en chiffres les 
tensions dans tout le massif du toit et du mur. 

Pour un cas spécial nous calculons la distribution des tensions principales 
selon l’axe des Y et l'axe des Z, 

Profondeur au dessous de la surface h = 800 m. 

Poids spécifique des terrains pee) 5, 

Epaisseur de la veine de charbon 2a — 1 m. 

Largeur du vide créé dans le charbon 2b — 4 m. 

Tension maximum de cisaillement (condition de plasticité) 
ps = 10 kg/cm2. 

Alors on obtient c—b—10 m. 


1p Waele: 
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Le charbon se déforme, coule, plastiquement sur une distance de 10m 
de chaque coté du creux (figure 15). 


, ae 


Fig. 15. La matiére se comporte comme plastique de 26 jusqu’a 2c et en dehors 
de cette limite se comporte comme élastique. Nous pensons a la flexion légére 
du toit et du mur. 


Les tensions (toutes des pressions) sont représentées a échelle dans la 
figure 16 pour tous les points le long de l’axe des Y et de l’axe des Z. 


| 
| 


axe aes y 


compression 100k 9/ cm? 


axe des z 


Fig. 16. Distribution des tensions normales oy et oz et de la tension de 

cisaillement t le long de l’axe des Y et des tensions principales Py et Pz le long 

de l'axe des Z pour de la matiére plastique resserrée entre toit et mur de 
matiére dure. 
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On remarque une inadmissibilité dans le cours des o, selon l’axe des Y. 

Ces tensions ne peuvent devenir infiniment grandes pour y=—b et 
y—c. Cette absurdité est due a une acceptation qu'on a du faire pour 
pouvoir intégrer les équations des tensions. L’influence de cette acceptation 
est discutée dans l’annexe 3. En vérité l’anomalie de o, dans les points 
y—b et y=c nexiste pas. On peut raccorder les trois parties de la 
courbe par des arrondissements légers. 

La base de tout le calcul est que la tension de cisaillement tmax. = ps 
est constante (nous prenions p,=10 kg/cm?) sur toute l’étendue des 
surfaces de contact entre toit et mur. Nous supposons qu'il existe la plus 
grande résistance possible contre le glissement de la masse plastique, 
considérons ces surfaces entaillées comme celles d’une lime. 

I] résulte de ce calcul, de la représentation des résultats en la figure 16 
gue quand le toit et le mur resserrent une masse plastique ils sont épargnés. 
Au bord du front les tensions sont presque nulles et les tensions (pressions) 
élevées ne se développent qu’a une certaine distance du vide, la ou elles 
sallient avec la pression dans tous les sens de maniére a devenir in- 
offensives. 

Si le toit et le mur étaient séparés par de la matiére vraiment plastique, 
on naurait besoin d’aucun souténement, du moins si le rocher n'est pas 
trop fissuré et disloqué par les poussées tectoniques. 


§ 6. La pression du toit et du mur sur la couche de charbon, cette 
matiere considérée comme pulvérulente. 


Mais le charbon n'est pas une matiére parfaitement plastique. Il est au 
contraire trés fragile et prés du front il est fissuré par les plans de clivage. 
Le charbon ne devient plastique qu’a des pressions de quelques milliers 
d’atmosphéres. C’est seulement a ces pressions qu'on peut en fabriquer 
des briquettes sans agglutinant. A 3000 atmosphéres méme on n’‘obtient 
qu'une quasi cohérence, les briquettes ainsi obtenues se désagrégent quand 
on les met dans I’alcool. 

A la pression rencontrée dans nos exploitations, au plus quelques 
centaines d’atmosphéres, le charbon écrasé, méme disloqué en gros mor- 
ceaux se comporte comme de la matiére incohérente. Mieux que d'appliquer 
les lois de la plasticité on fera appel a la théorie de la poussée des terres, 
on appligquera les lois de l'équilibre dans les massifs a frottement interne, 
ce gue nous faisons dans l’annexe 4. 

Si l’on remplace le monolithe de charbon par du menu charbon, par du 
sable, il est possible de calculer par des méthodes connues, la distribution 
des pressions dans le rocher et dans la masse pulvérulente, quand on 
comprime celle-ci, jusqu’a ce que le front s'avance. On trouve alors que la 
matiére sécoule d’elle-méme si l’on n’exerce pas une certaine pression 
sur le front. Mais si on prend cette précaution, si l'on soutient le charbon 
par du boisage on trouve que le frottement le long du mur et du toit 
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fait accroitre la pression vers I'intérieur selon une loi exponentielle, Pour 
l'exemple calculé dans l'annexe 4 cet accroissement est représenté dans 
la figure 17. Nous avons construit notre figure avec une pression po 
minime et avec un petit coéfficient de frottement. En réalité. l’accroisse- 
ment de la pression, quand on s’éloigne du front et entre dans la couche, 
est tellement impétueuse qu'un peu de pression contre le front, qu'un peu 
de cohérance dans la matiére suffit a faire s’élever la pression du toit et 
du mur trés rapidement a des valeurs extremes. 

Ainsi le frottement est capable de réaliser la distribution des pressions 
dans le rocher, représenté dans la figure 1, sauf la pression infiniment 
grande au ras du front du charbon,. Celle-ci sera beaucoup atténuée. 

Mais quand la pression du toit et du mur atteignent des valeurs tellement 
élevées, quand grace au frottement il se développe dans le toit et le mur 
enserrant la matiére incohérente une distribution de pressions comparable 


Fig. 17. Accroissement de la pression du toit et du mur sur de la matiére 
incohérente selon une loi exponentielle et notations pour le calcul dans l’annexe 4. 


a celle qui se produit sur une couche de matiére parfaitement dure et 
indéformable, quand les pressions et surtout les tensions de cisaillement 
dans le toit et le mur surpassent de beaucoup la résistance a la rupture, 
il n’est plus permis de considérer le rocher comme incassable. 

Nos calculs, basés sur les lois de l'élasticité, perdent alors leur validité. 

Dans nos exploitations on observe des déplacements du front de charbon, 
fissuré, clivé, brisé en gros morceaux. Les géométres des Mines de ]'Etat 
ont fait maintes observations et relevés de ces déplacements, tant des 
avancements du front de charbon relatifs au toit et au mur que des trans- 
lations absolues. Les déplacements relatifs, irréguliers sont de l’ordre d'une 
dizaine de cm pour une couche d’un métre d’épaisseur. Dans la litérature on 
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trouve un grand nombre d’observations qui confirment nos constatations !). 

Mais le charbon en site n'est pas précisement de la matiére incohérente. 
Aussi quand on l'abat, le front n'est pas muni d'un boisage. Le fait qu'il 
ne s'écoule pas de lui méme, qu'il faut l’abattre au marteau-pic pneumatique 
gqu’on introduit dans les plans de clivage, démontre que nous n’avons pas 
encore la loi exacte qui commande les déplacements a J'intérieur de la 
couche de charbon exprimé vers le vide. En effet toute la masse n'est pas 
pulvérulente. I] y a de gros morceaux qui se réfractent. Pour les matiéres 
plastiques la résistance a la traction et a la compression sont égales, mais 
le charbon ne supporte que trés peu de traction et la résistance a la 
compression est grande quand on ajoute de la pression dans tous les sens. 
Dans l'annexe 4 nous appliquons 4a la fin la loi, la condition de rupture 
pour cette espéce de matiére et trouvons de nouveau la distribution, 
l'accroissement rapide de la pression, selon la formule exponentielle. Ce 
que nous avons dit de l'influence du frottement qui fait croitre la pression 
prés du front outre mesure, s'applique a plus forte raison pour le charbon 
encore cohérent, sauf que le front n’a méme pas besoin de souténement. 

Par ce long chemin nous sommes arrivés 4 la conclusion qu’aux profon- 
deurs ou nous exploitons le charbon l’accroissement des pressions prés du 
front est tel que le toit et le mur ne peuvent pas rester indemnes. Il s'y 
produit au moins des crevasses, mais si l'on examine le toit avec attention 
on trouve des dérangements de blocs qui certainement sont dus a l'abatage 
du charbon. 

Jusqu’a présent nous avons toujours supposé que le rocher demeurait 
intact, mais de cette maniére nous ne pouvions pas expliquer le comporte- 
ment du chantier pendant l’exploitation. Nous espérons traiter dans un 
troisiéme chapitre la distribution des tensions dans le rocher et dans la 
veine autour de la taille quand on aura dépassé la résistance des deux 
matiéres a la rupture. 


1) WEISSNER, Gebirgsbewegungen beim Abbau flachgelagerter Steinkohlenflésse, 
Gliickauf 22 Okt. 1932, p. 945. 

LOFFLER, Zur Abbaudynamik bei streichendem Blindortbetrieb, Der Bergbau 9 Juni 
1938. 


Hydrodynamics. 


On the application of viscosity data to the determina- 
tion of the shape of protein molecules in solution. By J. M. BurGErs. 
(Mededeeling N°. 38 uit het Laboratorium voor Aero- en Hydro- 
dynamica der Technische Hoogeschool te Delft.) 


(Communicated at the meeting of February 24, 1940.) 


1, In Ch. III of the ‘Second Report on Viscosity and Plasticity” 
formulae have been given for the resistance of small particles of elongated 
form, and for their influence upon the effective viscosity of the liquid in 
which they are suspended; and the application has been discussed of these 
formulae to the results obtained with suspensions of polystyrenes by 
STAUDINGER and SIGNER!). A discussion of their application to suspensions 
of methyl cellulose has been given by SIGNER and v. TAVEL 2). PoLson 
has applied the formula for the influence of such particles upon the 
effective viscosity to the analysis of data obtained with suspensions of 
proteins 3), and a report of this work recently has been given by PEDERSEN 
in Part I, Ch. B, of SVEDBERG and PEDERSEN’s new book ‘‘The Ultra- 
centrifuge’ 4). 

The way in which the formulae for the viscosity and for the frictional 
coefficient have been applied by POLSON and the discrepancy which has 
been found between certain calculated and observed values, calls for some 
remarks which will be collected in the following pages. In connection with 
these remarks a few data also will be supplied for some cases not treated 
in the “Second Report”, viz. for disk-shaped particles (oblate rotational 
ellipsoids), and for certain systems consisting of rigidly connected spheres. 


1) ji. M. BURGERS, Ch. III of the “Second Report on Viscosity and Plasticity”, 
Verhand. Kon. Nederl. Akad. v. Wetenschappen (le sectie) Vol. 16, No. 4 (Amsterdam 
1938), pp. 122—126 (resistance formulae), 145—153 (influence upon the effective 
viscosity), 168—181 (application to polystyrenes). 

2) R. SIGNER und P. Vv. TAVEL, Die Form und Grésse von Methylcellulose-Molekeln 
in Lésung, Helv. Chim. Acta 21, 535 (1938). The subject also has been treated in the 
chapter contributed by R. SIGNER to TH. SVEDBERG and KAI O. PEDERSEN’s book ‘The 
Ultracentrifuge” (Oxford 1940), pp. 431—442, whereas cellulose acetates and some other 
linear high polymers are considered by E. O. KRAEMER and J. B. NICHOLS, ibid. 
pp. 416—431. — See also footnote 20) below. 

3) A, POLSON, Ueber die Berechnung der Gestalt von Proteinmolekiilen, Kolloid- 
Zeitschr. 88, 51 (1939). — See also: Nature 137, 740 (1936). 

4) TH, SVEDBERG and KAI O. PEDERSEN, The Ultracentrifuge (Oxford 1940), 
pp. 38—44. The results of the calculations also have been given by TH. SVEDBERG in 
the “Opening Address to a discussion on the protein molecule’, Proc. Roy. Soc. (London) 
B 127, 9—10 (1939). — It may be mentioned that the problem of the determination of 
the shape of tobacco mosaic virus particles in solution has been discussed by J. R. 
ROBINSON, Nature 143, 923 (1939). 
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The procedure followed by POLSON in order to obtain the molecular 
weight of a protein exclusively from viscosity and diffusion data, without 
having recourse to sedimentation measurements, was as follows: The 
length-diameter ratio L/d of the molecules, which were assumed to have 
the form of elongated rotational ellipsoids, was deduced from the specific 
increase of the viscosity of a solution, 75) —= 1/9 —1, making use of the 
known value of the partial specific volume V of the dissolved protein. The 
result was applied to derive the absolute dimensions of the molecule from 
its experimentally determined frictional constant fp»), which can be 
obtained from the diffusion constant D (measured by means of an optical 
method 5)) through the equation f/x, = RT/D. From the dimensions and 
V the molecular weight then can be calculated. 

However, on comparing the molecular weight found in this way with 
that deduced from observations on the sedimentation equilibrium, as is 
done e.g. in table 5 of “The Ultracentrifuge’ ®), it appears that the 
calculated values are much too low: by about 30 % when KUHN’s formula 
for 7sp7) is used, and by about 50 % when the more exact formula, derived 
from JEFFERY’s calculations by BURGERS 8), is applied. A correct result 
could be obtained only when instead of these theoretical formulae, an 
empirical expression, given by POLSON 9), is taken. 


2. It would appear to the present writer that a more convenient basis 
for a discussion of the cause of the discrepancy can be obtained by 
arranging the calculations in a different way. 

As explained in the “Second Report on Viscosity and Plasticity’’ 10), 
molecular weights always should be deduced from direct measurements, 
in the present case preferably from the sedimentation equilibrium (M,), 
or else from the sedimentation velocity in combination with the value of 
fexp as obtained from diffusion measurements (M,) 11). These values of 


°) See A. POLSON, Nature 137, 740 (1936), where it is stated that the diffusion 
constants were measured by the refractometric method of O. LAMM, Zeitschr. f. physik. 
Chem, A 138, 313 (1928) and B 143, 177 (1929). 

6) “The Ultracentrifuge”, Table 5, p. 44, — The same table occurs in TH. SVEDBERG, 
Proc. Roy. Soc. (London) B127, 10 (1939) and in A. POLSON, Kolloid-Zeitschr. 88, 
59 (1939). : 

") W. KUHN, Zeitschr. f. physik. Chemie A 161, 24 (1932). 

8) Second Report, pp. 152—153. 

®) This formula can be written in our notation as follows: 


Nsp/c = V [4,0 + 0,098 (L/d)?], 
as the quantity G used by POLSON is equal to cV. See A. POLSON, Kolloid-Zeitschr. 
88, 56 (1939); K. O. PEDERSEN, The Ultracentrifuge, p. 43. 
10) Second Report, p. 184. 
1) TH. SVEDBERG, The Ultracentrifuge, p. 9. — The fact that the two methods 
in those cases where both can be applied, practically lead to the same values for the 


molecular weight, proves that diffusion and sedimentation velocity are both governed 
by the same mean frictional coefficient, 
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the molecular weights will be assumed as trustworthy here. When now 
the following assumptions are made: (1) that the volume taken in by the 
molecule in the solution can be deduced from M by means of the value 
measured for V; (2) that the molecule has the form of an elongated 
rotational ellipsoid, the value of 7, can be used to find the value of L/d 
(under the circumstances of the viscosity experiment). The formula to be 
applied is: 


Heol C = V AGiacs Sree fe 2 ed) 


where 7s,/c is obtained from PoLson’s data for d(y/7,)/dn 12); V from 
data given in table 48, p. 406, of “The Ultracentrifuge’; while A, is a 
function of L/d, tabulated in the ‘Second Report” 13), 

Making use of the equation: 


ran Om VIR VANN Ges oe 6 ie Ga te (2) 


(where Ny is AVOGADRO’s number — 6,06. 10 2%), it then becomes possible 
to calculate the dimensions L and d of the molecule. These now can be 
used to find the frictional constants, f, for the motion in the direction of 
the axis of the molecule, and f. for the motion in the direction perpen- 
dicular to the axis; from these the mean frictional constant is derived by 
the equation: 1/f ,—1/(3f;,) + 2/(3f.). When we write: 


pour WHEN Ala” Ga) Pee eee Be, 3) 


then for elongated rotational ellipsoids 14): 


jim log (Lid + VV L?/d? —1) 


VAS eee (4) 


(the log being the Napierian logarithm). Finally the sedimentation con- 
stant 15) S = (dx/dt)/w2x is given by: 


vik oye al 


a= Bir UN ae (5) 


When the calculations are executed in this way for the proteins men- 
tioned by PoLson and PEDERSEN, and the results are compared with the 
experimental values of the sedimentation constant, it appears that the 


12) A. POLSON, Kolloid-Zeitschr. 88, 58, Tab. VI (1939). 

13) Second Report, p. 153. 

14) The formulae for the resistance of an ellipsoid have been derived by OBERBECK, 
and are given eg. in C. W. OSEEN, Hydrodynamik (Leipzig 1927), pp. 186—189, — 
See also: J. PERRIN, Journ, de Physique et le Radium (WAIN) Sy, 402) (iia) arava 77, 
Nal (Gg), 

15) See ‘The Ultracentrifuge”’, p. 5. 


Proc. Kon. Ned, Akad, v. Wetensch., Amsterdam, Voi. XLIII, 1940. Dil 


310 


correct order of magnitude is obtained, but that on the average: Reyer es ees 
— 1,3 Scaic, aS will be seen from the following table: 16) 


TABLE I. 
(Elongated ellipsoids, unhydrated). 
p= < lisa = 13 | 1913, 
Name of the protein M ee a F = 108 L ed + er ee Som 
Gliadin 27000 32200 ee ois) AOS) | SOON Mee ail See |) alee] Bal 
Lactoglobulin 38000 | 47200 5.98) 1021) 210) 20-81) 2.645)" 2238) 3.12 
Ovalbumin 40500 | 50100 Sn (h0)\| WOLO | AO ||P 5 | 2 || BeBe! Sis5)D 
Haemoglobin 68000 84100 5.38| 9.1 | 237|26.1| 2.66 | 3.67| 4.41 
Serum albumin 68000 | 84000 G57 || 1s | 269 | 24.5) 2.67 | Sigsee|) S656) 
Serum globulin 150000 184400 OO | Mes | SAAN) POLO Ae || spate) 7/51! 
Amandin 330000 407000 7.04| 11.8 | 476|40.4| 2.69 | 9.8 |12.5 
Thyroglobulin 650000 770000 New i Mera.) YANE)| SS Gy |) PR || WS. 5) || WEIL 
Homarus haemocyanin | 800000 977000 HESS) || WOK} |) WOES.) Hei US. | 22.6 
Octopus %) 2800000 3420000 Sy! WS) VU Oona |) Zot | SIO” 8) 33 
Helix pomatia _,, 6700000 8160000 OHO WOW | NANO) IS) |) BET | AS 98.9 


3. The discrepancy between calculated and observed results thus again 


turns up before us, but now in a form in which a better judgment can be 
made. 


In order to explain this discrepancy SVEDBERG, PEDERSEN and PoLsoNn 
have advanced the hypothesis that the protein molecules in the solution 
might be hydrated to such an extent, that their actual volume is equal to 
about 1,59 times the value MV/N 4,17). It is logical to assume that the 
consequent increase in dimensions then also must be taken into account 
in calculating the value of the frictional constant. It has been pointed out 
by KRAEMER that the hydration practically does not affect the driving 
force acting on a protein molecule in an aqueous solution during the 

16) The data used in the calculations mostly have been taken from table 48, p. 406, 
of “The Ultracentrifuge”. For the molecular weight the value of M, has been taken, 
with the exception of that of Octopus haemocyanin, where M, is used. The values of 
Ns ple have been derived from POLSON's data, Kolloid-Zeitschr, 88, 58, Tab. VI (1939); 
in the cases of serum globulin and Helix pomatia haemocyanin, however, the values of 
NsplcV have been derived from POLSON's Tab. III, Lic. p. 56, as it was not evident 
how the most suitable mean value should be obtained from the numbers given in Tab, VI. 

17) See A. POLSON, Kolloid-Zeitschr, 88, 56 (1939); K. O. PEDERSEN, The Ultra- 
centrifuge, p. 43, — According to POLSON the so-called electroviscous effect can be 


neglected under suitably chosen conditions (see also PEDERSEN, The Ultracentrifuge, 
p. 26). 


Syl! 


sedimentation process 18); nor does it affect the number of molecules 
present in a solution of a given concentration of c grams of dry weight 
per unit of volume. Hence we can use the equation: 


Hsp _ Na aL d? 
a a cee eee 416) 


together with equation (5) also in the present case. Instead of eq. (2), we 
now, however, take: 


GL PN Nee NN) VAIN 


When the calculations are repeated in this way, a smaller value is 
obtained for Lid; L decreases slightly, whereas d increases, but the value 
of the sedimentation constant, as calculated from (5) with the new values 
of L and d, practically remains unchanged, so that the discrepancy is not 
removed in this way. The results have been given in Table II. 


TABLE II. 
(Elongated ellipsoids; volume increased by 59% in consequence of hydration). 
; Ib, 
Name of the protein a 108 L 108 d LO’ Sate 
or 

Gliadin Aah DY 18.8 TGS 
Lactoglobulin Dod 168 29.2 Dead 
Ovalbumin D4 a5) 163 3055 2.62 
Haemoglobin 4.80 181 37.6 3.79 
Serum albumin (5555) 222 33.9 Bil 
Serum globulin 9555 oval 381.9 313 
Amandin iD 400 315). 5 9.9 
Thyroglobulin LORS 636 60.7 Loe3 
Homarus haemocyanin 6.4 495 UU os 18.9 
Octopus o 6 985 103 38.3 
Helix pomatia _,, 6.3 995 158 79 


It is possible to arrange the calculations differently, making use neither 
of eq. (2), nor of eq. (7), but solving L and d directly from eqs. (6) and 
(5). The results arrived at in this way, however, appear to be extremely 
unprobable, as they lead to excessive values of L/d and to very low values 
of d, whereas the volume turns out even smaller than MV/N,. As 
examples may be mentioned: 

ghadineseLid=— 3005 se 410 10-8; de=1, 14 )0-8 
Gvalbuming wid ==—/10;" 23560 10-8; d==0,51.. 10-8 
amandin« —Ljd=— 265; “b= 730.10—8;. d==2,76, 10-8 


18) KE, O, KRAEMER, The Ultracentrifuge, p. 63. . 
oh 
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It must be concluded, therefore, that the hypothesis of a large increase 
of volume in consequence of hydration does not help to remove the dis- 
crepancy between the results of the sedimentation and the viscosity 
measurements, so long as the assumption of an elongated ellipsoidal form 


is retained. 


4, When the values of the ratio L/d for the proteins are calculated 
from the observed values of the sedimentation constant S, as has been 
done by PoLson 19), results are obtained which are more than 50 % smaller 
than those found from 7jsp, as given in table I above. It might be suggested, 
therefore, as a possibility for an explanation of the discrepancy, that the 
molecules should have different shapes in the two types of experiments to 
which they are subjected. 

This possibility has been considered in the “Second Report” in con- 
nection with the application of the formulae to the results obtained with 
suspensions of polystyrenes, where a similar discrepancy between the 
calculated and the observed values of the sedimentation constant was 
found 20), It was shown there that the forces exerted by the liquid upon a 


19) A, POLSON, Kolloid-Zeitschr. 88, 56, Tab. III (1939). 

20) Second Report, pp. 176—178. The differences between the calculated and the 
observed values of the sedimentation constant S in the case of polystyrenes is seen from 
the table given at p. 176, where it must be kept in mind that the calculated values refer 
to the motion in the direction of the axis only, and consequently must be divided by 
about 1,5 in order to give the mean values for all directions of the axis in space. It 
will be observed that in the case of the polystyrene with M = 270000 the discrepancy 
is much greater than those which are found for the proteins. In the case of the suspensions 
of methyl cellulose, on the other hand, which were investigated by SIGNER and Vv. TAVEL 
(see footnote 2) above), the discrepancy between the calculated and the observed 
sedimentation velocities is much smaller. In ““Vhe Ultracentrifuge’, p. 437, Table 54, the 
values of V and 7,,,/c have been given; the values of L/d and of L and d calculated 
from these are collected in Table 55 at the same page. For the calculation of the 
sedimentation constant SIGNER has used a formula somewhat differing from the one 
applied in the present text (see Helv. Chim, Acta 21, 542 (1938) and “The Ultra- 
centrifuge’, p. 435); when the values are re-calculated with the aid of eqs. (5) and (4) 
given above, the following results are obtained: 


Fraction IV Ii] | II 
Molecular weight (My) 14100 24300 38100 
L/d 77 109 139 
LOSES yr 0.70 On7e 0.85 
IOS o 5s 0.83 0.79 0.89 


Here again the calculated values are less than the observed ones, although the difference 
is becoming smaller for the higher molecular weights, SIGNER remarks that the lengths of 
the molecules, calculated from Nsp exceed the lengths calculated from the number of 


os) 


molecule in consequence of the shearing motion which exists in a viscosity 
experiment, tend to produce an elongation of the molecule: on the other 
hand, during a sedimentation experiment there is practically no tendency 
towards a deformation of the molecule. Nevertheless, the fact that the 
viscosity of a protein solution does not appear to be markedly dependent 
upon the velocity gradient 21), makes it unprobable that large deformations 
are caused by the shearing motion of the liquid, unless it might be supposed 
that the molecule should possess two relatively stable forms, one of which 
would appear during the sedimentation experiments, while the other would 
appear in the viscosity experiments, even if the shearing velocity would 
be small. 


5. There remains, however, an important point which needs con- 
sideration, viz. that the formulae used in the calculations have been deduced 
for particles of ellipsoidal shape, and that numerically different results must 
be expected for particles having other shapes. Unhappily there do not 
exist exact formulae for particles of other shapes; the formulae given in 
the “Second Report” for cylindrical particles are approximations, which, 
although useful for great values of L/d, are not sufficiently precise for 
application in the problem here before us. 

It might be supposed that the molecules of the proteins should have the 
form of oblate rofational ellipsoids; for particles of such forms exact 
expressions can be given. Equations (1), (2) and (5) can be used also in 
this case; the values of Aj; and 2 then must be calculated anew. When it 
is assumed that the Brownian movement is sufficiently strong to make all 
directions of the axis in space equally probable, then for A; we can start 
ivom eq. (14.17) ot the Second Report (p. 151), provided for Cy-Ge, 
C, we now substitute expressions which are valid for oblate ellipsoids of 
revolution, and which also have been given by JEFFERY 22), like those for 
the elongated ones. A few results are collected in the following table, 
which is a counterpart to the table given at p. 153 of the “Second Report’’ 
for elongated ellipsoids (as before, we have written d= 2b for the equa- 
torial diameter, and L = 2a for the axial diameter of the ellipsoids): 


glucose residues present in the molecule. — In the case of the cellulose acetates, considered 
by KRAEMER and NICHOLS (‘The Ultracentrifuge’’, pp. 426—431), the opposite relation 
is found: here the lengths calculated from Wego remain below the maximum lengths 
deduced from the number of structural units, while the calculated sedimentation velocities 
exceed the observed ones, For other linear high polymers the calculated values of S 
again are too low; it is mentioned that here the molecules may be coiled, so that similar 
deformations might be possible as in the case of the polystyrenes. 

21) A. PoLsoN, Kolloid-Zeitschr. 88, 57, (1939). 

22) G. B, JEFFERY, The motion of ellipsoidal particles immersed in a viscous fluid, 
Proc. Roy. Soc. (London) A 102, 174—175 (1922—1923). The expressions for Cy, C2, C3 
can be derived from eq. (61), p. 174, in which the values given at p. 175, eqs. (68), 


have to be substituted for a9, 64, 65 (5 =o ay — Bg/2, according to p. We 


BAe 


Values of 4,,; for oblate rotational ellipsoids. 


bla = a|L Ary ba ai Ary 
1 a 2.50 12 6.39 

5 3, XO 15 7.64 

6 3595 20 9.74 

8 iit 25 11.85 

10 5.56 30 | 13.96 


From bla — 


is sufficiently accurate. 


35 onward the approximate formula 


Ar & 4b/3 wa + 1,19 


When the frictional constant in the present case is written in the form: 


fea Sac rid Nala (8) 
and the formula for the sedimentation constant as: 
0 a 
= 3a Nad (9) 
then the expression for 4 becomes, for oblate ellipsoids of revolution 
APM LD a) 
»— arcigl d Hep il (10) 


1 — L?/d? 


When the calculations are performed in the same way as was indicated 


in section 2, again assuming that the volume of a molecule is unaffected 
by hydration, the curious result is obtained that practically the same values 
of the sedimentation constant are found as had been derived upon the 
assumption of an elongated ellipsoidal shape. This will be seen from 


Table III. 


ABE ESI 


(Oblate ellipsoids, unhydrated). 


Name of the protein + 108 d 108 L 10S 5 
Gliadin Sil Fe 125 3.97 1.62 
Lactoglobulin iil 0) 100 9.1 2.47 
Ovalbumin 10.3 100 9.6 2.65 
Haemoglobin 9.6 116 12 85/3 
Serum albumin 2,3 125 10.2 Bo ay/ 
Serum globulin 18.3 186 10.2 238) 
Amandin 13555 219 16.2 10.0 
Thyroglobulin 20.3 310 1553 18).5) 
Homarus haemocyanin NP) (0) 282 2359) 19.2 
Octopus iM iss. 3} 493 PHL 50) 38.9 
Helix pomatia _,, 11.9 570 48.0 79 


MS) 


One still might wish to repeat the calculations for a different value of the volume of 
a molecule, in order to see whether the assumption of hydration might help us now. It 
could be attempted again to calculate L and d directly from eqs. (6) and (5). It is 
found, however, that upon the supposition of an oblate ellipsoidal form no solution of 
these equations can be obtained. By way of example in the following table the 
cases of gliadin and amandin have been considered. Starting with the value of 
Ld A,y,= (6a) . (75 p/¢) . (MIN ,) various assumptions are tried for the ratio d/L; for 
each of these values the corresponding values of L and d can be found, and the calculated 
value of L/k can be compared with the experimental value of this guantity, which is 


equal to (1—eV)/(3a) . (MIN4) . (iS pee): 


gliadin amandin 
experimental value of L d? 4,: 891000. 10—24 5450000. 10—24 
L/a : 2/47 lO=8 1, WO 
oblate ellipsoids 
value assumed for d/L calculated value of L/2 
400 SUROSS 149, 10-8 
100 81.10—-8 148.10—-8 
25 81.10-8 a MO 
5) TNO O83 
sphere 1 FN MO =8 130.10-8 
elongated ellipsoids 
value assumed for L/d 
2 T3eNOSS 134.10-8 
10 82.10—8 130), O83 
100 70.108 128.10-8 
200 65,6.10—8 EUS) ), O28 
400 Gil So IOs Wale, MO 


Hence it is found that the only simultaneous solutions of the equations are obtained 
for elongated ellipsoids with L/d between 200 and 400, which solutions already have been 
given at the end of section 3. 


(To be continued.) 


Mathematics. — Ueber lineare Linienkomplexe bei vier Geraden im Rj, 
Von R. WEITZENBOCK. 


(Communicated at the meeting of February 24, 1940.) 


Uebersicht. Ich bestimme im Folgenden den allgemeinsten linearen 
Linienkomplex, der kovariant zu vier Geraden in vierdimensionalen Raume 
gehért und beweise, dass durch vier Geraden allgemeiner Lage und ihre 
vier Transversalen, also durch eine sogenannte ,,Doppelvier” 7? lineare 
Komplexe gehen, unter denen drei spezielle vorhanden sind. Die Brenn- 
punkte der iibrigen fiillen eine Ebene aus. 


Ue 

Ich habe vor einiger Zeit bewiesen!) dass sich jede ganze rationale 
projektive Invariante von fiinf Geraden a, a, p, m und a im Rj, die wir 
auch durch die Ziffern 1, 2, 3, 4 und 5 andeuten, ganz und rational 

ausdriicken lassen durch die Typen 
A1,23,45 = (12? oid) (42 3) eae (S23 945) ik Zoe Tithe We (1) 

ik 

wobei z.B. S23 der Verbindungsraum der Geraden 2 und 3 ist. Ai,23,45—=0 
driickt also aus, dass die Schnittebene der beiden dreidimensionalen 
Raéume S23; und Sis von der Geraden 1 geschnitten wird. Die Invarianten 


Ai, 23,45 a Ai, 45,23 


sind zyklisch-symmetrisch in 123 und ebenso in 145. 

Die Frage nach dem allgemeinsten linearen Strahlenkomplex K, der 
zu vier Geraden 1 bis 4 kovariant ist, lautet dann in algebraischer 
Formulierung: es ist die allgemeinste projektive Invariante K der fiinf 
Geraden | bis 5a zu ermitteln, die in den 2x linear ist. Nach Obigem 
ist demnach K linear aufzubauen aus Invarianten Am, ix,rs, worin genau 
einer der Indizes ein a ist mit Koeffizienten, die Polynome der Am, ix, rs 
sind, wo alle fiinf Indizes aus 1 bis 4 genommen sind. Von diesen zuletzt 
genannten Invarianten habe ich bewiesen, dass sie auf vier algebraisch- 
unabhangige zu reduzieren sind, z.B. auf 


A, = Agi3.2 
A, = A123,24 
A; = Ao13,34 cutee, ey Ee (2) 
Ag = A314,24 


1) Proc. Kon. Ned. Akad. v. Wetensch., Amsterdam, 42, 245—252 (1939), 


ene 


Die Invarianten A mit einer Reihe 2 verteilen wir auf zwei Klassen: 
erstens Invarianten Jn Am,ix,-s in der alle fiinf Geraden vorkommen. 
Solche J, gibt es 15, die sich aber vermdge der zyklischen Symmetrie 
linear, mit konstanten Koeffizienten durch sechs linear-unabhangige 
ausdriicken lassen’), z.B. durch die folgenden sechs: 

= A1,23,47 = A, 13,40 
Jo Aisa sn b= eek Ae) 
| rae Ie SS Aas te 

Zur zweiten Klasse rechnen wir diejenigen Invarianten K;;—= A, in 
denen eine der Geraden 1 bis 4 fehlt. Macht man wieder von der 
zyklischen Symmetrie Gebrauch, so gibt es deren zwolf, namlich: 
Kg = Agia Ko, = As 14,40 K3, = Agia14 Ky = Ai, 23,24 
KGa == Ago K,, == Aa ie te Key = Aina 27 K43 = A2,13,32 


(4) 
Ky == A4o3. 30 5G@e == Ajai aa a= AQT 4 ct Ke= A312,12 


Hier gibt der erste Index r von K;; die Gerade an, die nicht in A 
vorkommt; die zweite Index s gibt die Gerade s, durch welche die 
Ebene mit der Gleichung K;;—=0 geht. K;;—=0 ist also die Ebene der 
beiden Geraden s und r*. Dabei ist z.B. 4° die Transversale von 1, 2 
und 3. 

Die Invarianten der Gestalt (4) mit als erstem Index sind durch die 
K,s von (4) ausdriickbar. So ist z.B. wegen der zyklischen Symmetrie 


Az,23,24 —— ta Ai423,22— A2,23,42 —= —K,, 


da Ao2342 verschwindet. 
Das allgemeinste K hat demnach die Gestalt 


6 
Ka Di Jp Dees Kas, ee aa : . 5 (5) 
1 rs 


wobei die Koeffizienten wu; und m;; Polynome der vier Invarianten A; 
von (2) sind. 


Se2, 


Wir beweisen jetzt, dass die Invarianten (3) auf die K;s von (4) 
reduzierbar sind. 
Es sei b mit a=1 Aquivalent; dann haben wir 


Pe As 74 75 sae = (42 ag 47) (b3* 47) (b32 n°). 
Hier bringen wir die Reihe b des dritten Faktors in den ersten: 


P= (b2? 3%) (a3? 42) (a2? 42)(b3. 2) + 2 (ab 23%) (a2? 4?) (23? 4?) (b3?2r”). 


2) Le. p. 251: 
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Der erste Term rechts ist hier Aj,34,24. Ai,23,32; beim zweiten vereinigen 
wir die beiden Reihen a im ersten Klammerfaktor und erhalten 


Zab 23") (a2? 42) =—taeo 4b 2 4), 
da (a?6..) zu Null fiihrt. Somit kommt, da 


Ad? 03,34 = Ai,23, 34 
ist 


A\1,23,24 Ai1,34,30 + Al, 24,34 A1,23,32 = — Aj,23,34 Ai, 24,32 een (O) 


oder, in den A; und K;s ausgedriickt : 


A 
L2a3 wet As 43 A; ( ) 


Durch Permutation der Ziffern 1, 2, 3, 4 ergibt sich hieraus fiir alle 
Invarianten (3) Reduktion auf die K;s von (4). 

Es ist fiir das Folgende niitzlich die Wirkung der sechs Transpo- 
sitionen i><k der vier Ziffern 1, 2, 3, 4 auf die Invarianten A; und 
K,s zu tibersehen. Wir stellen dies in einer Tabelle zusammen, wobei 
wir links die Invarianten und oben die Transpositionen als Eingangs- 
reihen setzen. Bei den K;s geben wir in der Tabelle nur das Indexpaar, 
lassen also den Buchstaben K weg. So findet man also z.B. in der mit 
K,; beginnenden Zeile unter der Transposition 23 in der Tabelle 
— 42 stehen. Das heisst: vertauscht man in K,,; die zweite mit der 
dritten Geraden, so entsteht — K4p. 


12 es 1 4 ee 24 3 4 
A, —A, ek As Ay —A, —A, 
A, —A, —A, —A, —A, —A, —A, 
A; —A, —A, —A, —A, —A, —A, 
A, =A, —A, —A, —A, =A, ae 
Kis “23 —31 143 a2) (3 =) 
Ky, +24 —34 +41 —14 —12 =13 
Ky say —32 +42 —13 —14 ENE 
Ko, qel4 —24 —21 +34 —42 =93 
Ky, =i2 ——20 —24 +31 a4 25) 
Kos IBN 2 2 ae 
Ks Ee —13 st +21 —31 +41 
K3> =i)! —]2 —32 “23 —34 4:49 
Ky —4] —42 +12 —43 —24 +32 
K43 —43 —4] +13 — 42 —23 +34 
Ky =42 3 +14 =41 —21 34 


ol) 


Durch Vertauschung von 1 mit 3 ergibt sich so aus (7) 


A A 
A z——K. ZTE K 2 
3, 24,1 27 1 aA, 


Addieren wir dies zu (7), so entsteht wegen 


Ax,24,13 = Ai,24, oC A\3 94 100? 


LA get oe — K,, a + Ky, = —K,, A (8b) 
und hieraus, wieder durch Vertauschungen: 
Abe a KG Kn e o (8a) 
Ly = Anan = —Ka 5? — Ke & Ks. $4 Kn 3. (8c) 
53 


Wir konnen jetzt im allgemeinen Ansatz (5) fiir den Komplex K die 
Invarianten J weglassen und 


Ks A la A eK ge. oc te Ay ee Os + ’ . (9) 
rs 


setzen. 
Substituieren wir hierin der Reihe nach ~=—1, 2, 3 und 4, driicken 


also die Bedingungen dafiir aus, dass die Geraden 1 bis 4 dem Komplex 
K angehGren, so ergeben sich die Gleichungen 


eee) eA ee Aiea Oe (10,) 
Ae a Ae ee (0. (10,) 
Pe ee Aa 0 ee 0) 
hip As hy As ele Ay 0. (10,) 


Wir nennen wieder 1* die Transversale der Geraden 2, 3 und 4. 
Dann lauten die Bedingungen, dass die Geraden 1” bis 4° den Komplex 
K angehoren: 


qn ay qin =0- 1 PS) 
hat a de qa=0. 2S Pies 
ao a hst ate =0. 7 oe ea 
Alu qelat giao... ~ (108) 
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Diese acht Gleichungen (10) und (10°) sind, wie man leicht nach- 
rechnet, durch eine lineare Abhangigkeit verbunden. Es ist namlich, 
identisch in den zw6lf 4;s: 


1 1 1 
A. (10,) + (103) + 4 (10s) + (109) = / a 


= A, (10) )—A, (103) + Az (103 )—A4 (105). , 


ee 


Die Gleichungen K;;=0 stellen zwélf Ebenen des R, dar. Da zwischen 
elf Ebenen im Ry, stets wenigstens eine lineare Beziehung existiert, 
miissen zwischen den zw6lf Komitanten (4) wenigstens zwei lineare, 
identisch in allen Verdnderlichen 2;, geltende Beziehungen vorhanden 
sein. Diese wollen wir jetzt ermitteln. 

Zunachst gelten ftir die A;; mit 


KES Ki 0 “in-allen-waipie os 2 ee ee) 

die Gleichungen (10): sie entstehen ja dadurch, dass wir ax = 1 bis 4 
und = 1" bis 4° setzen. 

Der Schnittpunkt der beiden Geraden 1 und 4" ist gegeben durch 


a;(aS23), der der beiden Geraden 2 und 1* durch a; (a.S3,4): die Ver- 
bindingslinie dieser beiden Punkte hat daher die Koordinaten 


Tie = (a a)ie a S53) (AiSs4) 
Setzt man dies in (12) ein, so erhalt man 
bop AG Ay A — eA 
Hier gibt die Vertauschung von 3 mit 4 nach der Tabelle 
dy Ay Al ie Ag lA A) 
Hieraus und aus (103) kann man /,; eliminieren und erhdlt 


A, As; “ 
ea 


Die sechs Vertauschungen i><k fiihren nach der Tabelle zu fiinf 
weiteren Gleichungen: 


Ang =S= Ais 


A,A Ag 
Aig =—dgg os AV Men re 
va A;A AA 
We Ty bi x eee, 
A34 = Ay» eee 


ING 
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Setzt man dies alles in die Gleichungen (10) ein, so findet man, dass 
sich auch die iibrigen 4,; durch 2,;; und 4,3; ausdriicken lassen; man erhilt: 


: AY Ee A; AG 
= —_ foe Sai | ase I — /\ ae 2 i == 2 
ae 2 3a 34 13 AVA | 23 As 7M, 
A,A A A A 
A= Nie mee 3 A 3 a alee: 
21 13 A 23 ie 43 13 A, dys As (13a) 
A A,A 
Nan == — jh g ye aa ee 
32 23 A 41 23 Ave 


Mit diesen Werten werden auch die Gleichungen (10*) identisch in 2,; 
und 4,3; befriedigt. 

Setzen wir jetzt diese 4,;; in (12) ein, so entsteht, wenn wir alle 
Nenner beseitigen: 


K=A,A,A; A, (i, A>M,,—4,;, A; M,;)= 


woraus 


M,;=0 und M,;=0 fin allen wiz}. . . . (14) 


folgen. Hierbei sind die M;; Ausdriicke, die sich mit Hilfe der Abkiirzung 


aK 5 (15) 

wie folgt darstellen lassen: 
M,3—=21—12 + 13+ 31—24— 42 + 43—34. . . (16a) 
Ma; = 2112 = 14 = 41 + 32—23 + 43—34. . . (166) 


Bei den Transpositionen i<k bleibt M3; invariant oder geht in —M; 
oder in M,;—/M,, iiber. 

Man beweist leicht, dass die zwei Identitaten (14) von einander unab- 
hdngig sind und dass alle linearen Abhangigkeiten der K;s auf diese 
beiden Identitéten zuriickfiihrbar werden. Man kann also nach (16a) und 
(166) z.B. Ky, und K,3 linear-homogen durch die zehn iibrigen K;s von 
(4) ausdriicken. Die allgemeinste, in den 7; lineare Komitante lautet dann 


eG gg ee sa a fay a ee) 
TS 


und es ist, in der Bezeichnung (15): 
eee ee ee ia 04 (18) 
AR 104i 14-4 23 — 324-34. 2. (19) 
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§ 5. 


Wenn K von (17) durch die Geraden 1 bis 4 und 1* bis 4° gehen 
soll, gelten fiir die 4,; die Gleichungen (10), wobei noch 442 = 443 =0 
genommen ist. Aus (10) folgt dann auch 4,,=0 und die iibrigen 
Gleichungen fiihren auf 


1 1 
3 Az + Ag Ag+ 442 A, =0 ee a 
] I 
dog Ay + Ay, Ay + 4,3 Az = 0 Wiese cia tea =U 
A, A; 
; (20) 
dz, Ay + 432 Az + d34 Ag = 0 A dos Ao =e 
l ] 1 
As Ais B, dog 4 a A434 = 0, 


Hieraus folgt, dass sich alle 4,;; durch drei unabhangige, z.B. 4,2, 4 
und 4,3; ausdriicken lassen: 


A, 


. by = — Fh Be dys By = ly lg=0 
an yee Sai K A332 = — oh Ag @ 
K ng = Eh ly = — GSE Anes in Wei) 
Somit wird K: 
Kin) Ky 2 Ki 3 Ka 4 Ara Kad + 
td — SK +28 Ky — BA Kx tKa{+o. . (2 
+ is} Kis 52 Ky — S253 Kg +? Key 


Wenn wir in den Ausdriicken L; der Gleichungen (8) vermdge (18) 
und (19) die 42 und 43 eliminieren, so entsteht 


£,= 13 + 23—14—24, L,= Mae oe te ee 44 (22) 


und dies sind, bis auf Faktoren = die Koeffizienten der 2,5 in (21). Wir 
k 


erhalten somit den Satz: 
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Durch die acht Geraden 1 bis 4 und 1* bis 4* einer Doppelvier des 


R; gehen oc? lineare Strahlenkomplexe. Sie sind gegeben durch die 
Schaar 


My Ly oA; Late ds Lae A, Aig 34 a > A igo AG Ais — an 25) 
Hier ist z.B. 
Ala, 12,34 = (912) (22 S34) = 0 


ein spezieller Komplex, namlich die Schnittebene der beiden Réaume 
S12 und Ss4. 

Die 27 Strahlenkomplexe (23) haben oo? Brennpunkte, die auf vier 
Ebenen liegen. Drei davon sind die Ebenen L;—O selbst und diese 
Ebenen stellen die einzigen speziellen Komplexe der Schaar dar. Die 
vierte Ebene wird gegeben durch die Verbindungsebene der drei 
Schnittpunkte L; Lz, hat also die Gleichung 


(a0 Sis S04 Si4 S53) (x0 Si4 S93 Sio S34) (a0 Sia S34 Si3 S24) == () 
Nach einiger Rechnung fiihrt dies auf 
= 4-4) 32 30, se (24) 


wobei die Bedeutung von rs aus (15) zu entnehmen ist. Man rechnet 
weiters leicht nach. dass diese Ebene E weder die Geraden 1 bis 4 
noch die hierzu assoziierte schneidet. 

Wenn wir von Komplex K nur verlangen, dass er durch die vier 
Geraden 1 bis 4 geht, so gelten fiir die i,; nur die ersten drei Glei- 
chungen (20). Man kann dann z.B. 4)2, 423 und 434 durch die iibrigen 
Ars ausdriicken und erhalt fiir K die Gleichung 


dys (12-13) + Ag (1412) + dag (24-23) + doy (21-23) + ) 


(25) 
+ ds, (31—34) + dg (3234) = 0. | 
Hier ist z.B. 
A A 
12-13 =Kn a ine Anan a As222— FB Ar3432=0 ; (26) 
ein Komplex durch alle vier Geraden und nach (25) gibt es also — wie es 
sein muss — 0° solche lineare Komplexe. 


Es ergibt sich von hier aus auch eine Lésung fiir die Aufgabe: 
Ermittlung der Gleichung eines linearen Strahlenkomplexes K, der durch 
neun Geraden 


Aik, Aik» +. ++ Vik 
gegeben ist. In Gestalt einer zehnreihigen Determinante 


Dm Cy .q7. 05, 045) == 0 + * + + + . (27) 
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ist natiirlich diese Gleichung sofort hin zu schreiben. Wie kann aber K 
durch die quindren Invarianten dargestellt werden? Hier bekommen wir, 
wenn wir in (25) fiir die ;x die fiinfte bis neunte Gerade einsetzen, fiir 
die sechs homogenen /;s fiinf Gleichungen, kénnen diese 1,; also elimi- 
nieren und erhalten eine sechsreihige Determinante fiir K: 


(2-315) (1412) eee ae 
(12-13), (14212), 5 ee ee 

oe | ae eee - | =0 . (28) 
(1213) 4 120 easing 


Hier bedeutet z.B. (12—13);, dass fiir i, in (26) die Koordinaten der 
fiinften Geraden eingesetzt sind. Wir ersehen hieraus, dass 


(A, A, A; A,)®.K 


ein Polynom von Komitanten des Typus (1) wird. 


Mathematics. — Ueber assoziierte Geraden bei Regelflachen im R, 
Von R. WEITZENBOCK. 


(Communicated at the meeting of February 24, 1940.) 


Uebersicht, Zu vier erzeugenden Geraden G, bis G, einer zweidi- 
mensionalen Regelflache F im vierdimensionalen Raume R, lasst sich die 
wassoziierte’ Gerade G, konstruieren. Wenn die vier G; gegen eine 
Erzeugende a von F konvergieren, nahert sich auch die assoZiierte 
Gerade G; einer Grenzlage g. Man erhalt so — von Ausnahmen abge- 
sehen — zu jeder Erzeugenden a von F eine assoziierte Gerade Ga asic 
bilden in ihrer Gesamtheit die zu F ,,assoziierte Regelflache’”’. 

Ich ermittle hier die Gleichung von g unter der Voraussetzung, dass 
die homogenen Linienkoordinaten a;, der Erzeugenden a von F als 
regulare Funktionen eines Parameters gegeben sind. Die hierzu erforder- 
lichen Rechnungen sind nicht ganz miihelos, man hat bis zu den Gliedern 
dreissigster Ordnung in den benutzten Potenzreihen vorzudringen. 


et 


Die Geraden einer einparametrigen Schaar bildem im R, die Erzeu- 
genden einer zweidimensionalen Regelflache F. Wir denken uns F 
analytisch gegeben durch die zehn homogenen Koordinaten a;, einer 
Erzeugenden, wobei diese aj, stetige und geniigend oft differenzierbare 
Funktionen eines komplexen Parameters ¢ sind: 


ae ae) (eet, De 5) 


Da die aj, Linienkoordinaten sind, verschwinden die fiinf quadratischen 
Ausdriicke 


= 
2 ays Ags, 2) Ais Ay5,.4-52 412 434 


fiir alle Werte von t. Wenn b mit a Aquivalent ist, lasst sich dies 
symbolisch ausdriicken durch das Nullsein der Kovariante 


4. DX a3 a45 = (xa? 6b?) = Mn =0 1 oe Ok 


Hieraus ergeben sich durch fortgesetztes Differenzieren nach f eine 
Reihe von Gleichungen, die fiir das Folgende von grundlegender Be- 
deutung sind. 


Wir setzen 
d? Aik (t) 
(a, )ix = Tee , (ar)ix — Sore nee 


Proc. Kon, Ned. Akad. v. Wetensch., Amsterdam, Vol. XLIII, 1940. OY) 
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und deuten a,a,,a),... an durch die Ziffern 0,1,2,.... Dann ergibt 
sich, identisch in f¢: 


Moo = 0 
on) 
Mii + Mo2 = 0 
3 M12 + Mo3 = 0 
3 Mx + 4Mi3-+ Mos = 0 
10 Ma; + 5 Mis + Mos = 0 se al) 


10 M33 = 15 Mo = 6 Mis le Mos = 0 

35 M34 a5 75) Mos =e 7 M6 = Mo = 0 

35 Mas + 56 M35 + 28 Mos + 8 Mi7 + Mos = 0 
126 M4; + 84 Mz. + 36 Moy + 9 Mis + Moo = 0 


Hier sind die (a,)ix, (a2)ix,.... die Koordinaten von linearen Ebenen- 
komplexen K,,K,.... Der Komplex Kn ist speziell, stellt also eine 
Gerade dar, wenn Mj,=0 gilt. So ist z.B. der Komplex K, speziell 
fir Mi; =—Mw=0; wegen Mo. —0 enthalt K, jede Ebene durch 
die Erzeugende air. 

Wenn die Flache F eine Torse ist, d.h. ihre Erzeugenden die Tangenten 
einer Raumkurve 


Gi oe0) Fi aecitl Wa A te 5) 
sind, so findet man 
y )ik 
y" )ik 
: y" )ik = (y He AYE: 
(a3)ixk = ae yin Ely ol iz 


u. s. £. 


woraus sich ausser Moo = 0, Mo, 0 noch weiters ergibt: 
Mii == 0; Mo2 == 0), M;? == 0; Mos ==), 


Dagegen sind die Mix mit i+ k=4 im Allgemeinen nicht Null. 
Liegt F ganzlich in einem linearen R;, so ist Mo2 bis auf einen von ft 
abhangigen Faktor eine Linearform in x mit konstanten Koeffizienten: 


M07 =o 2) ee) (0) ee ee ee ee ( D) 


Wir wollen im Folgenden diese beiden Falle ausschliessen, also vor- 


aussetzen, dass die Flache F weder abwickelbar sei noch zur Gdnze in 
einem linearen R; liege. 


SPATE 


Se) 
Vier Geraden 1 bis 4 mit den Koordinaten ajp, aix, Pik und mix 
haben im R, die vier Invarianten 
ma’ p*) (ma? a?) = (m3) (mS4}) 
aa? p?) (aa? m?) = (aS23) (aS24) 


( 
(aa 

3 = (aa? p”) (ap? m?) = (aS13) (aS34) | 
( 


(3) 


pa?’ m’) (pa? m?) = (pS14) (pSx) 


und die eres der zu 1 bis 4 assoziierten Geraden G; wird gegeben 
durch '): 


G;=G,A,A;A,+ G,A, A;A,+G;A,A,A,+G,A,A,A;,;=0. (4) 


Wir nehmen jetzt als Geraden 1 bis 4 vier Erzeugende von F, ent- 
sprechend den vier Parameterwerten a, 6, y und 6. Es ist dann also 


a a2 a a? 
Gi = ik + Fj (ai)ie + 5; (a)ix one Op Typ lik toy 2ue +... 


2 5) 
Gr=an+ Fla) tS (a2)iz + 0+ Flt Stet .. 6) 


u.s.f. mit y und 0. 


Fir diese G; wollen wir (4) berechnen. Hierzu haben wir vorerst die 
Ai notig. 

Wir berechnen A,. Aus (5) finden wir zuerst fiir Sj2, wenn wir aus 
allen Gliedern den Faktor $(a— f)? absondern: 


S2=Mat| 3+ AMS] +] i+ A? M+ 3 300 Mis | + 


+ | Gye + Pe? + a8 + 8M +55 oP e+ A) Mi | + 


1 2 2)\2 17! 27! 1 Dh 
+ | 5p 8 + B+ PP Mis + 35 oP (e+ BP Mis + 35 PMs | + 
+ | 7p FA + 898 + a8 + Mert 
+ 355 08 a+ Ao +08 + #9) Mia + 755 Plat 8) Mis | + 
+ | gy REL PP Mie + 77 BO? +B + BP Mi + 


tf 1 t 
+ 79 FU A Ms 56q 0 Ms] + | 


1) Vogl. diese Proceedings 42, 248 (1939). . 
22 
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Wir haben hier die Glieder gleich hohen Grades in a, f zwischen 
eckige Klammern gesetzt. Aus (6) lesen wir ab, dass 
(x Mo2) = (xa? a2) = (x0? 27) —0 


die Gleichung des langs der Erzeugenken aj beriihrenden R, darstellt. 

Vertauschen wir in (6) f mit y, so ergibt sich Sj3 und damit lasst 
sich nach der ersten Gleichung der Formeln (3) A, berechnen. Wenn 
wir zundchst mj. allgemein lassen. so ergibt sich nach Abspaltung des 
Faktors +(f—y), also bis auf den Faktor 


THO) a OP) es eee ee 


fiir (mS{3)(mS12) die Reihe 
(m S13) (m S12) = 102,03 =i E (2 ha B =S 7) ™02, 04 Ir cen ata | 


+ fe (a + B) (a + y) mo3,04 + ; a? mo3,13 + 


fa? + a(B + y)} mo2,14 + 
l 
+ 9g 12a? + 2008-4 1+ BP By +77} monos | + 


ol) 
+ | app 2h FeO + N+P + B+ 7) manos + (8) 


1 

= 10 Nees 2a(p + y) 4 BP eyes mo2,15 -+ 
1 1 

“i 8 a? {8 > y) Mo2,24 + 60 (a =P) (a + 7) (es pt) ™o3,05 + 
1 

i fois (a + B) (a+ y) mo3,14 + = a (a*—By) mon se 

=O E ™Mo2,07 + @1 ™Mo2,16 + @2 Mo2,25 + = a? By mis | eee 


In (8) haben wir schliesslich nach (5) fiir mix die Reihe 
) OF 
Qik P77 lik + 57 2ik 1... 


einzusetzen um die verlangte Potenzreihe fiir A, zu erhalten. Die Rechnung 
ergibt, wenn wir fiir das Differenzenprodukt 


HT == {a—B) (a—y) (a0) (B=) (B= 0) (yee) ae) 
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die Abkiirzung I gebrauchen und 
Q = 013,25 = (aa? a) (aa? a?) 2 2. |. . (10) 
setzen: 
Ae Ol Gites ie aes 


Hier bedeuten die Punkte Glieder vom zehnten und hdheren Grade in 
dena, 8, y, 0. 

Wenn wir z.B. a mit £ vertauschen, werden nach (5) G, und G, ver- 
tauscht, A, geht dann in —A,, A, in —A,, A; in —A, und A, in 
—A, tiber. Auf diese Weise ergeben sich fiir die A; aus der zuletzt 
angeschriebenen Gleichung die Reihen: 


A, =— te Q. IT. (a—f) (a7) (a0) +. | 
PO a ACL EAC MM) (11) 
poms a ye NY ee 
Nee EG) 00 (aa aN 


Bevor wir nach Gleichung (4) G; berechnen, sind einige Bemerkungen 
beziiglich der Grésse Q von (10) zu machen. 

Sind ax, Gik, Pik, Mix und zx die Koordinaten von fiinf linearen 
Komplexen, so ist 


J = Gap, mx = (aa? p*) (am? 27) = 16. D aik Gi, i, Dizi, Mkyks Tkek 


eine projektive Komitante dieser Komplexe. Sie ist symmetrisch in a und p 
und in m und x, schief-symmetrisch in den Paaren ap und mz: 


J = aap, mx = —- Ama,ap. 
Schliesslich ist J zyklisch-symmetrisch in a,a,p und ebenso in a,m, a: 


aap,ma + @pa,mx + Pas,ma = 0 


und ebenso mit a, m,2 
Diese Symmetrieeigenschaften im Vereine mit den Glechnanen (1) 


ergeben eine Reihe von Reduktionen bei denjenigen Ausdriicken, die 
aus J entstehen, wenn wir an Stelle der 


ik» Aiky Diks +++-+ die Reihen aik — Oik, (a;)ix — 1 aes (a) ix a Qik, oo OY 


nehmen. Es sei dies bei der durch die Gleichung (10) eingefiihrten 
Grésse @ ndher ausgefiihrt. 
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Wenn wir in (8) 
6 oe 
Mik — Oi t Wi lik + 57 ik ++ 
einsetzen, so werden die Glieder vierten Grades in a, 6, y, 0, da 
Ooi, 0 


fiir alle i, r und s gilt (die Oj. sind Linienkoordinaten!), gegeben durch 


1 
O1= 12 a? By 013,14 | 
+o E 102,06 + 92 102,15 + 93 102,24 + %%4 103,05 -- 


ar * a (a + B) (a + 7) 103,04 + 3 a (a*—By) ls | i 
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+ 3j 


ie (a + B) (a + y) 203,04 + ; a? 203,13 + eet te) 


= j (a? aa (6 af y)) 202,14 + 65 202,05 | ot 


6? {1 

ap 3 E (2a =- 8 =. y) 50z04 4 «on + 
jt 

ei 4/ P 402,03. 


Die hier auftretenden Koeffizienten K: 
ijk,mn G+ j-k-- m+ n= 9) 


sind alle auf das Q von Gleichung (10) reduzierbar oder verschwinden. 
Sei z.B. K = 202,05. Die zyklische Symmetrie beziiglich 2, 0, 5 gibt hier 


202,05 —= — 02,25 — 5o2,02 = 0, 


denn 0o2,05 verschwindet, da 0; Linienkoordinaten sind und 502,02 ist 
Null wegen der schiefen Symmetrie in den Paaren 02,02. 
Wenn K = 102,15 ist, so haben wir zufolge der dritten der Gleichungen (1) 


K = 1o215 = — litis 


und dies ist Null, da allgemein r,,,:,—=0 gilt, wie aus der zyklischen 
Symmetrie in r, r, r sofort folgt. 


Bei K = 402,03 gibt die zyklische Symmetrie beziiglich 4, 0, 2: 


Ke 402,03 Ss 024,03 =s 204,03 = 204, 03+ 
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Aus (1) haben wir 
04 = — 3,22— 4.13, 


also wird 
402,03 = 3.222,03 ar 4.21303 = 4.21303. 
Hier ist, wieder wegen der zyklischen Symmetrie beziiglich 2, 1,3: 
213,03 == 123,03 — 312,03. 


Der letzte Term wird nach (1) wegen 03—=—3.12 gleich 3.32 12, also Null. 


Beim ersten Term der rechten Seite gibt die zyklische Symmetrie be- 
ziiglich 1, 0, 3: 


193,03 i 023,13 — 323,01 = 013,23 = @; 


da 323,01 wegen 01 nach (1) wegfallt. Also wird schliesslich 
213,03 = — Q und 4203 = —4Q. 


Auf diese Weise sind alle Koeffizienten in (12) zu reduzieren. Man erhalt: 


013,14 —= —2Q 

102,06 ===1(0); 102,15 ==), 102,04 10) 103,05 ==) 103,14 =, lo4.13 a2) 

Zee — 4) 20343 == ©), 209 4 = — 4 Q), 20795 == 0 (13) 
302,04 = 4), 302,13 —C 

402,03 =— 4. 


Setzt man dies in (12) ein, so wird 
Opa O(a—0)7(f--0) (7-0)... ss (14) 
was schliesslich mit (7) zusammen zu den Gleichungen (11) fiihrt. Diese 


letzteren ergeben dann 


Ay As Ay= — 755 Q I (87) (88) (78) +. 
A, As Ay= + 755 Q I (0-7) (0-8) 78) + 
(15) 


Ay A, As=— gps O(a B) (a—6) (B—6) +.... 


Ay Ay As= + 753 QI (af) (a—7) (B-7) + 


wo die Punkte Glieder andeuten, die wenigstens vom 28, Grade in 
a, 8, y, 6 sind. 
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Man zeigt leicht an einem Beispiel, dass die Invariante Q = 013,23 im 
allgemeinen von Null verschieden ist. 


Schliesslich setzen wir (15) in (4) ein, wobei die G; durch die Reihen (5) 
gegeben sind. Als Glieder niedrigster Ordnung ergeben sich in den 
Reihen, mit denen 

Mik = Vikn Liga cipeee os 
pe: 
2° 
Sm = a™ (B—y) (B—8) (y—6)—B™ (a—y) (a—6) (y—0) +- } 

+ y™ (a—f) (a—8) (p—8)—0" (a—f) (a—y) (87). J 


Man rechnet leicht nach, dass Sp, S; und S, identisch wegfallen, dass 
dagegen S,;0 wird, sodass die Reihenentwicklung fiir G; mit den 
Gliedern dreissigster Ordnung beginnt. Hieraus schliesst man, dass bei 


Q+#~0 


Q@? II* die Ausdriicke 


zu multiplizieren sind, bis auf den Faktor — 


(16) 


lim G; (a0, B 0, y > 0, 6 > 0) = giz 
eine Linearkombination von Ojz, liz, 2ix und 37, sein wird, wobei 3 ix 
mit dem Koeffizienten Q? versehen ist. 
Se 
Man wird so zu dem Ansatz 
Jie Op An lie ie oie en Le) 


gefiihrt. Die Bedingung dafiir, dass die gi, Linienkoordinaten sind, lassen 
sich dann nach (1) zusammenfassen zu 


My 2? + Mr we? + M33? + 24 Moo + 2» Mo3 + 
+ 2AuMi2 + 2,4» M33 = 2 nv Mp3 == 
Driickt man hier Mi, Mm und Miz nach (1) durch die iibrigen M ix 


aus, so bekommt man 

Moz (2 ¢—2) + Mos (2»— $4.) —$ 0? Mos + Mis(24 —§ 2) + ) 
; (18) 

+ 2»Mo;+»2 M3; —0. 9 


Dies ist eine lineare Beziehung zwischen den sechs Raumen R;, deren 
Koordinaten die Mj, sind. Nennen wir in der Matrix 


|| Moz Mos Mos M3 Mp; M3; | 
die fiinfreihigen Determinaten O;,, so gilt 
» Mit Oik = 992 Mo + @3 Mo3+...+0;;M3=0. . (19) 


und dies muss bis auf einen Proportionalitatsfaktor mit (18) identisch 
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sein. Es gilt also vorerst die Ausdriicke ©; zu berechnen. Wir setzen 
(Vgl. (10)): 
Q = 01323, R=04333, S = 093,33, Dee 204,33. rae te (20) 


Dann ist z.B. 


ee ae ie 


=4,0)3 33 20301 —=— 160 RQ. 


Auf diese Weise bekommt man nach einiger Rechnung: 


OQ. =4(4RT—S?) 0,;—16QS | 
Cyt Qi Rs) 0,;—16QR (21) 
Cy a Cr, — 16,0” \ 


Wahlen wir also als Proportionalitatsfaktor +7, so liefert die Ge- 
gentiberstellung von (18) und (19) die sechs Gleichungen 


2u—2? =o? (t RT—S?) 
—2y + 4ip=o?(QT+$RS) 


22 Dp? 
fad Q R (22) 
2Av—$ ww’? = 407 QS 
0) IX. 
—407Q? 


Hieraus folgen 
C= POR , w= 200 


womit auch die vorletzte Gleichung (22) befriedigt ist. Aus (22), und 
(22), kann man dann e@ und 4 berechnen: 
= 30 '@? jue 12ZOW3OS-- R?) 
eo 2RLI2ZORS+IQ?T *° 2R?+12Q0RS+9O°T 
Mit diesen Werten wird auch (22), befriedigt, sodass sich schliesslich 
unter der Voraussetzung 


Dene ORS 9.O Wea0 alee ae (23) 


fiir die gesuchte Gerade g ergibt: 
Jik — ik (2R3+12QRS+9Q?T)—(ai)ix 12 Q(3QS+R?’) + cS (24) 
+ (a,)iz . 36 Q? R—(as)ix 72 Q? J 


Diese Gerade gix ist die ,,assoziierte’’ der Erzeugenden air. 


Chemistry. — Der Einfluss des Dispersitatsgrades auf physikalisch- 
chemische Konstanten. (Neunte Mitteilung.) Von ERNST COHEN 
und J. J. A. BLEKKINGH Jr. 


(Communicated at the meeting of February 24, 1940.) 


Polarographische Analyse der bei 25.00° C. gesattigten Losung 


des Bariumsulfats. 


Finleitung. 


In unserer siebenten!) und achten?) Mitteilung wurde der Nachweis 
gefiihrt, dass der Léslichkeit des Bariumsulfats in Wasser bei bestimmter 
Temperatur und bestimmtem Drucke (1 Atm.) ein bestimmter Wert 
zukommt, welcher nicht abhangt von dem Dispersitatsgrade des Boden- 
kérpers. Dieser Beweis wurde erbracht durch Messung des elektrischen 
Leitvermégens der gesattigten Losung des chemisch und physikalisch reinen 
Sulfats. In der vorliegenden Abhandlung beschreiben wir unsere Versuche 
zur Erhartung dieses Ergebnisses auf polarographischem Wege. Zu diesem 
Zweck verwandten wir das polarographische Analysenverfahren, welches 
uns auch in stand setzte, die Ba~lonenkonzentration der gesattigten Lésung 
quantitativ zu ermitteln, und aus dieser die Léslichkeit des Bariumsulfats 
in Wasser zu berechnen. Samtliche Versuche wurden bei 25.00° C. (und 
1 Atm. Druck) ausgefiihrt. 

Ueber die Bariumbestimmung auf polarographischem Wege haben 


J. HEYROVSKY und Mitarbeiter 3) bereits einige Untersuchungen angestellt. 
Indes entsprechen die Schliisse, welche sie daraus betreffs der Léslichkeit 
des Bariumsulfats und des ,,Alterns’” des Niederschlags ziehen, dem 
Tatbestande keineswegs, da erstens das von ihnen verwendete Barium- 


sulfat nicht chemisch rein, zweitens aber ihre Lésungen nicht vorher 
filtriert waren. 


A. Die Apparatur. 
1. Eine ausfiihrliche Beschreibung des auch von uns verwendeten 


Verfahrens gibt J. HEYRovSKY in dem Kapitel ,,Polarographie’’ des Werkes 
Physikalische Methoden der analytischen Chemie” von W. BOTTGER, 


1) Proc. Kon. Ned. Akad, v. Wetensch., Amsterdam, 43, 32 (1940). 
?) Proc. Kon. Ned. Akad. v. Wetensch., Amsterdam, 43, 189 (1940). 
3) Collection 1, 36 (1929); Phil. Mag. (6) 45, 303 (1923) 


. 
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Leipzig 1936, S. 260—323 1). Fiir nahere Angaben verweisen wir auf die 
Dissertation von J. J. A. BLEKKINGH Jr., Utrecht 1938. 


Der von uns verwendete Apparat von HEYROVSKY war von der Firma 
V. und J. NEJEDLY in Prag verfertigt. Das benutzte Spulengalvanometer 
war nach JULIUS 2) erschiitterungsfrei aufgehangt. Es liessen sich bei 
maximaler Empfindlichkeit Stromstarken bis 10-9 Amp. bestimmen. Die 
weiter unten abgebildeten Polarogramme sind in ihrer urspriinglichen 
Grosse reproduziert worden. (Empfindlichkeit des Galvanometers 2/s). 
Unsere Tropfelektrode bestand aus einer mit doppelt destilliertem Queck- 
silber gefiillten Kapillare aus Jenaer Glas, welche zu einer feinen Spitze 
ausgezogen war. Wir verwandten bei samtlichen Versuchen dieselbe 
Kapillare, wahrend das Quecksilbergefass des Apparates sich stets 64 cm 
iiber dem unteren Ende derselben befand. In dieser Weise war es méglich 
reproduzierbare Ergebnisse zu erzielen; dies ergibt sich aus unseren 
Figuren, welche sich auf Versuche beziehen, in denen stets zweimal 
dieselbe Lésung polarographiert wurde. Das Quecksilber aus dem Behalter 
str6mte langsam durch einen griindlich gereinigten, 1 m langen, dick- 
wandigen Schlauch aus transparantem Pararubber zur Kapillare. 


B. Tetramethylammoniumjodid als indifferenter Elektrolyt. 


2. Bei den Bestimmungen entfernten wir den Sauerstoff aus unseren 
Loésungen indem wir wahrend 1 bis 2 Std., in manchen Fallen auch wahrend 
langerer Zeit, reinen Wasserstoff durch dieselben leiteten. Es wird namlich 
der Sauerstoff eher abgeschieden (bei —0.4 und —1.1 Volt) als das 
Barium (bei —1.95 Volt), was zu einer erhéhten Stromstarke Anlass 
geben wiirde. 

Als indifferenten Elektrolyten versuchten wir zunachst (nach J. 


HEYROVSKY) Lithiumchlorid bezw. Lithiumchlorat (die Abscheidung des 
Lit erfolgt bei —2.3 Volt). Es gelang uns indes nicht hiermit befriedigende 
Resultate zu erzielen, wohl infolge der in den Lithiumsalzen vorhandenen 
Verunreinigung. Aber auch ohne Zusatz eines indifferenten Elektrolyten 
gelang es uns nicht unser Ziel zu erreichen. 

Dagegen stellte sich heraus, dass eine Losung des Tetramethylammonium- 
jodids (dieselbe zeigt erst bei -—2.6 Volt eine ausgesprochene Leitfahigkeit) 
sich bei unseren Bestimmungen vorziiglich als indifferenter Elektrolyt 
eignete. Zwar enthielt das Praéparat (pro analysi) von KAHLBAUM, auch 


nach sorgfaltigem Umkristallisieren aus L. W.?) noch spuren iee—n(c) 


und NH? — Ionen, wie sich aus Fig. 1 aus der Abweichung von einer 


1) Neuere Literatur: ,,De polarographische methode met de druppelende kwikelektrode 
ten dienste van het pharmaceutisch onderzoek” von J, MAAS, Dissertation Amsterdam, 
1937, sowie auch ,,Chemische Analyse mit dem Polarographen’ von H. Houn, Berlin 
1937. 

2) Ann, Physik (3) 56, 151 (1895). 

3) Vergl. unsere achte Mitteilung, § 1. 
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geraden Linie ergibt, welche bereits bei —1.6 Volt anfangt, sowie aus der 
schwachen Welle bei —2.1 Volt. 

Wir stellten uns eine 0.01 n Lésung des Jodids in L.W. her, und gaben 
stets 45 cc der zu analysierenden Lésung zu 5 cc dieser Jodidlésung. 
Von dieser Mischung brachten wir etwa 20 cc in das W.G. und leiteten 
sodann wahrend wenigstens 1 Std. Wasserstoff durch diese Lésung. 

Bei samtlichen Versuchen ist die Endkonzentration somit 0.001 n an 
Tetramethylammoniumjodid. Unsere Fig. 1 stellt das Ergebnis des blinden 
Versuchs mit dieser Lésung dar. 


s, <ps —i¥ | =I V 


Fig. 1. 0.001 n. Tetramethylammonium- Fig. 2. Wasser, filtriert durch ein nicht 
jodidlésung bei 25.00° C. vollig gereinigtes Glasfilter, durch wel- 
ches vorher eine Bariumsulfatlésung 

filtriert war. 


Ein Ausserlich ziemlich sorgfaltig gereinigtes Glasfilter gab beim 
Filtrieren von L. W. durch dasselbe noch Ba++-Ionen an das Wasser ab, 
welche von vorher auf diesem Filter filtriertem Bariumsulfat herriihrten. 
Dies zeigt sich deutlich in der Fig. 2 an der schwachen Welle bei 
—1.95 Volt. Wird das Glasfilter aber mit warmer, konzentrierter 
Schwefelsaure behandelt, und sodann griindlich mit L. W. ausgewaschen, 
so erhalt man ein mit Fig. 1 véllig identisches Polarogramm. 


C. Die Eichung mit Bariumchloridlésungen. 


3. Bevor wir das polarographische Verfahren auf gesattigte Barium- 
sulfatlésungen anwandten, analysierten wir auf diesem Wege einige 
Bariumchloridlésungen bekannter Konzentration. Zu diesem Zwecke 
stellten wir uns Lésungen desselben her, welche 10—5n bezw. 4 X 10-5 n 
waren. Wir benutzten BaCl,.2H,O (KAHLBAUM, p.a.), von welchem wir 


SoH 


zunachst festgestellt hatten, dass es genau 2 Mole Kristallwasser enthielt. 
Zwar enthalt dieses Praparat 0.1% Sr++, sowie Cat+, dies hat indes 
fiir unseren Fall keine Bedeutung. 


Die Fig. 3 und 4 zeigen die erhaltenen Polarogramme in ihren wahren 
Dimensionen. 
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Fig. 3. 10—° n. Bariumchloridlésung. Fig. 4. 4 < 10—° n. Bariumchloridlésung. 


D. Versuche mit filtrierten, gesattigten Bariumsulfatlosungen. 


4. Wie bei den Messungen des elektrischen Leitvermégens dieser 
Losungen beschrieben wurde (achte Mitteilung), stellten wir uns Suspen- 
sionen von nicht zerriebenem, chemisch und physikalisch reinem BaSO, 
her, sowie solche des Sulfats, welches im Achatmorser zerrieben war. 
Im letzteren Falle verwandten wir stets 0.7 g des trocknen Sulfats, 
welches wir in 150 cc der gesattigten Lésung gaben. Sodann saugten wir 
die Suspension auf einem Glasfilter ab (achte Mitteilung § 6), und setzten 
zu 45 cc des Filtrates so schnell als méglich 5 cc der 0.01 n Tetramethyl- 
ammoniumjodidlésung. 

Sollte die Konzentration des Bat+ im Filtrat des zerriebenen Sulfats 
eventuell eine gréssere sein, so wird dieselbe durch dieses Verdiinnen nicht 
mehr verschwinden kénnen, da die Lésung ungesattigt geworden ist. 

Unsere Fig. 5 und 6 zeigen die Polarogramme der Lésungen, welche in 
der beschriebenen Art und Weise dargestellt wurden mittels des nicht 
zerriebenen, bezw. des zerriebenen Bariumsulfats. 

Eine Differenz zwischen der Bat+-Ionenkonzentration in diesen und 
vielen andern Versuchen, war nicht nachweisbar, auch nicht dann, wenn 
die Lésung wahrend einiger Zeit ruhig gestanden hatte. 


5. Aus den Mittelwerten des Grenzstromes (der Hédhe der Welle) 
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und denen bei den Bariumchloridlésungen (Fig. 3 und 4), berechneten wit, 
dass die bei 25.00° C. gesattigte BaSO,-Lésung 2.08 X 10~°n ist, somit 
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Fig. 5. Nicht-zerriebenes BaSO4 Fig. 6. Zerriebenes BaSO« 


Gesattigte BaSO4-Lésungen in Wasser bei 25.00° C. 


2.43 mg gelésten BaSOy, pro Liter enthalt. Dabei wurde angenommen, dass 
das sich in Lésung befindliche BaSO, vollig elektrolytisch dissoziiert ist. 
Fiir das Loslichkeitsprodukt (Aktivitatsprodukt) Bat+ x SO”, bei 
25.00° C. ergibt sich hieraus der Wert 1.08 K 10-19. 

Dieser Wert stimmt iiberein mit demjenigen, welchen A. C. MELCHER 1) 
mittels des elektrischen Leitvermégens gefunden hat. 


E. Versuche mit nicht-filtrierten Bariumsulfatsuspensionen. 


6. Auch einige Versuche mit nicht-filtrierten Suspensionen des Barium- 
sulfats, welches zuvor nicht, bezw. wohl zerrieben war, haben wir 
ausgefiihrt. Nachdem wahrend 1 Std. Wasserstoff eingeleitet war, nahmen 
wir die Polarogramme auf, welche in Fig. 7 und 8 reproduziert sind. 

Auch hier ergeben sich gleiche Ba++-Ionenkonzentrationen. Die zum 
Bodenkérper gehérigen Gegenionen werden offenbar nicht bei derselben 
Potentialdifferenz ausgeschieden, wie die freien Bat+-Ionen der Lésung. 

Handelt es sich darum die Konzentration des Sulfats mit der in den 
friiheren Bestimmungen (in welchen 45cc gesattigter Lésung verwendet 
wurden) zu vergleichen, so haben wir den 45/50 Teil des Ergebnisses in 
Rechnung zu ziehen, da jetzt ein Ueberschuss an Bodenkérper vorhanden 
ist, so dass sich 50 cc einer gesattigten Lésung bilden. Wir finden dann, 
dass diese Lésung 2.31 X 10-5 n ist, somit 2.71 mg Bariumsulfat pro Liter 


1) J. Am. Chem. Soc. 32, 50 (1910). 
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enthalt. Dabei ist wiederum angenommen worden, dass vollige Dissoziation 
vorliegt. 


Die Léslichkeit des BaSOy, in einer 0.001 n Tetramethylammoniumjodid- 
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Fig. 7. Nicht-zerriebenes BaSO4 Fig. 8. Zerriebenes BaSOg 


Suspensionen von BaSOx4 in einer 0.001 n. wasserigen Tetramethylammonium- 
jodidl6sung bei 25.00° C. 


lésung ist somit bei 25.00° C. (2.71 2.43) —0.28 mg pro Liter, hdher 
(etwa 10%), als die in Wasser. Auch hieraus ergibt sich, dass die Schliisse 


von J. HEyrovsky 1) dem Tatbestande nicht entsprechen. 


Zusammenfassung. 


In unserer 7., 8. und 9, Mitteilung wurde sowohl mittels Leitfahigkeits- 
messungen, wie auch auf polarographisckem Wege der Nachweis erbracht, 
dass sich ein Einfluss des Dispersitatsgrades auf die Léslichkeit eines 
kristallisierten, chemisch und physikalisch reinen Stoffes (BaSO,) nicht 
feststellen lasst, falls man durch geeignetes Filtrieren der gesattigten 
Lésung desselben dem Auftreten des ,,Dispersitatseffektes” vorbeugt. 

Mittels polarographischer Analyse erfolgte dann weiter die Feststellung, 
dass die Léslichkeit des Bariumsulfats bei 25.00° C. 2.43 mg pro Liter 
betragt, und dementsprechend das Loslichkeitsprodukt 1.08 x 10~10. 


VAN 'T HoFF-Laboratorium. 
Utrecht, Februar 1940. 


1) Collection 1, 36 (1929); Phil. Mag. (6) 45, 303 (1923). 


Geology. — Remarks on the geology of Colombia and Venezuela. I. 
The age of the non-fossiliferous slates and of the metamorphic 
schists. By L. RUTTEN. 


(Communicated at the meeting of February 24, 1940.) 


The author has been in the necessity of studying, during the past months, 
the geologic literature on Colombia and Venezuela. He has got, thereby, 
some ideas which differ materially from the conceptions generally held 
with regard to the geology of these republics. He regrets that he has paid 
only short visits to the republics — the longest having been an excursion 
ta NW. Venezuela and Cucuta, where he travelled as a guest of the 
Caribbean Petroleum Company — but he is of the opinion that this lack 
of local experience needs not impede him to form and to publish personal 
conceptions on the geology of the area. There are especially three subjects 
which the author would like to treat: the age of the non-fossiliferous slates 
and sandstones and of the metamorphic schists, the areal distribution of 
the cretaceous, and the historic geology and the tectonics of the tertiary. 

The metamorphics of Wenezuelan Guyana, which so clearly belong to 
the Guyanan-Brazilian shield, are generally regarded as prepalaeozoic, 
although this age jhas not been proved with absolute certainty (45). In the 
following we shall not discuss them; we shall occupy ourselves exclusively 
with the metamorphics of the cordilleran region. 

The slates and schists of the Cordilleras have been regarded by many 
authors as palaeozoic or pre-palaeozoic. W. BERGT (5, 6, 7) called the 
slates palaeozoic and the schists prepalaeozoic; W. SIEVERS (36, 37, 38) 
considers slates and schists both to be archaean; E. GROSSE (13), under 
the influence of G. STEINMANN, regards all the slates and schists of 
Antioquia as prepalaeozoic and thinks that their age increases with 
increasing metamorphism; E. HUBACH (17) and E. SCHEIBE (30) consider 
the Quetame-schists to the East of Bogota as praecambrian and H. GERTH 
has, in the most general sense, pronounced that slates and schists in 
Northern South America are prepalaeozoic. He writes: 


,.Die Hauptmasse der fossilfreien regionalmetamorphen Gesteine, die durch die ganzen 
Kordilleren ...... immer wieder an die Oberflache kommen, halte ich fiir algonkisch. 
Wiederholt sehen wir ...... palaozoische oder mesozoische Schichten mit den meta- 
morphen Gesteinen verfaltet, ohne dass die erstere ihre Fossilfiihrung eingebiisst 
INBIDYESM soonor Wohl werden auch die jiingeren Sedimente bei intensiver Faltung dynamo- 
metamorph verandert, aber dann nehmen sie doch nicht den Habitus der hier behandelten 
regionalmetamorphen Gesteine an” (11, p. 80—91). 


34] 


Two years later, GERTH is somewhat more prudent in his opinion, at 
Jeast with regard to the schists of the Caribbean Coast Range, where 
LAMARE has found — amidst the schists — true radiolarites (24). GERTH 
writes: 


Wenn diese (the radiolarites) auch nicht gerade als absolut beweisend fiir ein 
paldozoisches Alter der betreffenden Gesteine angesehen werden konnen, méchten wir 
uns nicht mehr so entschieden gegen das palaozoische Alter gewisser metamorpher 
Gesteine in der nérdlichen Kette des karibischen Gebirges aussprechen, wie wir es 
friiher getan haben” (12, p. 346). 


Whereas the authors hitherto mentioned apparently try to give the rocks 
in question a very high age, there are others who have regarded them as 
palaeozoic or, partly, even as mesozoic, EUGSTER (10) calls the meta- 
morphics of the Cordillera de Bogota palaeozoic; R, SCHEIBE (32) equally 
considers some schists of the Cordillera Central of Colombia as palaeozoic; 
De CIZANCOURT and SCHNEEGANS (9) do the same with respect to the 
Venezuelan metamorphics. REINHARD (28), who is strongly convinced of 
the postalgonkian age of the schists of the Coast-Range in NE. Venezuela, 
discusses the possibility of their mesozoic age. 

Various authors, who described slates and schists from Colombia and 
Venezuela, have not given an opinion about their age. Among these we 
may mention AUBERT DE LA RUg (3), BENDRaT (4), HETTNER-LINCK (16), 
RUTTEN (29) and STUTZER (40). 

In general, the hitherto mentioned authors did not give thorough 
arguments for their opinions; we may say that their opinion has been 
based chiefly on a certain ‘‘feeling”. Two arguments, however, have 
been put forward from time to time. Firstly that there are known from 
Colombia and Venezuela non-metamorphic, fossiliferous, carboniferous, 
devonian, and, according to some geologists, non-metamorphic ordovician 
deposits, and that therefore all the metamorphics must be predevonian 
resp. pre-ordovician. Secondly that the lack of fossils in the slates and 
schists under discussion proves their great age. It is desirable to argue 
that these two “proofs” have only very small value. As to the first 
argument we may cite that, on the island of Aruba (44) there have been 
found hornblendeschists of cretaceous age which are strongly metamorphic 
and which, according to their metamorphism would f.i. be regarded by 
GROSSE (13) as to be of archaean age. Moreover, there occur in many 
countries at the side of non-metamorphic, older rocks younger ones which 
are more or less metamorphic. As to the second argument it is sufficient 
to indicate that f.i. the Cambrian of the Ardennes and the Silurian of 
Brabant are slightly metamorphic and that they are so poorly fossiliferous, 
that, if they occurred in a country so slightly investigated as Colombia or 
Venezuela, they certainly would not have yielded fossils until now. 

We shall now pass to the discussion of the different regions, for which 
there have been given somewhat better arguments for the age of their 


metamorphics. 


Proc, Kon. Ned, Akad, v. Wetensch., Amsterdam, Vol, XLIII, 1940. 23 
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1. W. KeHRER (22) found in the southernmost Colombia two 
metamorphic series: ‘‘Glanzschiefer’’, phyllitic slates and quartzitic sand- 
stones are resting with an angular unconformity on sericite-chlorite-schists, 
phyllites and quartzitic conglomerates. He proves with these interesting 
observations that there exist certainly two independent formations of schists 
in this part of the country, of which the older one presents the strongest 
metamorphism. He regards the older formation as Algonkian and the 
younger one as possibly Cretaceous, but for this interpretation he does 
not give convincing arguments. As a matter of fact the algonkian age 
remains absolutely unproved, as palaeozoic rocks in this part of the country 
have not been found until now. 

2. R. ScuetBe (31) has found in Antioquia a formation of porphyritic 
lava’s and tuffs, the cretaceous age of which has been proved by some 
intercalated fossiliferous layers. In sections, where this porphyritical 
formation outcrops, he finds also a series of green schists, and he feels sure 
that these schists are of the same age. E. GROSSE (13) regarded some 
years later, on the base of their metamorphism, these schists as of middle- 
algonkian age. I should say that SCHEIBE’s opinion is by far more 
convincing, but that he has not proved the cretaceous age of his green 
schists. As he has not shown — by careful sampling and by microscopical 
study of the samples — that there exists a gradual transition between the 
members of the porphyritic formation and the green schists, it remains 
possible that there are two formations of quite different age, which, in 
consequence of intense folding or thrusting, have been brought into 
seemingly normal contact. 

3. HETTNER (15) described from the area E. of Bogota a series of blue 
clayshales, quartzites and quartzitic conglomerates which he called 
“Quetame-series’”” and which he regarded as probably pre-cretaceous: they 
occur below fossiliferous, cretaceous rocks, STUTZER found to the W. of 
this zone of ancient rocks (40) fossiliferous, non-metamorphic carboniferous 
strata, which contain, however, according to W. KEHRER (21) true 
quartzites. In the neighbourhood have been found also non-metamorphic, 
devonian rocks (33). The Quetame schists are, therefore, very probably 
pre-devonian. W. KEHRER (22) calls them even pre-palaeozoic, and, 
although this is quite possible, it is far from being proved. We have fi. 
to consider that there have been found in Colombia ordovician fossils in 
a series of rocks, which, according to SCHUCHERT (33), are “strongly 
folded quartzites, slates, marbles and black limestone’, whilst, in 
Venezuela, ordovician fossils have been found in ‘‘a partially metamorphic 
shale and sandstone series, in some places approaching a micaceous schist” 
(TERRY, teste SCHUCHERT, 33, p. 693). L. KEHRER, equally, calls the rocks 
of the Venezuelan Ordovician “‘semimetamorphic’’ (20). It is, therefore, 
very well possible that the semimetamorphic Quetame-schists are silurian. 

4. The strongly metamorphic rocks of the Eastern Cordillera N. of 
Bogotaé may very well be pre-palaeozoic. It is, however, necessary to 
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consider that there have been found in the Cordilleras during the last twenty 
years more and more outcrops of “hercynian"’ rocks, and that it is not 
impossible that the strongly granitized rocks N. of Bogota are palaeozoic, 
and metamorphosed by hercynian granites. LLERAS Copazzi (26) even 
regards the granites of this region — without proving his opinion — as 
“Andengesteine’’, i.e. as old-tertiary or as youngest-cretaceous. 

5. The silver-white phyllites which have been found as float in rivers 
on the eastern side of the Sierra de Perija (25) are very probably 
pre-devonian, because in these same rivers the devonian is fossiliferous 
and non-metamorphic. It can, however, not been proved for the present 
that these phyllites are prepalaeozoic. 

6. The gneisses of the Sierra Nevada de Santa Marta may very well 
be pre-palaeozoic, but they might equally well be palaeozoic, their meta- 
morphism being due to hercynic granitization or to hercynic dynamometa- 
morphism (5, 36). 

7. The schists of Goajira peninsula (40) are certainly pre-cretaceous, 
their pebbles having been found in cretaceous conglomerates: it is, however, 
impossible to say, even approximately, how much older than the Cretaceous 
they are, 

8. The area in Venezuela where the metamorphics have been best 
studied, is the Cordillera de Merida. In the past thirteen years there have 


been published four stratigraphical surveys, which are reproduced here in 
a diagram. This diagram shows, that great differences in interpretation are 
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still in existence. There are only two pre-mesozoic formations, the Palma- 
rito series and the Caparro series, as to which no disaccord exists among the 
four authors, their age having been proved by carboniferous and ordovician 
fossils. The Bellavista series, united by three authors with the Caparro, 
is separated from it by KUNDIG (23) and regarded as older. The 
Mucuchachi phyllites are taken by CurisT (8) for Archaean, by OPPENHEIM 
(27) for Cambrian, by KUNDIG (23) for eo-palaeozoic. The fact that 
L. KEHRER (20) has found traces of fossils in this series — he mentions 
even the occurrence of ammonites — renders a pre-cambrian age highly 
improbable. KiUnpiG states that the rocks of the Mucuchachi series 
pass gradually into those of the Bellavista series, into micaschists and 
into true contactrocks on approaching a contact with post-Mucuchachi 
granites. It is clear that, if this is true, it becomes improbable that the 
Bellavista schists and the micaschists are pre-palaeozoic. As possibly 
pre-palaeozoic rocks there remain then only the rocks of the Iglesias series 
(KUNDIG): gneisses, migmatites and arterites. KUNDIG regards them as 
archaean “Unterbau-Gesteine’. He does not prove, however, that the 
Mucuchachi-Bellavista-Micaschists-series and the Iglesias-series really 
belong to different geologic cycles. His sections nowhere show a normal 
contact between the rocks of the first and second group, nor does 
he mention the existence of conglomerates in the Mucuchachi, containing 
pebbles of the Iglesias series. As long as no unconformity has been 
found between the Iglesias and the Mucuchachi series nor pebbles of 
Iglesias rocks in Mucuchachi conglomerates, we must take into con- 
sideration the possibility that both series have the same (?palaeozoic) 
age, the rocks of the first series having been metamorphosed more intensely 
by stronger granitization and migmatization in deeper zones of the 
palaeozoic geosyncline. 

9. LIDDLE has tried to prove that a great part of the schists of the 
Coast Range of Venezuela is of silurian age (25). His arguments, however, 
are by no means convincing, as he departs from the wrong supposition 
that limestones of the Coast Range have yielded silurian fossils (29). 

10. LIDDLE has been the first to publish data on the occurrence of 
metamorphic, cretaceous rocks in Northern Venezuela (25) between the 
meridian of Barquisimeto and the eastern part of the Coast Range. He 
states that in this area cretaceous sandstones, shales and limestones have 
been converted into quartzites, schists and marbles. His contention that the 
metamorphosis is due to the influence of an intrusion of augite-porphyrite 
(25, p. 68—69) is certainly wrong. SCHURMANN (34) equally states that 
in the neighbourhood of Barquisimeto even the eocene rocks show signs of 
beginning metamorphosis, and that rocks of cretaceous age have been 
changed into silvery-white micaschists. L. KEHRER (19) says that the 
cretaceous limestones in NE and E. Lara have been changed into marbles, 
and the shales, marls and shales with organic substance into sericite-schists 
and graphite-schists. The statements of these three independent and 
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competent geologists are so positive that it seems impossible to avoid the 
conclusion that in this part of Venezuela there occur on a large scale 
metamorphics of cretaceous age. This conclusion, however, is so important 
that we can not feel satisfied by the general statements, published until now. 
Only a map, on which the transition from normal to metamorphic sediments 
in the direction of the strike of one or more structures is indicated, 
accompanied by a detailed description and by microphotographs of samples 
from crucial localities might satisfy completely our curiosity. It may be 
added that the author is not the only one who is longing for more evidence 
in this important matter. GERTH (11) has cited with some suspicion 
LIDDLE’s opinion about the cretaceous metamorphics, and, in the second 
volume of his “Geologie Siidamerikas” (12, p. 347), when SCHURMANN’s 
publication had already appeared, he still regards the schists of the 
Serrania del Interior as parts of the basal complex of Venezuela. 

I may add, that, during an excursion in 1930, | collected myself some 
cretaceous rocks in Venezuela. The description of three samples follows 
here: 


Lower Cretaceous Tomon Sandstone from Agua Caliente near San Antonio, SW. 
Venezuela. The rock is a white, somewhat porous, quartzitic sandstone, The slide presents 
in some places a very course, quartzitic structure, In other parts of the slide the rounded, 
original quartzgrains are still visible, each quartzgrain having, however, been enlarged 
by apposition, so that the general structure has become equally quartzitic. There are 
also places, where the pores between the original, rounded quartz grains have been filled 
with quartz which has an orientation of its own. 


Lower Cretaceous Tomon Sandstone from La Cuchilla between Carora and Valera. 
A grey, quartzitic rock. The slide shows quartzitic structure throughout. In some places, 
however, the original “‘sand-grains” are still recognizable. Between the quartzes there 
are to be found many small flakes of sericite, sometimes passing into semi-coarse 
muscovite; moreover there occur some idiomorphic crystals of tourmaline. The rock is 
a sericite-quartzite with tourmaline. 


Lower Cretaceous Tomon Sandstone from La Cuchilla between Carora and Valera. 
A red-yellow mottled, feebly schistous rock. The slide shows an irregular fabric of quartz 
grains and sericite flakes; the latter sometimes passing into muscovite. Small grains of 
limonite are responsible for the red colour. The rock is a slaty quartz-sericite schist. 


From the foregoing description it is to be seen that the Cretaceous 
presents traces of metamorphism, even far to the W. of Barquisimeto. The 
described samples may be reckoned to belong to the uppermost part of the 
epizone. These facts seem to support the views of LIDDLE, SCHURMANN and 
L. KEHRER. 

11. If I rightly understand the data of the geologic literature, the rocks, 
treated-with sub 10 all belong to the Serrania del Interior. SCHURMANN 
(35) as well as L. KEHRER (19) are of the opinion, that there occur equally 
cretaceous metamorphics in the northernmost range of the Caribbean 
mountains, and, according to LIDDLE (25, p. 75), there are to be found on 
the island of Margarita at the side of “silurian” schists also cretaceous 
metamorphics. As LIDDLE himself states that it is hardly possible to 
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differentiate these two complexes of schists, we may neglect for the moment 
his contentions. 

The arguments of SCHURMANN and L. KEHRER are the following. The 
radiolarite described by LAMARE (24) from the northern part of the Coast 
Range might very well be a cretaceous rock, the Cretaceous of Venezuela 
abounding in many localities with radiolarites; the marbles, so common in 
the northern parts of the Caribbean Mountains might be metamorphosed 
cretaceous limestones; the limestone of Valencia — from which LIDDLE 
has erroneously mentioned silurian fossils — might be a cretaceous lime- 
stone. As to the last argument the authors might find support in a very old 
publication of KARSTEN (18), who, in a diagrammatic section, indicates the 
Valencia-limestone as “Exogyrenkalk’’ and who, in his text, points to the 
similarity of this limestone with that of the Morros de San Juan — on 
the southern side of the Serrania del Interior — in which, according to 
LIDDLE, Exogyras have been found. TRECHMANN’s discovery (42, 43) of 
rocks with cretaceous fossils in the Northern Range of the island of 
Trinidad, and SPATH’s description of tithonic ammonites from the same 
area (39), which seem to indicate that the metamorphic rocks of this range 
are cretaceous, might equally support the views of SCHURMANN and 
KEHRER, HEDBERG (14), mentioning TRECHMANN’s discovery, has indicated, 
on the other hand, that it is difficult to understand that in the easternmost 
part of the Coast-Range of Venezuela there should occur at the side of 
strongly metamorphic, cretaceous rocks other rocks of the same age which 
are absolutely non-metamorphic. Dr. ‘TRECHMANN was so kind as to provide 
me with some samples from the neighbourhood of his fossil-locality in 
NE. Trinidad. Their description follows here. 


Conglomerate, W. of Jetty, Toco, NE. Trinidad. It consists of pebbles of more than 
1 cm with a sandy cement. The elements of the conglomerate are chiefly quartz and 
quartzite. Some small pebbles of a limonitic rock which contains quartz-splinters and 
muscovite-flakes, The cement contains quartz, muscovite and some calcite. 


Dark, calcareous shale, E. of Jetty, Toco, NE. Trinidad. Very dark rock with an 
unmistakeable echinid spine. The slide shows numerous organic remains (echinid frag- 
ments, calcified sponge spiculae etc.) in a brown, dense, isotropic ““groundmass”. There 
is no trace of metamorphism, 


Dark, calcareous shale, W. of Jetty, Toco, NE. Trinidad. The slide shows remains of 


organisms in a dark “groundmass’’ with quartz splinters and mica flakes. There is no 
trace of metamorphism. 


Conglomeratic sandstone, S. of Toco, NE. Trinidad. The slide reveals the presence 
of the following elements: 1. very much quartz, partly cataclastic, 2. finegrained quartzite 
to quartzschist, sometimes with muscovite, 3. grains of perthite, 4. grains of orthoclase, 
4. coarse quartzite. The cementing material consists of much quartz and sericite, The 
rock contains metamorphic material and has itself not suffered from metamorphosis. 


Black limestone E. of Jetty, Toco, NE. Trinidad, There is — in the slide — not a 
trace of metamorphism in this very finegrained limestone; it contains siliceous sponge 
spiculae, and grains and some crystals of pyrite. 
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Conglomeratic, calcareous sandstone, 4 miles S. of Jetty, Toco, NE. Trinidad. The 
elements of the conglomerate are: 1. microcline, 2, abundant guartz, partly cataclastic, 
4, grains of an intergrowth of quartz and plagioklase, 4. quartzite with some muscovite, 
5. very dusty plagioclase. 


Apparently there occur in this part of Trinidad non-metamorphic 
limestones and grits, which contain pebbles of metamorphic rocks, It has 
been known for a long time, that there occur also in the Northern Range 
of Trinidad true metamorphics. It is, therefore, not impossible that, in 
N. Trinidad, we have to do with a folded basement with in-folded jurassic 
and cretaceous rocks. 

AGUERREVERE and ZULOAGA (1, 2) have, in two very recent publications, 
described the Serrania del Interior and the Coast Range s. str. They 
distinguish to the S. of Caracas two different groups, both metamorphic: 
the Caracas group and the Villa de Cura group. They are of the opinion 
that there occur many cretaceous rocks in the Serrania del Interiér. Their 
conclusions are, however, somewhat contradictory. In their first publication 
they say that the Caracas group is cretaceous, in the second that the 
Caracas group is cretaceous or older, the Villa de Cura group cretaceous 
or younger. With regard to the rocks of the Coast Range s. str. they say 
at some place that they might equally be cretaceous, at other places, that 
they are possibly archaean. It seems better, for the moment, to leave out 
of this discussion the work of AGUERREVERE and ZULOAGA. 

Whilst many observations seem to indicate that probably a part of the 
metamorphics of the “Serrania del Interiér” is cretaceous and that possibly 
some of the metamorphics of the Coast Range s. str. may be cretaceous, 
there is one publication which indicates that there occur probably also 
much older elements in the Coast Range. IJZERMAN (45) has studied a 
series of granitodioritic rocks of the Venezuelan cordilleras and comes to 
the conclusion that: 


“the granitodiorites of the Andean folding area reveal less constant features than 
those of the old core of South America ...... On the one hand they are characterized by 
greater variabilities, on the other hand many tynves cannot be differentiated from those 
belonging to the old core ...... It is a fact of more significance that the typical orthite 
and the primary epidote we met with in many countries of the old core, seem to be 
absent in the Andean folding area’’ (45, p. 492—493). 


Now, LAMARE (24, p. 142) has described from the Silla de Caracas 
in the Coast Range s. str. a granite, which, according to his description 
can hardly be distinguished from the granites of the Old Core; it contains 
large microclines, primary epidote and orthite. This discovery indicates 
that there may exist in the Coast Range, at the side of possibly much 
younger rocks, very old, precambrian elements. 

The facts, discussed in the foregoing, allow for the following conclusions. 

1. The age of by far the most non-fossiliferous slates and schists in 
the Columbian-Venezuelan Cordilleras is unknown; it is not permitted to 


348 


speak — as has been done very frequently — of an extensive precambrian 
basement. 

2. There is only one case, where the precambrian age of rocks is 
rather probable: the granite of the Silla de Caracas. 

3. Many schists are certainly precretaceous (Goajira, Sierra de Perija, 
Sierra Nevada de Santa Marta, Cordillera de Merida, Cordillera Oriental) 
and probably predevonian (Sierra de Perija, Cordillera de Merida, 
Cordillera Oriental); it is quite possible that many of these schists are 
precambrian. 

4. The fact that W. KEHRER (22) has found in S. Colombia two 
different formations of schists, separated by an unconformity, indicates the 
possibility that there exist in both republics complexes of metamorphics 
of very different age. 

5. It is possible that the green schists of Antioquia) are cretaceous. 

6. It is very probable that in some parts of the Caribbean Mountains 
(s.1.) there exist metamorphics of cretaceous age. Every generalization, 
which would lead to the conception, that all the metamorphics of this 
area are to be regarded as Cretaceous, is to regard with the greatest 
suspicion. The granite of the Silla de Caracas, the data of HEDBERG for 
the NE. part of Venezuela and the data procured by the Northern Range 
of Trinidad may induce us to the greatest prudency. 

I cannot deny that the harvest of this critical study on the metamorphics 
of Colombia and Venezuela is somewhat meagre. We have seen that 
there exists the greatest uncertainty about the age of these 
rocks. It is, however, better to know this uncertainty than to feel satisfied 
with a not-founded knowledge. Happily, the remedy for gaining better 
knowledge can be indicated. Careful fieldwork, combined with intensive 
sampling and laboratory-study of the samples will with great probability 
lead to more definite knowledge of the age of the rocks under discussion. 
It is chiefly the green schists of Antioquia, the metamorphics of the 
Serrania del Interiér, parts of the Coast Range and the Northern Range 
of Trinidad, which promise good results. 
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Physiology. — Muscle sounds of a single twitch. By G. VAN RIJNBERK 
and H. D. BouMaAN. 


(Communicated at the meeting of February 24, 1940.) 


Introduction. 

HERROUN and YEO have described for the first time that a single twitch 
produces a muscle sound (1885). They stimulated the muscles of the upper 
arm with single shocs and listened with a stethoscope. 

With stimuli strong enough to give a clearly visible twitch, a low sound 
was heart comparable to the first heart sound. 

BERNSTEIN in 1890 found (also in listening experiments) that a single 
shoc given to the gastrocnemic of the rabit produces a muscle sound, It 
was of an explosive caracter comparable to p or ¢, or the sound of lighting 
a gas lamp. The frequency could not well be determined because of the 
extreme shortness of the sound. 


I. Experiments on dogs muscles. 


Our experiments were performed on dogs. Motor nerves were stimulated 
with single condenser discharges. The sound has been recorded by the 
same technique that has been used in our experiments on voluntary 
contractions in man. Dog and experimenter are enclosed in the sound proof 
room of the laboratory, the amplifier and recording string oscillograph are 
outside this room. 

We have used the gastrocnemic and some segments of the rectus 
abdominis muscle. In all experiments a muscle sound has been found. 
Experiments have been performed with the nerve intact and with the nerve 
cut. The results in both cases shows no differences. In some experiments 
the muscle (gastrocnemic) was left completely in situ, in other experiments 
the tendon has been cut and the muscle freed from its surroundings and 
placed on cotton wool. (The blood supply was not damaged by the 
operation.) The reason for this was the following. In some experiments on 
muscles in situ we were struck by the extreme regularity of the records, a 
regularity which was completely unknown in our previous records of 
voluntary contractions in man. The records had the caracter of being due 
to the tapping of some structure which can easily be set in vibration, The 
records in fact is those of an aftervibration. The question arose whether 
this vibrating structure is the muscle itself or the underlying bone. 

The answer to this question is obvious from two records reproduced here. 
One is from a muscle in situ while the other is from a muscle isolated as 


much as possible from the bone but with the blood supply undamaged. The 
muscle was placed on cotton wool. 
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The difference between the curves is clear. The vibrational caracter has 
entirely disappeared by the isolation of the muscle. Therefore these 
experiments are the first records of the true muscle sound of a twitch, 
HERROUN an YEO’s and BERNSTEIN’s results have been due to the under- 
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Muscle sound of a twitch Muscle sound of a twitch, muscle isolated 
Muscle in situ. with undamaged bloodsupply and placed 


on cottonwool. 


lying bones. A sound as recorded in fig. 1 can doubtlessly give an aural 
impression as described by BERNSTEIN. The sound recorded in fig. 2 
would sound entirely different, if at all audible. 

In all experiments the main feature of the muscle sound is a single 
diphasic effect: the af effect. In many curves this is followed by some 


Fig. 3. Rec. 84. 


The yd complex of the muscle sound of a single twitch of an isolated gastrocnemic. 
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irregular vibrations of low amplitude, which are not comparable to after 
vibrations. In some records where they are exceptionally big they have 
been labelled y 6 complex (see fig. 3). 

The curve does not change on cutting the motor nerve, the complex 
therefore cannot be described as a reflectory after discharge following 
the main twitch. As the nerve has always been stimulated close to the 
muscle, differences in conduction time of different nerve fibres cannot 
account for the after effect. 

An increase in frequency of the stimuli gives little or no change in the 
recorded sounds. (See figs. 4 and 5.) 

The phenomena recorded at higher frequencies when tetanic contraction 
sets in will be the subject of a seperate paper. 

The experiments show therefore that a single twitch does produce a 
muscle sound. If the muscle is in situ above the underlying bone the 
caracter of the sound is determined mainly by the bone. It is comparable 
to after vibrations. 

In an isolated muscle however the sound is much more simple. The 
main feature of the sound is a single diphasic vibration (af complex) 
followed in many experiments by a second diphasic vibration of much 
smaller amplitude (y 6 complex). It is impossible to discuss the pitch of 
such a phenomenon. 

We thought of finding a good muscle for these experiments in the 
rectus abdominis muscle of the dog. We thought it reasonable to suggest 
that if a middle segment of this muscle contracts, the surrounding segments 
will isolate it from its bony insertions (sternum, pelvis). The curve 
reproduced in fig. 6 shows that this is indeed possible. 
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Twitch of M3 of the rectus abdominis. 


The necessity for great care however becomes obvious from a comparison 
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of this record with the record reproduced in fig. 7. The contracting segment 
was situated one segment more caudad, and therefore closer to the pelvis. 
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Twitch of M4 rectus abdominis. 


In this curve the regular vibrations of an after vibration can already be 
recognised. 

The sharp excursion in this experiment preceding the complex is a 
signal which recorded the moment of the stimulus. 


II. Experiments on frogs. 


In order to verify the results obtained with dogs muscles, some 
experiments were performed on the more easily available gastrocnemic of 
the frog. Excised nerve muscle preparations have been used throughout. 

A simple apparatus was constructed which enabled us to keep the muscle 
under any tension necessary and at the same time eliminated vibrations 
not due to the muscle itself. It consists of a heavy lead bloc. On one 
end a support was attached (also of lead) in which the femur could be 
clamped. A small lead carriage moving on ball bearings and provided with 
another clamp served to secure the tendon. For isotonic contractions a cord 
attached to the carriage could be loaded with different weights. For 
isometric contractions the carriage could be fastened to the supporting 
lead bloc. 

Records of the acoustical phenomena of a single isometric twitch showed 
a constant picture, consisting in one single diphasic vibration. There is no 
after effect (y 6 complex). The results therefore appear to be in close 
agreement with those obtained on dog muscles. 

Yet we do not expect to continue our experiments with frogs as the 


much stronger and bigger dog’s muscles give much larger excursions in 
our records. 
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Conclusion. 


The acoustic phenomenon of the muscle twitch suffers from serious 
interference from sounds produced by the bones if the muscle is stimulated 
in situ. The interference appears in the form of a common type of damped 
vibration. If the muscle is isolated from its bony insertions and placed on 
cotton wool (blood supply undamaged), records of a much simpler type 
are obtained. 

The acoustical effect then consists of a single diphasic excursion 
(called af complex) followed by another diphasic excursion of much 
smaller amplitude (called yd6 complex). Records obtained from frog's 
muscles are not essentially different from those obtained from dog's 
muscles. 


Acoustics. — The residue, a new component in subjective sound analysis. 
By J. F. ScHouTen. (Natuurkundig Laboratorium der N.V. Philips’ 
Gloeilampenfabrieken, Eindhoven, Holland.) (Communicated by 
Prof. G. Hoist.) 


(Communicated at the meeting of February 24, 1940.) 


§ 1. Introduction. 


In a former paper 1) we described some experiments performed with an 
optic siren which permits to realize periodic sounds of any prescribed wave 
form. It was found that the non-linear distortion in the human ear is by 
far less pronounced than is stated by some authors and in particular that 
if a complex sound of moderate loudness objectively misses the funda- 
mental tone, this tone is not generated in any appreciable amount in the 
ear. It is a fact that a pitch equal to that of the fundamental tone is 
ascribed even to those sounds in which the fundamental tone is not present. 
The almost generally accepted hypothesis to account for that behaviour 
consists in the assumption that this fundamental tone is generated within 
the ear by means of non-linear distortion. Our experiments thus proved 
this hypothesis to be invalid. 

We were finally led to speculations as to how the ear might ascribe a 
pitch to a complex sound and suggested that the periodicity of the wave 
form rather than the distance of the harmonics in the Fourier spectrum 
might be the physical property determining this pitch. No mechanism 
presenting itself to enable the ear to perceive that periodicity, the question 
was left there unsolved. 

Further investigations and in particular the studying and repeating of 
experiments almost a century old2) led us to the conclusion that the 
fundamental problem underlying these and similar paradoxical phenomena 
is not a question of perception of pitch, but rather a question of subjective 
sound analysis. 

Once a radical change is made in OHM's seemingly trivial acoustical 
law of sound analysis’) the explanation of the “case of the missing 
fundamental” 4) and similar problems follows quite naturally, 


1) J. F. SCHOUTEN, The perception of subjective tones. Proc. Kon. Ned. Akad. v. 
Wetensch., Amsterdam, 41, 1086 (1938). 

*) A. SEEBECK, Beobachtungen iiber einige Bedingungen der Entstehung von Tone. 
Pogg. Ann. 53, 417 (1841). A. SEEBECK, Ueber die Sirene. Pogg. Ann. 60, 451 (1843). 

3) G. S. OuM, Ueber die Definition des Tones, usw. Pogg. Ann. 59, 513 (1843). 
G. S. OHM, Noch ein Paar Worte iiber die Definition des Tones. Pogg. Ann, 62, 1 
(1844). 

#) S. S. STEVENS and H. Davis, Hearing, New York 1938, p. 99. 
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This extension of OHM’'s law involves an important consequence both as 
regards our conceptions of subjective sound analysis, as well as regards 
those of the mechanism of sound perception, 

We shall, therefore, although the essential clues for solving the problem 


can be found on page 1092 of our former paper, reintroduce the matter 
from the very beginning. 


§ 2. Subjective analysis of a periodic impulse. 


Periodic sounds containing a great number of higher harmonics present 

themselves to the untrained ear as one sound of a certain sharp tone 
quality with a pitch equal to that of the fundamental tone. Since HELM- 
HOLTZ’ careful investigations 5) we know that the suitably trained ear is 
able to perceive the lowest dozen of harmonics with ease separately in the 
sound. This confirms OHM’s acoustical law which states that sinusoidal 
vibrations only are perceived as a pure tone, that a complex sound is 
analysed by the ear into its different sinusoidal components and that these 
components will be perceived as pure tones having a pitch determined by 
the respective frequencies, 
As a first experiment we listen to a periodic impulse of width = q 
which contains the lowest score of harmonics in slowly decreasing ampli- 
tude (Fig. la). The fundamental frequency used was 200 cycles/sec. Our 
conclusions to be made are naturally restricted in that respect. 

In this periodic impulse a strong and sharp note of pitch 200 is imme- 
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Fig. 1. Elimination of components, a. Periodic impulse of width a b. Periodic 


impulse without second harmonic. c. Periodic impulse without fundamental tone. 


5) H, HELMHOLTZ, Die Lehre von den Tonempfindungen, 1862, Kapitel IV. 
3) jf, 12, SCEOUWMAN, Ie. jos MOS 


Proc. Kon, Ned. Akad. v. Wetensch., Amsterdam, Vol. XLII, 1940. 24 
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diately perceived, whereas some higher harmonics, although by far weaker 
than this gravest note, can be heard separately one after another if the 
attention is fixed upon them. 

A very critical method of drawing the attention even of an untrained 
listener to a particular harmonic consists in adding first (by means of the 
second holder of the optic siren) that harmonic in the same amplitude but 
in opposite phase to the impulse. The harmonic, e.g. the second, (Fig. 16) 
is thus cancelled from the sound. By then covering the second holder the 
harmonic is made to reappear and can be heard in often surprising clear- 
ness. The loudness of the harmonics decreases rapidly with increasing 
order. The twelfth harmonic is about the last one which can be heard 
separately, the higher harmonics are not separately perceptible. 

So far the experiments, but for the great loudness of the gravest note 
(generally identified with the fundamental tone), present nothing essen-~- 
tially new. If, however, we now cancel the fundamental tone (Fig. 1c) 
we find that the sharp note of pitch 200 remains unchanged present in the 
perceived sound. Moreover, if thereupon the fundamental tone is again 
added to the sound, it is heard separately as a pure tone of pitch 200 of 
low loudness comparable to that of the second and third harmonic. After 
some training this fundamental tone may even be heard without any help, 
although with more difficulty than the next harmonics. 

The crucial point is thus that, as to subjective analysis, the sound 


O% contains two components of pitch 200, 
140 one of which, having a pure tone-quality 
tae is identical with the fundamental tone, 
100 whereas the other, having a sharp tone 
80 quality and great loudness, is of different 
60 origin. 

Kh Measurements of the loudness of the 
se various harmonics were taken by com- 
; paring in successive contrast a particular 


RoE Se 45 6 7 8 harmonic in the periodic impulse and a 
lumber of harmonic 

a, pure tone of the same frequency. In 

eli emea ae cenere Fig. 2 the relative amplitudes of the test 

Genres, Cees 2, WA\. frequencies, necessary for the match, are 

given. It will be seen that the funda- 
mental tone in the periodic impulse is heard slightly too strong, the second 
harmonic in almost its true loudness and the higher harmonics gradually 
weaker and weaker, 

If the harmonics are eliminated one by one, starting from the lowest, 
the sharp note does not, at first, materially change either in character or 
in loudness. If, however, the same is done starting from the highest the 
sharp note gradually loses in sharpness as well as in loudness. This 
behaviour suggests that the sharp note is connected with the presence of 
high harmonics in the complex sound. 
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If the intensity is raised, the pitch of the sharp note remains practically 
unchanged. The fundamental tone, however, exhibits the well-known fall 
up to about half a tone. 

The totally different nature of the fundamental tone and the sharp note 
of pitch 200 is best brought out by adding the fundamental tone in 
increasing amplitude to the impulse. The increasing loudness of the funda- 
mental tone is easily perceived, whereas the sharp note does not materially 
change in character and merely seems slightly to decrease in loudness. 

Thus the fact that the lowest note in a complex sound is more easily 
perceptible than e.g. the second, third and fourth harmonic is not due to 
a particular enhancement of the loudness of the fundamental tone or to a 
general impairment of the harmonics in question, but to the presence of 
a hitherto unknown additional subjective component 7) of almost the same 
pitch as the fundamental tone. 


§ 3. The hypothesis of the residue. 


Hitherto, as expressed in OHM’s law, it was generally accepted that the 
ear analyses a complex sound into components of a pure tone-quality each 
of which corresponds with one frequency of the inner-ear sound field 8). 
The difficulty remained that the highest harmonics, although not separately 
present in subjective analysis, add materially to the loudness as well as 
to the tone-quality of the complex sound. 

We now find that, apart from those pure components, an additional 
component of sharp tone-quality may exist which cannot be correlated with 
any single frequency of the sound field. 

We propose to call such an additional subjective component a ‘‘residue’”’. 
It is very well possible that in a complex sound several residues are present. 

Oum's law of subjective sound-analysis may now be extended as follows: 

1. The ear analyses a complex sound into a number of components 
each of which is separately perceptible. 

2. A number of these components corresponds with the sinusoidal 
oscillations present in the inner-ear sound field. These components have a 
pure tone-quality. 

3. Moreover, one or more components may be perceived which do not 
correspond with any individual sinusoidal oscillation, but which are a 
collective manifestation of some of those oscillations which are not or 
scarcely individually perceptible. These components (residues) have an 
impure, sharp tone-quality. 


§ 4. The pitch of the residue. 


The loudness of the residue in a periodic impulse is greatly diminished 


7) The alternative hypothesis, consisting in the supposition that the observed sharp 
note is not a component of the sound, but the total sound itself leads to a great many 
difficulties and can, we believe, be discarded. 

8) Thus including the components generated within the ear by non-linear distortion. 
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if the higher harmonics are cut off, e.g. by means of a low-pass filter. 
The residue is thus a collective manifestation of those higher harmonics. 

Which physical property of these harmonics might determine the pitch 
of the residue? Two possibilities present themselves: either the distance 
between the harmonics, or the periodicity of the total wave form of the 
harmonics in question. In the case of the periodic impulse the two properties 
lead to the same pitch, the distance between the harmonics and the periodi- 
city both being 200 cycles/sec. 

Comparison of the two wave forms represented in Fig. 3 enables us to 
answer this question. The first wave form contains the even harmonics 


WAVE FORM SPECTRUM 


a 


24 6 8 10 12 14 16 18 20 
Number of harmonic 


Ube 
13579 18 15 17 19 


Fig. 3. Two wave forms having the same distance between the harmonics 
but a different periodicity. 


only and is thus the octave of the periodic impulse used hitherto. The 
distance between the harmonics is 400 cycles/sec., the periodicity of the 
wave form, as well as of any group of adjoining harmonics is also 400 
cycles/sec. The pitch of the residue is found to be 400. The second wave 
form contains the odd harmonics only. In this case the distance between 
the harmonics is again 400 cycles/sec., the periodicity, however, be it of 
the total wave form or of any particular group of adjoining harmonics is 
200 cycles/sec. We found that a residue is present, although less pronounced 
than for the first wave form, and that the pitch of this residue is 200. No 
indication of any subjective component of pitch 400 was found. 

The ear thus ascribes a pitch to a residue by virtue of the periodicity 
of the total wave form of the harmonics which are responsible for this 
residue. 

Summarizing we are thus led to the following view upon subjective 
analysis of a periodic sound containing a great number of harmonics. The 
lower harmonics can be perceived individually and have almost the same 
pitch as when sounded separately. The higher harmonics, however, cannot 
be perceived separately but are perceived collectively as one component 
(the residue) with a pitch determined by the periodicity of the collective 
wave form, which is equal to that of the fundamental tone. We are thus 
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confronted with the surprising result that the harmonics highest in frequ- 
ency, are perceived as a subjective component almost lowest in pitch. In a 
periodic impulse the residue must be the most prominent component due 
to its richness in harmonics and to the comparatively high sensitivity of 
the ear to these harmonics. 

The pitch of a complex sound now follows quite naturally: The pitch 
ascribed to a complex sound is the pitch of that component to which the 
attention, either by virtue of its loudness, or of its contrast with former 
sounds is strongest drawn. 

Therefore the pitch of a complex sound may be different depending upon 
the circumstances under which it is heard. An example of this behaviour 
was given on page 1092 of our former paper. 


§ 5. SEEBECK’s experiments. 


It should be definitely stated that the views proposed here, although 
rather radical in comparison with the modern acoustical point of view, are 
no more than an extension of the work of the admirable and ununderstood 
SEEBECK. In 1841 SEEBECK (l.c.) described a great number of experiments 
performed with his acoustic siren. This lead OumM (l.c.) in 1843 to the 
formulation of his famous law. Far from being satisfied SEEBECK in the 
same year published a second paper in which he proved that OuM’s law, 
although qualitatively accounting for the observed phenomena, could by 
far not describe these quantitatively. Therefore SEEBECK proposed an 
extension of OHM’s law. The discussion is closed by OHM in 1844 with a 
second paper, which, though humorous, added nothing essential to clear 
the difficulties. 

HELMHOLTZ (l.c.) too, while discussing the controversy between 
SEEBECK and OHM impresses the opinion upon the reader that SEEBECK, 
though otherwise a very keen observer, must have been mistaken here. 

We shall describe three of SEEBECK’s experiments. In the first experi- 
ment air was blown from a pipe against a revolving disk containing a 
number of concentrically arranged equidistant holes. SEEBECK thus produced 
a ‘periodic impulse’. If a second pipe was placed against the same side of 
the disk at half-interval distance from the first the pitch was heard to 
jump an octave upwards (Figs. 5a and 5b). The number of impulses per 
second is thus doubled, the spectrum will contain the even harmonics in 
double amplitude and no odd harmonics. SEEBECK’s objection to this 
formulation as a means of quantitatively describing the observed pheno- 
mena is the following: in the first sound the loudness of the gravest note 
is very great compared to that of the second and third harmonic. If the 
odd harmonics are compensated and the even ones doubled, the second 
harmonic (now the fundamental tone of the new sound) should be a little 
stronger than in the first sound but not so strikingly strong as is actually 


heard. 
HELMHOLTZ suggests that SEEBECK, while not using the proper means 
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of drawing his attention to the second harmonic of the first sound, may 
have underrated its loudness. The experiment, however, is easily repeated 
with our optic siren and completely confirms SEEBECK’s description. In 
view of the hypothesis of the residue the explanation is simple, The strong 
gravest note heard by SEEBECK is not the fundamental tone (which must 
have been of almost the same loudness as the second harmonic) but the 
residue. Doubling the frequency of the impulse not only doubles the 
amplitude of the second harmonic but also shifts the residue an octave 
upwards, thereby seemingly increasing the loudness of the second harmonic. 

SEEBECK’s second experiment is more immediately convincing. He 
compares three wave forms (Fig. 4). The first (Fig. 4a) is again obtained 
by one pipe. The second (Fig. 4b) by two pipes at a distance equal to the 


WAVE FORM SPECTRUM 


2 


12345 10 15 20 
Number of harmonic 


el ey) FP TIPE SIE i) 

Fig. 4. SEEBECK’s experiment 2. The amplitude of the fundamental tone in 

b and c is twice that in a. Yet as to subjective analysis the gravest note in c 
is barely stronger than in a and strikingly weaker than in b. 


interval of the holes. This merely results in doubling the amplitude of the 
periodic impulse. The third wave form (Fig. 4c) is obtained by placing 
the second pipe against the opposite side of the disk at half-interval 
distance from the first pipe. 

It is easily verified that the amplitude of the fundamental tone (which 
has the same pitch for the three wave forms) in the second and third wave 
form is twice that in the first one. 

Contrary to OHM’s expectation SEEBECK finds that in the third wave 
form the gravest note is barely stronger than in the first and strikingly 
weaker than in the second wave form. 

Again the explanation is simple once it is realized that SEEBECK did not 
observe the fundamental tone, which must have been comparatively weak, 
but the residue. This residue must be much weaker in the third than in 
the second wave form since the former contains only half the number of 
frequencies, 


Lastly we describe a third very elegant experiment. If, in the case of 
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two pipes placed against the same side at half-interval distance, this 
distance is slightly changed, the pitch of the sound immediately jumps an 
octave downwards. For demonstration of this effect SEEBECK constructed 
a disk containing four concentric rows of holes. In the first row the 
distance between the holes was 20°, in the second 10°, in the third alter- 
nately 94° and 104° and in the fourth 9° and 11°. 

The second row, of course, gave a pitch an octave above that of the 
first row, in the third the octave was still the most prominent although the 
lower note was distinctly audible, in the fourth this lower note was more 
prominent than the octave. It is seen from Fig. 5 that even for the last 
wave form the fundamental tone is still very weak. 

We repeated this experiment by putting one impulse in the first, another 


WAVE FORM SPECTRUM 
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Fig. 5. SEEBECK’s experiment 3. The pitch of wave form 6 (400) is an octave 

above that of wave form a (200). In wave form c a component of pitch 200 is 

distinctly audible. In wave form d it is even more prominent than the component 

of pitch 400, The fundamental tone, however, is very weak in wave form c as 
well as in wave form d. 


in the second (turnable) holder of the optic siren. If the two impulses 
were exactly half a period apart the pitch was 400. A slight turning of 
the second holder, however, was sufficient to make the pitch drop to 200. 
The effect is very striking indeed. Using again the impulse of width = 
we observed the same behaviour as that described by SEEBECK. 

In terms of the residue the explanation of this effect is that this slight 
change causes especially the higher odd harmonics to reappear which 
results in the change of pitch of the residue from 400 to 200. The funda-~ 
mental tone itself does not play any part in the phenomenon and, in fact, 
can be heard very weakly after a slight change and more strongly after 
a great change in the distance between the two impulses, quite in accor- 
dance with its objective intensity. 
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To account for these phenomena, all of which show a discrepancy 
between the theoretical and experimental loudness of the gravest note in 
subjective analysis, SEEBECK suggested that those higher harmonics, which 
cannot be perceived separately and which have a common period equal to 
that of the fundamental tone, in some way or other enhance the loudness 
of this fundamental tone. He even went so far as to suggest that this 
enhancement is not merely due to non-linear distortion in the ear (l.c. 
1843, p. 480). 

The essential feature of our investigations consists in stating that it is 
not the fundamental tone itself which is so much stronger than should be 
expected, but that an additional subjective component of almost identical 
pitch and of often great loudness is present in the sound. 

It is OHM's merit to have indicated the general principle underlying 
the phenomena of subjective tone analysis and to have formulated a law 
which up till now was considered to be so trivial that a renewed critical 
testing seemed superfluous. 

We cannot but immensely admire SEEBECK to have realized, by means 
so simple as his acoustic siren the short-comings of OHM’s law and even 
to have suggested an extension of this law which was scarcely considered 
seriously until now, almost a century later. 

In a following paper we hope to investigate the theoretical consequences 
of the existence of the residue and to show that, as regards the physical 
part of the analysing mechanism in the inner-ear, the existence of this 


new component is by no means so improbable as one might expect 
beforehand. 


§ 6. Additional remarks, 


It is of interest to reconsider now some other older and newer experi- 
ments on subjective sound analysis. 

KONIG 9), experimenting with SAvaRT’s siren, found that while holding 
a stiff wooden peg against the cogwheel (128 teeth, time of revolution 
1 sec.) a rattling sound as well as the tone c (128 cycles/sec.) are heard. 
If however, the edge of a paper card is taken, the rattling sound is scarcely 
perceptible whereas the tone c can be heard very clearly. 

Our interpretation of this behaviour is that the first sound contains a 
great many harmonics and thus brings about a residue (the rattling sound) 
whereas the second sound is comparatively pure. 

A highly remarkable consideration can be found in STUMPF’s work. 
STUMPF 1°) experimenting on synthetic vowels obtained by sounding 
together various harmonics of a given frequency finds that the pitch is 
equal to that of the fundamental tone even if that tone is not objectively 
present. He ascribes this throughout to a difference tone generated within 


®) R. KONIG, Ueber den Zusammenklang zweier Tone. Pogg. Ann. 157, 226 (1876). 
10) C. STUMPF, Die Sprachlaute, Berlin 1926, p. 185. 
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the ear by non-linear distortion. At one place, however, he wonders: “Yet 
I have doubted sometimes whether it (the difference tone) is, even sub- 
jectively, really present within the ear. One might imagine that in a sound 
consisting of the objective tones c2g2c3e3, the pitch cl only presents itself 
to the listener, without the fundamental tone c1 entering subjectively into 
the sound”, 

This paradoxical formulation is very similar to our interpretation in 
our former paper, which we now, however, prefer to abandon in favour 
of the supposition that an alien component, the residue, is present in the 
sound, This component may, if loud enough, determine the pitch of the 
total sound. 

The hypothesis of the residue may also be of importance in connection 
with the strike note of bells. We hope to deal more explicitly with that 
fascinating problem later. It may suffice now to draw attention to its chief 
characteristics. 

In well-designed church bells a prominent note: the strike note, is heard. 
No partial corresponding to the pitch of that note is present in the unhar- 
monic spectrum. An experimental rule, still adhered to by those workers 
not sufficiently misled by its theoretical improbability, asserts that the pitch 
of the strike note is an octave below that of the 5th partial of the bell 
spectrum. More physical or technical minded people, however, having 
found that the presence of the strike note is linked up with the presence 
of at least the 5th and the 7th partial, and having stated that the difference 
in frequency of those partials is sometimes almost equal to the pitch of the 
strike note, interpret it as a difference tone (formed within the ear) of 
those partials. In its essential features we are confronted with the “case of 
the missing fundamental” all over again. 

On account of our hypothesis of the residue one might expect the 
following behaviour: if, in an unharmonic spectrum, by chance a number 
of partials are integer multiples of a certain frequency, an additional sub- 
jective component of corresponding pitch having a sharp tone quality and 
eventually great loudness should occur. If no such grouping is possible 
no extra component of definite pitch should be present; if several groupings 
are possible more than one additional component might occur. All of these 
expectations correspond to essential features of the strike note of bells. 

We are thus led to suppose that the strike note of bells should, indeed, 
not be interpreted as a difference tone and moreover that it is another 
example of a residue. Conversely one might call the residue the ‘strike 
note’ of periodic sounds. 


The author wishes to thank Ir. R. VERMEULEN for his valuable criticism. 


S.i.H., January 1940. 


Mathematics. — Ueber BeESSELsche, STRUVEsche und LOMMELsche 
Funktionen. (Zweite Mitteilung). Von C. S. MEIER. (Communi- 
cated by Prof. J. G. VAN DER CORPUT). 


(Communicated at the meeting of February 24, 1940.) 
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§ 4. Aus (36), (55) und (56) geht hervor 


K, (2z)=42 Go? (z2|0, -) = 4 2 Gis 


Ebenso folgt aus (35), (55) und (56) 


Jy (22\=2 Gis (= 


und 
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Wendet man nun (26) auf die rechte Seite von (74) an, so erhalt man 


UG, (2 z) = roy! ot ( u? 
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Nimmt man hierin a=v¥, so bekommt man mit Riicksicht auf (49) 7°) 


Ki vad Pee ere 28 (yes 


ae Tice 


diese Beziehung gilt fiir |argz| << 47 und R(r) > 4. 
Auf analoge Weise findet man mit Hilfe von (75) und (45) 


Deve: i 
J (22) — FES | J, (zu) Y, (zu) (u2—1)'-! u'-2" du; . (78) 
1 


hierin ist z>0 und R(v) >4. 
Aus (76) und (46) folgt schliesslich 


22) sin yt ————— 


(—y— 4) 


ef Meth JF (z u)} (u2—1)-"-! u!+?" du; (79) 


hierin ist z > 0 und R(r) << — 2, 


Afyeeey ap\. . : 
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§ 5. Fir die Funktion Y,(z) gilt wegen (26) und (37) 
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aus (83) mit a—+»7-+A2 ergibt sich also 
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Diese Beziehung gilt fiir z > 0, R(v) >—# und beliebige Werte von 
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Die den Relationen (81) bezw. (82) entsprechenden Integraldarstellungen 
der Funktion J,(z) sind 
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Te J_, (z) + Y, (z) sin vx 
cos v7 
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auf analoge Weise folgt (88) aus (58) und (82). 


27) Siehe (32). 
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Die Relationen (81), (82), (84), (85), (86), (87) und (88) waren 
bekannt ”°), 


§ 6. Die Anwendung von (7) auf (42) liefert 


Qu 1 av 
Ss »(z)= EF ieee 2s a gy | 
ba! resecaTianrera) OCH ‘ ince) (89) 
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Die Funktion S,,,(z) besitzt daher wegen (89) (mit a=314-+ 406 und 
B=+4A—H4o) die Integraldarstellung 


22 Gee ie 3 bEA+ $e) (° | 
ty (Za Ze Lu—dty yo 2 ) [ss (zu) 
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Hierin wird z~0, |argz|< 2 und R(w+y) <1 vorausgesetzt; 4 und 
o sind beliebige Zahlen mit (4) >O und u—A +o F1,3,5,.... 

Fiir oc—yv—ZJ findet man eine den Formeln (58), (63) und (84) ent- 
sprechende Beziehung, namlich 
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2 2 


diese Beziehung gilt fiir |Jargz| <2, R(w+y) <1 und beliebige Werte 
von A mit 8 (4) > 0. 

Man bekommt eine mit (69) verwandte Integraldarstellung fiir Syu,»(z), 
een neman te te und w= cosif setzt in (90). Denn durch diese 
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28) NIELSEN, [16], 222—223; WATSON, [17], § 6.13, Formel (4). Fir (86) vergl. 
man Mayr, [5], 230, Formel (3a). 
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liber: in dieser Relation ist |argz|< 2, R(m+y»)<1 und o beliebig 
mit R(o)< 1 und }u+oFl. 

Setzt man o——4 und u=cosht in (90), so erhalt man mit Riicksicht 
auf (73) (ich benutze S,,»= Sy,-») 


oh () =e ee 
wr Psu +9) 


Pe (93) 
x | S,,-2,3 (Z cosh #) Dee (cosh £) sinh’ t cosh! t dt; | 
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hierin ist jargz|<a, R(wtr)<1 und 4 beliebig mit RAZ) >0 und 


an ie 
Der Spezialfall von (92) mit o=4 liefert wegen (70) 


Qu AT a 
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diese Beziehung, giiltig fiir |argz|< 2 und R(u+v)<1, kommt auch 
zum Vorschein, wenn man A=# setzt in (93). 


§ 7. Aus (41) und (42) geht hervor 


H. (2) — ¥. (@) == yeaa on a er) 


Mit Hilfe dieser Beziehung kann man aus (90), (91), (92) und (93) 
Integraldarstellungen fiir H,(z)— Y,(z) ableiten. Aus (94) und (91) 
folgt z.B. 
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hierin ist |arg z| <a, R(v) <4 und 2 beliebig mit KR (2) > 0. 


Ist | arg z|< a und R(a—4)>4, so findet man durch Anwendung 
von (26) auf (41) 
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Nun hat man wegen (55) und (50) 
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aus (96) mit a=» ergibt sich somit, falls |argz|<a und K(v) 4 ist, 
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Ebenso folgt aus (39) und (48) 
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diese Beziehung gilt fiir |argz| << $2 und R(v) >t. 
Schliesslich folgt aus (40) und (47) 


{J_,(2z)—L, (2z)} sin va = 


ah 


hierin ist |argz|<($a und R(v) > 4. 
Die Relationen (97), (98), (9 
(79); Formel (98) war bekannt ?’). 
Mit Hilfe von (26) kann man auch Umkehrformeln *°) fiir (97), (98), 


(99) und (100) ableiten. Die Umkehrformel von (98) lautet z.B. 


Sa s F 
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hierin wird z >0 und | (v)| <4 vorausgesetzt. Diese Beziehung kommt 
zum Vorschein, wenn man (26) mit a=+4 auf (44) anwendet und (38) 
benutzt. 
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und (100) entsprechen (77), (78) und 


29) WATSON, [17], 417, Formel (9). 
30) Die Umkehrformel von (26) 1st 
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Aus (26) und (40) geht hervor 
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2 B 
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ersetzt man nun a durch A—4¥, so erhadlt man wegen (40) die mit (58), 
(63), (84) und (95) verwandte Formel 


{L(2}Ls(2)} cos (>a) x= = — SS i) (1 yaa (ze) —Lya(za)} (2 — 1)! v4 du; 


hierin ist |argz|< $2, R(v) <4 und 4 beliebig mit Ht (A) > 0. 
Die entsprechenden Integraldarstellungen der Funktionen H,(z) und 


I,(z)—L,(z) sind 
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(wo K(v) >— 2 und J beliebig mit O< 8 (2) << RH (v) + §) und 
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(wo Ki(v) >—1 und 4 beliebig mit O< H(A) << KR (r) + 1). 

Der Beweis von (101) bezw. (102) geht durch Entwicklung der Funk- 
tionen H,_,(zu) bezw. J,-,(zu)—L,-2(zu) (siehe (31) und (34)) und 
gliedweise Integration *'). 

Nun folgt aus (34) und (31) 
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31) Die (101) und (102) entsprechende Integraldarstellung fir J» (z), namlich 
Gjiaah 2 : 
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ist schon lange bekannt; man vergl. WATSON, [17], 373, Formel (1). 
32) Ich nehme an, dass R (v) > — + ist. 
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Diese Beziehungen waren bekannt 33). 


Ibs (z) = L, (z) == 


§ 8. Durch Anwendung von (7) auf (49) erhalt man 


Nie. is 7 )/ 
Ae le ae 
Re nd TE eae | . (103) 


X oF, (a—v, a + v3; a+ B; 1 —u?) (u2@—1)"+ #1 y!-28 dy, 
Nun folgt aus (56), (55) und (36) 


l= Gos (C2| 8, 0) = 20? Ke (20). 


Setzt man a=+4 und u=cosht in (103), so findet man also mit 


Riicksicht auf (68) 
K> (z)=228 Va | Kg (2z cosh ¢) Piss (cosh 2 t) sinh®+#tcosh?tdt; (104) 
0 


hierin ist |argz|<(42 und f beliebig mit K(f) >— 4. 
Mit Hilfe von (46) und (35) beweist man ganz analog 


8 5 
SEA Ga te as Ja(2zcosht) P =? (cosh 22) sinh’ +? ¢ cosh? t dt; 


diese Beziehung gilt fiir z > 0, |W (v)|< 1 und beliebige Werte von f 
mit —2< RK (6) << —2|R()| +3. 
Die Anwendung von (7) auf (50) ergibt 


a) Q) 4 cos vat z | 
A,” (z) Hy” (z) = Gana rete { 8 uw el . (105) 


x< oF, (a—v,a +»; a + B; 1—n?) (u2@—1)**F 1 a7? da. \ 


Nun hat man wegen (55) und (41) 


4—48 ) 
4—48,48,—48 


33) Man vergl. WATSON, [17], 330, Formeln (1) und (2). 


Proc, Kon, Ned Akad, v. Wetensch., Amsterdam, Vol, XLII, 1940. 95 
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Aus (105) mit a=4 und u=cosht folgt also mit Riicksicht auf (68) 


H°” (z) Hy (z eos tees fi {H_3(2zcosh t)— Y_(2 zcosh#)} 


< P!=? (cosh 2#) sinh?*? ¢ cosh! ¢ dt; 


2 


in dieser Beziehung ist |argz|<2, |R(%)|<4 und f beliebig mit 
(0) Bates 
In ahnlicher Weise findet man mit Hilfe von (47) und (40) 


22° sin 2 


ln af Ig(2z cosh t)—L_g (2 z cosh #)} 


(0) 
>< P!—# (cosh 2 £) sinh?+! t cosh? t dt; 


{1?,(z)\—I; (z)} cos z= 


in dieser Beziehung wird |argz|<4$2, |R(»)|< 4% und R(~) > —Z 
vorausgesetzt. 
Eine verwandte Integraldarstellung fiir J,(z) J_,(z) ist 


_ 2 B pi! A 
ACA Fst a rae (2 zcos¢) )Pi_f(cos2¢)sin®+! pcos! y dy; (106) 


hierin ist PB beliebig mit —$< K(f) <1. 
Denn die zugeordnete LEGENDREsche Funktion erster Art wird fiir 
—1<w<1 definiert durch **) 


m We iether 
Ea le) — I'(1—m) 


oF, (—n, 1 +n; 1—m;4—4w). 


Die rechte Seite von (106) ist also mit Riicksicht auf (31) gleich *°) 


1 
1 
AIGA ag | OP et 2 eb 1a) (a) ®t? 
6 


BS pity 


Ist —4£< K(f) <1, so besitzt dieser Ausdruck den Wert 3°) 


1 
P(i—v) F(t +») 1F, (4; l—», 1 + »; —2?), 


und diese Funktion ist gleich *”) J,(z)J_,(z), so dass Formel (106) 


bewiesen ist. 


34) Man vergl. HOBSON [2], 227. 


) 
35) Ich ersetze cos? g durch u. 
36) Siehe [13], Formel (16). 

37) WATSON, [17], 147, Formel (7). 
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§ 9. Fir das Produkt K,.(z) K,(z) gilt wegen (7) und (52) 


IGZ)tG (2) = rant Ge (2 Fe 


MaF (a—F eter atsu—ty;a+B; 1—v?)(u2—1)*+8-1 y1-28 du, 


One )| 
a,f.tu+4r,—hu—$r) | 107) 


Nun folgt aus (56) und (36) 


= Go2 (0 tu +4», —tu—4r) = 2Kys, (20): 


weiter hat man *°) 


De AWA py eee cinlie 1) -—="cosh 2e. 8 ee ye (108) 
Aus (107) mit a=0, 6=4 und u=cosht ergibt sich also 
Ka Kh = if Ga Or eaniicsineniede . | 610e) 


diese Beziehung, giiltig fiir |argz|< 42, war schon lange bekannt 29). 
Durch analoge Anwendung von (7) auf (51) bekommt man 


Tul | Jez cost) cosa (uy) tt; (110) 


hierin ist z>>0O und |R(u—»)| << 2. 
Mit Hilfe von (7) und (54) findet man fiir die Funktion V,,,,(z) 


oy 18) 
Q* 
“| a, B, rH oY, —% fe 2) (111) 


ae (4 —ty) 


oat (Gt a Or ee 4s a Bs Pn?) (1) ae dia, 


Nun hat man wegen (56) und (42) 


il 
t ) 
14, 
tute, —tu—F? 


1 
Dy , 
2 T(E -Fp— =”) LS ap oa aa a) Seon. (20. 


pole 


38) Man vergl. [6], Formel (42). 
39) WATSON, [17], 440. Siehe auch MEIJER, [8], 16—17; [9], 483. 
2 


376 


Setzt man a—=0, B=+ und a=cosht in (111), so erhalt man also mit 


Riicksicht auf (108) 


a? 


ae rp ey 2 
iW = Sco iis | So,u+»(2z cosh #) cosh (u—) ¢ dt. 
0 


< 1, ist mit (109) 


Diese Beziehung, giiltig fiir |argz| <a und | R(u—7) 
und (110) verwandt. 

Durch Anwendung von (7) auf (53) findet man die entsprechende 
Integraldarstellung der Funktion U,,,(z), namlich 


ML 4 


fey pani meee 
Use. (ies EE +H) | Suen (22 Cosh t)cosh (i) t de; 
0 


hierin ist |argz|<a und |R(u—»)| <2. 


§ 10. Die Anwendung von (7) auf (54) liefert nicht nur (111), sondern 
auch 


__ 4cos($u+4)acos(t u—ty)a - 
Weyl eee u ies ae : [oe (2 


X oF, (a—4u—tr,a—4u+ 4»; a+ B—p; 1—v?) (w?—1)4 72-1 8 du. 


Nun gilt wegen (56), (55) und (42) 


Sie Ge 2 e+ee = i-H ee (1 — 4 ; ya ; 
133 eG sat ey (—tu—4r)T($—Fu+ 47) Sy» (20). 


Aus (112) mit a=0, 6=+4 und u=cosh ¢ ergibt sich also mit Riicksicht 
auf (73) 


V. (z) ey 2-24 +s Za 
¥, co — 
PEteet+a)lG+zu— ya 


al Sy,» (2z cosh t) P*} (cosh t) sinh'-“ t cosh—" t dt; 
0 


hierin ist |argz|< 2, R(u) <4 und R(utr)>—1. 


Kul Kole) = aya ia. mf (2 cotgh? t 
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Auf analoge Weise findet man mit Hilfe von (53) und (42) 
DATE Z aa 


U,, »(z) = 
Deets) Pgu—4n) a 


<f Su-1,» (2z cosh t) P!*? (cosh t) sinh?~“ t cosh—¥ ¢ dt; 


diese Beziehung gilt fiir |argz| <2, R(u) <1 und R (utv) > — 2. 
Die (112) entsprechende Integraldarstellung der Funktion K,,(z) K, (z) 
ist wegen (7) und (52) 


Xo (a—gu—hy, a— gut dy; a+ $—; 1—uv?) (w2@— 1) 8-11 pe 


Nun hat man infolge (56), (55) und (36) 
G3a(¢ 


Aus (113) mit a=0, B=+ und u=cosht geht also hervor mit Riick- 
sicht auf (73) 


0, 4 
. = Goo 2 i 1 1 1 
L 5 , 
OS ead aie ae sv 2(C 2M, g | 5) 


So" Go 2K? he, he) = 20 Ko (2). 


TGAZ) G2) == 2* va fm 2. zcosh t) P*i (cosh ¢) sinh!" t cosh— t df ; 


diese Beziehung, giiltig fiir |argz| << $2” und R(u) <4, habe ich friiher *) 
auf andre Weise abgeleitet. 
Ebenso beweist man, et von (51), 


Da= Pz —lU 
iVa 


hierin ist z > 0, R(u) <4 und R(2u+r) > —3. 
Setzt man a—4—tu-+7, p= —tp 
(113), so erhélt man 


Lez) J, (2. z cosh t) P!*} (cosh £) sinh!-" tcosh—" t dt; 


i 


< oF, (4—4,4—uty; 1—2u; 1—cotgh? t) sinh?“-”~? t cosh“*” ¢ de. 


40) MEYER, [9], 482, Formel (51). 
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Nun ist *’) 
Ff, (4—m, 1—m--n; 1—2 m; 1—tgh?¢) 
= 2-2m '({—m) sinh 7 ¢ cosh? 7" t Px (cotgh 2 t); 


diese Beziehung geht mit Hilfe von (13) in 
oF, (L—m,1—m + n; 1—2m; 1—cotgh? £) 
—2-?m ['(1—m) sinh"-?"*1 tcosh“"“! t P7 (cotgh 2 t) 


iiber. 
Fiir die in (114) vorkommende Funktion ,F; gilt daher 


oF (4—u,4—p+; 1—2 4; 1—cotgh? ¢) / a 
: 1X) 
— 2-# P'(1—p) sinh’-2"+4 t cosh~’? ¢ Py", (cotgh 21). ) 
Ferner hat man wegen (55) und (43) 
0, + 
4,0 2) + 2 
osi(¢ ied ae doe ae 


(116) 


Oe 
~~ 
oe 

| 
tol 


=e Gos (© 


eel eet 1 
Bort ice sven iar oa ies 2 


Jartemiyaet Wi 20 


Die Funktion K,.(z) K,(z) besitzt also mit Riicksicht auf (114), (116) 
und (115) die Integraldarstellung 


Qu-4 ov -t at 
Kulz K, Ne eZ cotgh t 
12) K(= Tay | 


Xx Wien, (22 cotgh t) P.-4 (cotgh 24)snh = feosh> cat: 


hierin ist |argz|< $a und KR (u) <4. 
Genau so findet man, wenn man a= —}—tu+4y, fp=— 
und u=cotght setzt in (113), 


+i ,—u—-i 3 
K, (z) Ky (= = = s 3 {face 
2 


0 


* W_ 3-4,» (2z cotgh t) BP’: (cotgh 2) sinh?” t dt; 
diese Beziehung gilt fiir |argz| <4 und R(u)< —i. 


41) Man vergl. HOBSON, [2], 236, Formel (74). 


Mathematics. — Développements en série de polynomes d' HERMITE et 
de LAGUERRE 4 l'aide des transformations de GAUSS et de HANKEL. 
Iif *°). Par ERvIN FELDHEIM. (Communicated by Prof. J. G. vAN 
DER CORPUT). 


(Communicated at the meeting of February 24, 1940.) 


§ 5. Quelques généralisations ultérieures. 


Cette partie du travail contient des généralisations pour certains 
résultats de I et II, et surtout pour la formule (21) qui a été établie 
dans un travail actuellement sous presse *’”). Nous allons montrer aussi 
que les développements (18), et une inversion de (60), permettent de 
donner une nouvelle interprétation pour quelques résultats connus. 

1°. Généralisation de l’équation intégrale (21). Si nous remplacons, 
dans la formule (46), y et z respectivement par u+v et u—v, cette 
relation nous fournit la transformée de GAUSS Gx(a,u) du produit 


ite a | Edn a u et v étant des valeurs quelconques de I'inter- 


valle (—o,-+ 0), les autres paramétres restant assujettis aux conditions 
imposées pour (46). Si nous posons, dans la relation obtenue, a=A=uw=1, 
nous serons conduits a l’équation intégrale 


+ 0 
—= tf pe (bea, 9) ce = | 
aU . 


—o 


(61) 
—=2"Al(vo+iu™—" LY (2u2+2v2) (m=n) \ 


Cette équation résulte immédiatement de (20), et ne peut pas étre consi- 
dérée par suite comme une généralisation de celle-ci. Mais son inversion: 


+0 
annie fect (o-Lin)?—? 1” (20? +20?) du= 
aU 


=H, (x--0) p(x —v) 


= 
IV 


26) Les parties I et II ont été publiées dans les Proc. Kon. Ned. Akad. v. Wetensch., 


Amsterdam, 43, 224—239; 240—248 (1940). 
27) FELDHEIM, [8]. Les nombres entre crochets renvoient a la Bibliographie placée 


a la fin de la 2© communication. 
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constitue bien une généralisation de la relation (21), qui ne peut pas 
atre déduite du développement (16). Nous allons donner un peu plus 
tard le développement équivalent a (62), mais démontrons d’abord cette 
équation intégrale a l'aide de (61). On en tire 


+ 0 


2m n! Ls J ecteroe (v+iu)™-" LY 40 = oid 
Tt e 


+00 +00 
= mall (ees eak da | ecto Hn (y tv) Ha (y —v) dy. 
IU 


—o 


Si l'intégration par rapport a u est d’abord effectuée entre les limites 
U one ae 2 

== et > le produit des deux intégrales précédentes pourra étre 

considéré comme une intégrale double qui est sirement absolument 


convergente. On aura ainsi 


+o 


lim || Leo ey Beene eee 


x (yx) 


U=—>a 


== feo) 


ce qui admet, comme il est bien connu, la valeur Hm(x-+v) Hn(x—»v), 
et par cela la démonstration de |’équation intégrale (62) est achevée. 
Pour v=0, les relations (61) et (62) se réduisent a (19) et (21). 
Observons que ces deux relations (61) et (62) permettent d’écrire 
léquation intégrale suivante: 


+o +00 
1 a * 4-\2. f4y)2 . (m—n) 
2 eat ixP—(v tiv? (y + jf y)m—2 [7 (2u? + 20?) du dv 
FA y) . (63) 
A 


SCE Nr eel tes 
et cette relation pourra étre simplifiée de fagon a redonner une forme 


particuliére de l'équation de WILSON (12c). 
Une équation du méme genre est 


+co +0 


1 + 4-\2, PENA es mT 
Qm pr [_[ewreroein (iu)m™—" cL (2uv) du dv=(—1)" Hinn(x) . (64) 


0) =——{2s) 


gui peut étre démontrée directement a l'aide de (24’). 
Si nous remplagons, dans (61), les variables auxiliaires u et v par 


381 


e@cos@ et @sin @, nous serons conduits, eu égard a (19), a la relation: 


+ 0 
| ele? HL (x) Hy (x) dx = 


—o 


+ 0 


== ial e (x le cos 6)? ih (x— o sin 6) H, (x—o sin A) bbe 


ou @ et @ sont arbitraires. Pour 00, on retrouve la propriété 
d’orthogonalité des polynomes d’HERMITE; pour 0=0, (65) devient une 
identité. Remarquons que le second membre de (65) est indépendant du 


paramétre 6. En particulier, pour d= nous retrouvons la relation 


aU 
Oy 
(20’) du § 3. Nous allons voir aussi l’analogue de (65) pour les 
développements en série. 

Remarquons que les formules inverses (19) et (21), écrites avec les 
fonctions de WHITTAKER correspondantes **), prennent une forme trés 
symétrique qui est peut-étre susceptible de généralisations. Par exemple, 
si m-+n est un entier pair, on aura les formules inverses 


m+n m+n+1 3 
(u) du = (—1) an 23 2a a, *Wintas m—n(20), 


u 
4 2 2 


+ 


Cue ar 2 


if 
4? 


3 
[revee) u *1Wm 1 _,1(u) Wr 
0 


m+n _ mtnrl a Lae _2 _3 
(—1) ? 2 2 [TseVae)e VW ere enced ou Wan 
0 


2 2 a 


Si m--n est impair, on aura des équations intégrales analogues, avec 
la fonction de BESSEL d’ordre + 3. 

2°. Applications. Nous avons vu’) que les polynomes de LAGUERRE 
vérifient la relation 


SEP Ley Gale ey) = Ps Li ety). . « (66) 


k=0 


Cette formule, par spécialisation des paramétres a et f, nous donnera 
des résultats intéressants. Posons, par exemple, a=f—=— 4. Si nous 
rappelons les relations (15) qui lient les polynomes d’HERMITE et de 
LAGUERRE, nous aurons la généralisation d’une formule antérieure: 


E ({,) Hane) Hanae a) = (1 md La (22+ 99 - + (67) 


k=0 


28) Modern Analysis (cité sous 14)). 


29) FELDHEIM, [8]. 


SoZ 


En appliquant aux deux membres la transformation G,.(2,—iz), les 
relations (3c) et (62) donneront: 


et cette formule est une conséquence immédiate de la formule de 
MEHLER 2°). La relation (67), si l’on y applique la transformation 
G,(2,—iz), redonne en vertu de (3c) et (68), l’équation intégrale (62). 
Nous avons supposé évidemment, pour que cette remarque ait un sens, 
que (68) a été démontré directement, sans l'emploi de (62). Elle permet 
encore d’écrire une relation analogue a (65): 


x G Fox (e) Ha n—2k (0) =| 
k=0 


(69) 
= a i Ax (0 \“ 2 cos 9) An—2x (@ LY 2 sin 0) \ 


k=0 
et l’on voit que le second membre est indépendant de 0. 
Si l'on fait ensuite, dans (66), a—=—+4, /—0O, les mémes formules (15) 
donneront lieu a la relation 
Alon 41 (x) 


——_-— ) 
Zee = 


(—1)* 22 k! ( 2) Lx (x?—y?) Honan (y). . (70) 


(y étant arbitraire). L’application de (62) a cette égalité conduit a 

+0 — 

es e— (xt iz? Alans (x ae) dx = 
DG 


ae WAS) 
= 2,2 (1) Hela + 2) iy) Hanae iy V2) 


/ 


Le second membre represente ici encore une expression qui est indépen-~- 


dante de y. Si l'on fait y=0, on retrouve une relation établie déja 
dans le=n? adie). 


30) Voir § 2, n°. 2. La relation (68) a été récemment démontrée par F. SCHOBLIK, 


Ueber eine Funktionalbeziehung Hermitescher Polynome. Monatsh. Math. Phys. 47, 
333—337 (1939). 

31) Un grand nombre des relations liant les polynomes d’Hermite et de Laguerre, 
parmi lesquelles, par exemple, les formules déduites de (31), et la relation (38), ont été 
récemment données par L, TOSCANO, Numeri di Stirling generalizzati, etc. Comm. 
Pontif. Acad. Sc. 3, 721—757 (1939), Mémoire paru aprés que le présent travail a été 
rédigé, La méthode de l'article cité est entiérement différente de la nétre, et fait usage 
d'opérateurs différentiels, appliqués a 


d” (e-x* 
a == (1 )\ee ie eect 


at(xitte-*) 


2 —x 7 (@) 
oe fh ene Ax) 
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Continuons encore a considérer des conséquences de la formule (66), 
Ou nous supposons maintenant les paramétres a et f des entiers non- 
négatifs (que nous désignerons, pour distinction, par r et s), et remplacons 


x et y respectivement par 2x? et 2y?. 
simultanément pour x et y, la transformation (21), 


Il vient 


Appliquons 4 cette relation, 


(je (0) Hes (0) Hack 0) Hrctes (2) = | 
k=0 


— 2ntrts nt | || e(xtiuyp—(y+iv) (ix)" (iy) (es (2x?+2y7) dx dy \ 


=} = Co 


(72) 


r et s étant des entiers arbitraires, non-négatifs. I] ne sera pas superflu 
d’indiquer le cas particulier de (72) correspondant au cas le plus simple 
ou r—s=O. Alors (72) devient, en vertu de (63), 


fl 
— 


k=0 


n 


k 


( 


il 


k=0 


) FH an oes Sle oa, 


(73) 


et l’on vérifie que cette relation peut aussi étre déduite de la formule 


de MEHLER. 


Nous insérons ici un résultat pouvant étre déduit aussi de (67), en y 
posant x—=u-+v, y=u-—vp, et appliquant aux deux membres les trans- 


formations (61) et (62). 
simplicité, que dans le cas 


3 (-1) 


k=] 


@ ) (2m)! Lom (o2)-+2 


Les cas particuliers correspondant 4 0=0 et 0O= 


Nous 


Ol! == Aime 


2m 


e 


) (2m—2k)! o** cos4k0.L 


(4k) 
2m—2k 


n'écrivons ce résultat, pour plus de 


(ies (74) 


= Hom (@ cos 6) H2m(o sin 8). ) 


4 


sont a confronter aux relations établies dans I (§ 2, n°. 2). 


oe 


Tt . 
— de cette relation 


Développement du produit de deux polynomes d’HERMITE de 


degrés et d’arguments différents. Passons a la généralisation de la formule 


(68) en cherchant le développement analogue pour produits de polynomes 
d’HERMITE de degrés différents. Nous partons de la fonction génératrice 


(2) qui permet d’écrire que 


3 Hn (oe + y) Ha (x2) 
m=0 n=0 M:N: 
=) 
f=01s=0) 


r+s 


as Weta tey 


rls! 


y, 2 leks (x (72) lak (y VADs 


et il faut chercher le coefficient de ¢" 7" dans la somme double du 


second membre. 
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Indiquons immédiatement la relation obtenue pour esr 


eh, (*74)= x (ke [ite ae ee eT 


qui constitue la généralisation de la formule (8). 
En général, on a 


Hin (x + y) An (x—y) = 


_ mtn m+n 


= (—1)" YD Ne ps Aves Ax (x yy?) Am+n-k (y [Y 2) 


avec 


Remarquons que At” est le coefficient de x* dans le développement 


du produit (1 + x)"(1 — x)". Si m=n, ce coefficient n’est différent de 
0 que pour k pair, et AY,” =(—1)* (Z), de sorte que (76) redonne la 
relation (68). Si l’on fait x0, on aura le développement 


m+n 
2 2k i m,n / 
Em (y) Hay) = 2-4) (1) TE ALE” Hinan ae (y) - (76) 


gui est a comparer a la formule (16). 
La relation (76) peut encore étre généralisée, par exemple, de la facon 
suivante: 


Am (x cos a+ y sin a) Hy (x sin a—y cos a) = 
iD ee n sad k (m,n) 
=(—1\" sin™ a.cos’@, 2 cotg’ a. Dy gle) ieee ay) 


k 
(m,n) __ eye tf m or 
ata an ) ( r ) eS oe 


a étant quelconque. Pour ce 7. (77) se réduit a (76). Indiquons une 


conséquence de (77). Multiplions ses deux membres par e~*’, et intégrons 
par rapport a x de —o% a+, (L’intégration terme A terme, comme 
pour toutes les séries intervenant dans ce travail, est é6videmment justifiée), 
Le second membre, par suite de l’orthogonalité, n'est différent de 0 que 
pour k—0, et alors on retrouve la relation (52). 

4°. Quelques conséquences des formules (18). Remarquons que les 
relations (18) peuvent étre considérées comme une généralisation de la 


avec 
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formule (31)*). En effet, si nous posons, dans (18), n=0, la formule 
fondamentale (13a) permet d’écrire que 


(—1)" c, (m, 0, B—m, a+ m) = (—1)’ on ) 


it 4 
3 (ayn ( é ) 
Vv=0 m—v 


et ceci est équivalent a (31). 

Nous avons vu que l’inversion de (18), sous la forme (35), est impos- 
sible, et nous avons démontré, par contre, la possibilité du développement 
(60). On peut alors se proposer de chercher l’inversion de cette formule 
(60), c’est-A-dire une relation analogue a (18): 


de sorte que 


Dee (a) Ley), 


| 


(7 (2) eee (+8) 
Lm (y) La’ (y) = = hy (m,n; a, B) La" (2y). . . ~ (78) 
n= 
Nous pouvons démontrer, par la méme méthode qui nous a conduit au 
résultat exprimé par (18) *%), que les coefficients hy de ce développement 
vérifient la relation 


hy (m,n; a, 8) =(—1)™**-* hg (m,n; B—m-+n,a+m—n) . (78a) 


(L’application de la transformation de HANKEL a (78) redonnera, en 
vertu de (12c), l’équation intégrale (22)). 

En particulier, si mn, les coefficients Ay sont nuls pour les valeurs 
impaires de k. Si nous prenons, par exemple, m—=n, a=—f=z2, et 
y =v’, la formule (78) devient 


Nn 
2 = 1 
ries (Vv) = 20 3 Peels (2r0\ nota te Berets ee 7 9) 
=0 
avec les valeurs correspondantes hy (n,n; 4,4) des coefficients. Or, si 
nous appliquons a cette relation la transformation de HANKEL de noyau 


Ji (Uuv) e-’ v4, l’équation (12c) de WiLson donnera 

{ Ii(2Vavje? 0 Hina (Vv) dv=2ut & hy LS} (2u)=e—" u Hons (Vu). 
=0 

0 


C’est une équation intégrale établie par S. C. MITRA 34) On déduit 


82) Cette formule (31) a été déja donnée, comme nous venons de le voir, dans le 
Mémoire de E. KOGBETLIANTZ, Recherches sur la sommabilité des séries d’Hermite. 
Ann. Sc. de I'Ec. Norm. Sup. 3, 49, 137—221 (1932). 

33) FELDHEIM, [8]. 

34) MITRA, Journ. London Math. Soc. 11, 252 (1936). 
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ensuite de (78), par dérivation, et en posant a=f=—+, le développement 
= = nl) 
FAlon—-1 (Vv) Aon (Vv) 2 us 2, hy Loesi(2v), . - . (79a) 


qui fournit, par application de la transformation de HANKEL de noyau 
Jo(2 [Yuv)e-’ v-, l'équation intégrale analogue a celle de MITRA, et 
trouvée par A. ERDELYI*). (Ces équations intégrales sont des cas parti- 
culiers de (22), et leur déduction précédente rentre dans la remarque 
faite au sujet de la formule (78) ). 

Remarguons, pour terminer, que (18) et (78) donnent lieu aux déve- 
loppements 


= n n k 
Han) = 3 he LE = (Hl? Fey (80) 
K=0 = Kk! 
(avec les mémes hx dans les coefficients!), et 
— bl, 
HA OS hy Le a ee Os) 
k=0 


L'application de la transformation de HANKEL a (80) (ainsi qu’a la 
relation que l’on en déduit par deérivation) ou a (80’) conduira a des 
équations intégrales démontrées également par MITRA, et MEIJER *°), et 
formant des cas particuliers de (22). 


Budapest, le 12 février, 1940. 


2) Ison, ||7I\. 
36) MITRA, Math, Zeitschr. 43, 205—211 (1938); C. S. MEIJER, Proc. Kon. Ned. 
Akad. v. Wetensch., Amsterdam, 41, 744—755 (1938). 


Mathematics. — On the thermo-hydrodynamics of perfectly perfect 


fluids. 1. By D. VAN Dantzic. (Communicated by, Prof, |e Ae 
SCHOUTEN). 


(Communicated at the meeting of February 24, 1940.) 


Summary. 


The equations of motion of a perfectly perfect (in particular of a 
relativistically perfect) fluid are brought into a general invariant form, 
independent of metrical geometry (§ 1). They are shown to be derivable 
from a simple variational principle. It states that the integral of the 
pressure over an arbitrary fourdimensional domain in space time, under 
a deformation “dragging along’ (cf. § 1) the chemical parameters 2” 
and the temperature vector 9, hence also the congruence of macrosco- 
pic wordlines, is equal to 6 x* times the virtual heat of the deformation, 
flowing through the boundary into U (§ 2). In § 3 some other variational 
relations are derived. In § 4 a result due to EISENHART and used by 
SYNGE is obtained by metrical specialisation from a one dimensional 
variational principle. 


§ 1. The equations of motion. 
The equations of motion of continuously distributed matter are 
according to EINSTEIN’) 
V7 = 0) ?) ee eg eer 


where \/j is the symbol for covariant derivation, whereas for a relati- 
vistically perfect fluid +) 
h A h 
yh O46 h 
i Oe ie Ay, )2)°) a oe  ) 
Cr Prk Fi, —4 A} Fj, Fit 4 
j= V—g( ik 444i + jk ) oF oe (4) 

!) E.g. A. EINSTEIN, [1], in particular § 17, 19, 20. The numbers between square 
brackets refer to the bibliography at the end of the paper. 

2) The suffixes h,i, j,k, 1 run independently through the range 1, 2, 3,4, corresponding 
with the space-time coordinates x’; r,s,t run through the range 5,...,4-+n, correspon- 
ding with the n chemical components /''; x,4, 4,” through the range 1, 2,3,...,4-+n, 
and A, B, C through the range 0,1,2,...,4-+ 7. 


3) The remark that (3) is valid for a perfect fluid in adiabatic motion only, was made 
by EINSTEIN [1] already; it seems to have been neglected by several later autors. 


4) In order to avoid superfluous factors V— g (gq = det gj), @ itself instead of o V4 


stands for the proper energy density. All densities except 0, 09, @ are denoted by Gothic 
letters. 


5) AG is the unit-tensor of space-time : Ae = ‘ pee The unit-tensors Ale E;, AR 
play the same role as Ay with regard to their respective ranges. 


h . . . * 
it = ce is a unit-vector along the macroscopic worldlines, i” gc elie SU 


8) 
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7 : ~h 
Usually R". is called the material part and ©S.; the electromagnetic 
part of the total tensordensity 3", of stress, momentum, and energy. 
According to MAXWELL’s equations 
z 
V7; 61 =— Fix 8 avid Wed teat Dear ae) 
where 8” is the electric current-density. Hence (1) becomes 


awl ak 
Vz CU aes | ra) . + . ‘ + + + . + (6) 


In order to bring our notations in accordance with previous papers a) 


we add to Ge the potential momentum and energy of the matter with 
respect to the electromagnetic field, which is 3" @;. In fact, let the 
substance be a mixture of n different components, in relative equilibrium 
with each other, so that they all have the same velocity. Further let 
m’ and e’ be the mass and charge of one particle (molecule, ion, 


electron, etc.) of the r* component, N° the number of particles®) of 
the rt component in a volume-element dV, MN" the corresponding 
particle-current, so that Nov =a Ni dBr%), dB, being the components 
of dV"°), and f; the potential momentum and energy of each particle 


of the r* component. Then gh— MN er, fi=e'g;. Hence the total 
potential momentum and energy in dV with respect to the field is 


No" fi — I dB, .e” Op ed Ge. Qi- 
The total stress-tensor of the matter is therefore 
nh h ( 
Viet ef =Rr te on eee 
whereas the remaining part is 
o) oh = CS ey eS) 


If the mixture is neutral, the second term in the right members of (7) 
and (8) vanishes. 


Introducing (7) into (6) we obtain 
Wey Sts Wasi oh ns Se OD 
where we used MAXWELL’s equations in the form 


Fir =2V ui pay Hae Vi pe—V ie pis Vi FV —9 = 8", whence V3/=0. 


”) D. vAN DanrTZzie, [2], [3], [4], referred to as Ph. Th., R. Th., and R.G. respectively. 


8) If we wish to avoid any molecular model, we may read ney as the number of 


gm. mol, of the rth component, contained in dV, and gr as the corresponding current- 
: r 

density. Then 2", ", x?, f7, m™, e” also h®Ve to be reckoned pro gm. mol. instead of 

pro molecule. N.B. Summation convention, for r, s also. 


®) By =ar and atr= we denote an equality defining its left and its right member 
respectively. 


10) Cf. Ph. Th. p. 684. 
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Equation (9) remains valid if other than electromagnetic forces work upon the 
matter, provided they are derivable from a vector potential, which is the 
same for all particles of each component of the mixture, and if we take 
for f7 the sum of e” y; and this vector potential. 

According to (3), (7) the mixture is what I called a “perfectly perfect 
fluid’ ''), the stress-tensor and particle-densities having the form !?) 


Oe A ee, eee (10) 

i oth alee gy aes weet 1) 

where J! = i? by) =ci"/k Ty") is the “temperature-vector”, T, being the 
proper temperature and k BOLTZMANNs constant, whereas®) 

eC MOct ep) Wie? fi. o aoe 2 os ea (12) 

According to Ph. Th. p. 697, R. Th. p. 604, generalised for a 


mixture of n components. p can for constant x‘ (i.e. in each given point 


of space-time) be considered as a function of n-+4 independent variables 
Oo and a’ with 


ee ee 
Ps = a5" row 


(13) 


If the dependance upon the coordinates is taken into account also, p depends 
upon them: 1°. through the 0” and 2’, 2°. through the f’, which enter into p 


in the combination A” + fj 9 only, 3°. through the gi;. If the curvature is 


negligible and no direction in empty space is given, p will have the 
form '*) 


Deda tof | 7 ==, Oy =ae 0 in = VV giy 8 OF =clkTp. (14) 


We suppose @ to be a function of its n+1 arguments alone (which 


nN, [2in, Ma, Os Noe I, Toy, jo, COX. 

12) A very thorough discussion of the relativistic stress-tensor on an axiomatic base was 
given by J. L. SyNGg, [1]. I find however some difficulty in accepting SYNGE’s hypothesis 
of “elementary impulses’. For either the “gas’’ consisting of these elementary impulses 
would behave in a way entirely different from any other gas with respect to its thermo- 
dynamic properties, or otherwise it could, like a photon-gas, be considered like one of 
the components of our mixture, admitted the alterations mentioned in R.G. § 10. In that 
case however it would be difficult to see, how the elementary impulses could have any 
appreciable influence, except at very high temperatures. In any case the complete inde- 
pendence of the interaction between molecules and the temperature would be hard to 
account for. 

13) In order to avoid factors c in the general (non-metrical) part of the theory, where 
they have nothing to do, we may take with respect to orthogonal coordinates x4 = ¢ 
(instead of xt=ct), g44=c?, gtt=c—%, etc. This leads to some unusual factors c, but 
makes comparison with the non-metrical as well as the classical theory easier. This system 
of notations was introduced by DE DONDER [1]. 

14) Cf, e.g. for an ideal gas R.G. (51). 
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themselves are functions of x”) and not to depend upon the x’ explicitly. 
Then, putting 9’ ae 97/00. Pr =ar Og/dd', 


pi Or Y—g, 
p= orf +e i)\V—s, 
EV See ee ey) 


Denoting by 0; the operator 0/dx! under constant ot and 4‘, and by 
0; the complete operator 0/0x', so that 


0 
0; =ae 0; — (07 7”) zs — (0; 2’) ples 


dor a 


a short calculation shows that (9) is equivalent with 


In this equation the metrical quantities do not occur anymore. Though 
its terms themselves are not invariant under arbitrary transformations 
of coordinates in space-time, the equation (17) is invariant as a whole, 
38"; having the form (10). At the other hand it can be shown that it is 
also necessary that SS"; has the form (10), ie. that the fluid is perfectly 
perfect, in order that (17) be invariant, as long as the operator 0; is 
defined by (16). We may therefore consider (17) as the equations of 
motion of an arbitrary perfectly perfect fluid, independent of metrical 
geometry. 

From here on, except when the contrary is stated, all equations are 
independent of any special assumption concerning the relation between 
energy and momentum. Hence they are independent of the axioms of 
relativity theory and also of metrical geometry. 

Inserting (10) and (16) into (17) we obtain 


830; pi: +7; 0:07 +7;0; 7 =—p-do;4".. . . . (18) 


Both sides of this equation are invariant quantities. 
The differential equations 


Fi sone) tt Mae ot memes (1.9) 


define a one-parametric group over a part of space-time. Its infinitesimal 


oe) 


transformation is determined by dx'—=y d@i", d6 being an infini- 
d ; 
tesimal increasement of the parameter @. Its LiE-symbol is 46 A DOr 
The infinitesimal transformation transforms any scalar function f into 
f+ df/d6. The corresponding invariant operator working upon a general 
geometric object X: (where the big points stand for rows of indices) 


was first defined by SLEBODZINSKI ). It will be called the Liz-derivative 
d 
and denoted by a It is defined as follows. The components of 
X: + d;,X* in a point (with respect to the system of coordinates under 
consideration) are the components of the value, Xi takes in a point 
x"+dxt—=xt+t o*d6 with respect to the system of coordinates, 
obtained by “dragging along” the original coordinates along dx", i.e. by 
ascribing to each point y as new coordinates y"’ the original coordi- 


nates of the point from which it was obtained by the infinitesimal 


d 
transformation: y’” =y"— #"(y)d6, By this definition — becomes an 


dé 
d 
invariant operator. Applied to a scalar f we get aia at Applied to 


a scalar density (or W-density '°)) p of weight + 1, to a contravariant 
vector v" and to a covariant vector w; we obtain !’) 


HP a, po = Ze + pa,os. ae ee 20) ) 

SLU = 913, ot —of 0, =F of, 9%, . (21) 9 

Fi — 2.9) Oy wy +1 wy = Ot 0,01. a 
Evidently 

Shee al WE Beton 60 sess, eke RIS 


1b) W. SLEBODZINSKI [1]. Cf. also J. A. SCHOUTEN and E. R. VAN KAMPEN [1], 
D. VAN DANTZIG [5]; J. A. SCHOUTEN and D: J. StTRUIK [1], p. 142. 

16) W-densities are quantities which under transformations of coordinates are multiplied 
by a power of the absolute value of the transformation-modulus. They were introduced 
by H. WEYL, e.g. [1] p. 98. Cf. also J. A. SCHOUTEN [1], J. A. SCHOUTEN and D. 
VAN DANTZIG [1]. 

17) Cf. J. A. SCHOUTEN and E. R. VAN KAMPEN, loc. cit., p. 4. 


oS etc.) works upon all following quantities 
0 

untill either a closing bracket belonging to an opening bracket preceding the operator, or 
a +, — or = sign, or the end of the formula is met. The symbol for a differential, a 
variation or an element (e.g. d, 0, d) belongs to the immediately following symbol alone. 


18) A differential operator (e.g. V jr dj, 


& ‘ ea) 
19) Introducing the Lie-symbols X4 =dt Go =at O/ Oj and Xz =a v/ 07, we have 


(2 v/) 0; = (Xj, Xo), where the right member is the ordinary LIE-bracket. 
do Ae 
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Applied to a covariant vectordensity (or vector- W-density) pi of weight 
+1 we obtain 


= Nd;p, +p; 01 +9: 0; %, « Aer Oat (24) 


dy Pi 
dd 


which is the left side of (18). Hence we have proved: 
Theorem 1. The equations of motion (17) are equivalent with the 
system of equations 


§ 2. Derivation from a variational principle. 

In this paragraph we show that the equations of motion (17) torca 
perfectly perfect fluid can be derived from a simple variational principle. 

Let U be a part of space-time, bounded by a differentiable threedi- 
mensional manifold B. The latter is provided with an exterior orientation *°), 
viz the direction outward from U. An element of U will be denoted by 
dU, its fourdimensional volume (as measured by the coordinates under 
consideration) by dll; an element of B by dB or dV; its components by 
dB; or d¥;. Let further a contravariant vectorfield z*—=z"(x') be 
defined in U, the components of which have continuous first derivatives 
with respect to the x!. 

In order to determine easily the variations of several quantities it is 
useful to introduce the Ligz-variation 6; with respect to the equations 
dx"|du=z" (x) in the same way as dp=— do ct was defined with 


dé 
respect to (19). 


The variation of an integral Wau { Wal under the infinitesimal 
U 


transformation 6 =a 6x/ 0; =ap Ou.z/ 0; is then 
WH dW | (6,W)dU= | 0) Wow) d= [ W bai ABy, (26) 
U U B 


as by definition 6,dU—0, the new element dU being obtained by 
“dragging along’ the original element. It is of course trivial that 6W 


reduces to a boundary integral as the variation consists of a displacement 
alone. 


20 : : P : ; : : 
) An exterior orientation of a simplex *' is an orientation for each simplex of com- 


plementary dimension having exactly one point in common with »; & itself need not be 


oriented. Cf. O. VEBLEN and J. H. C. WHITEHEAD [1], p. 56, and J. A. SCHOUTEN 
and D. VAN DANTZIG [1]. 


38) 


Now we take in particular YW =a. p =p (x", 0", 17), In U a congruence 
of curves is defined, viz the macroscopic wordlines of the motion, 
satisfying the equations (17). We now define the operator 6r, working 
upon a quantity X+ =a X%(x", 0,1"), depending upon 0 and i", by 
requiring that not only the coordinates, but also the values of 9” and J" 
shall be “dragged along’. More precisely: X:-+ 61 X% shall be the 
components with respect to the coordinates mentioned above of 


lee NTA “where "x" =x 1 dx" whereas “ and ‘1? are 
defined by requiring 6, 0"—0, 6° 47 =0. Hence 

One dee, epee OLN ee Ox: 

Of Xs S=a¢ ee + . + (27) 


In particular 
= CL Or OL — 7 Op Al 
= dp + pd; dx/—p;6 9 +p, 0/0; dx!—p, d=) . (28) 
== 0° p 4 5B, Oy dad, 
where 5° =a 5x! 0), 0; being defined by (16). 
Hence 
at | pat= | (Ox! dj p+ 2; 0; dx) dU= | 
U e 
: (29) 
= [wt ox'a%s + [ox @lr—a, waar | 
B U 


/ 


Hence we have proved: 
Theorem 2. The system of equations of motion (17) is equivalent 
with the variational principle 


tf rat={ Prox! See a» B0)' 24) 
U B 


for any part U of space-time and any continuously differentiable varia- 
tion 6x', and also with the variational principle 


at fp dWe=0 
U 


for any part U of space-time and any continuously differentiable **) 
variation 0x!, vanishing at the boundary B of U. 


21) pe =ard B, pt is the amount of energy and momentum contained in or flowing 


through the element dB. Analogously ney So”, etc. Ch Ph.- Th. p. 684. 


22) These conditions will not be mentioned explicitly further on. 
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This theorem seems to be sufficiently remarkable, to restate it in more 
explicite and intuitive terms. Let Iu [p] be the integral of the pressure 
» over a fourdimensional domain U in space time. A small deforma- 
tion x? > x =ap x"? +6x" deforms U into a new domain ‘U. Moreover 
the congruence of macroscopic worldlines in U is deformed into another 
congruence in ‘U. If the parameter 0, defined by 8% (x)=dx"/d6, is 
kept constant, the vectorfield 0”(x) passes into a new vectorfield 
“9h (x0) ap OF (x) + d(dx")/d6, which differs from GF ('sc) = OR (x) + 60" 
by dy, 0x", Leaving the values of 4” unaltered, i.e. taking CRA x ara X) a 
—j"(’x)—6,, the variation of the pressure integral is the difference of 
the integral of /p—ar p(’"('x), (47 (‘x), ’x") over ‘U and the original 
integral over U. This variation 6°J= Jy ['p]—Ju [p] is according to 
Theorem 2 (if the coordinates are such that x4 may be interpreted as 
the time) equal to dx‘ times the virtual heat belonging to the defor- 
mation, flowing through te boundary into U, if and only if the equations 
of motion are satisfied. 


§ 3. Other variational relations. 


The variational principle (30) can be brought into another form. 
Therefore let Sy be a part of a threedimensional manifold in space-time, 
which has at most one point in common with any (macroscopic) world- 
line, and which nowhere is tangent to such a worldline **), and let Oo 
be its twodimensional boundary. The threedimensional manifold Sy will 
be defined as the locus of all points x which have a difference of 
“thermasy” with Sp) equal to G6, i.e. such that the worldline through x 
intersects Sy in a point x) and that the integral 


[aga [kT a= [ kToas bess eH) 


taken along the worldline through x is equal to 6. Let U be the locus 
of all S, for 0O=G6=8,. Its boundary consists of So, Ss, and the locus 
B’ of the parts of all worldlines through Op) between Sp and S,,. 

We suppose U to be divided into elements dU of the same nature 
as U itself, ie. bounded by parts of two hypersurfaces Sj; and Sy + as 
and by oo? parts of worldlines passing through the boundary of an 
element dV of S,;. Then 


dU=dxt'dB,=0'dBrdd . . . . . . (33) 


6, 6, 
Jrat= [a0 [arorp= faa { 2, mee er 
U 0 Sp 


0 Sp 


and 
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or 
eS ‘pdll (35) 
dé rent) en 
Sq Ug 


if Uy is the part of space-time bounded by B’ and S,, Sy. 
Now let a deformation 6x" of Sy in itself be given, so that 


One ire ON te ie tet, (36) 
for any element dB; of Sy and let 5x” vanish at the boundary Oy of Sp. 
Let every element dB; of Sy as well as the displacement dx” be 
»dragged along” by the transformation dz, so that 

OO AC k=) eee 37 ee} 


and d,d¥;=0. Then (36) remains valid for every element of every Sy. 
Moreover 6x"=0 on B’. 
By (29) and (37) we have then OL p=d.p—D,; 6,4", hence 


Und oe to OUN con eee seer (38) 


hence, using (37) again, 


 frau= s frau ua [rst av pee (3c) 
0 Sg 


By (26) and the boundary conditions however 


Dy 
9 fpdu=[poxadi=| [pox ae] 06 oe (GU) 
U B 5 : 


6 


The right member of (30) need also be extended over Sp and Sy, only, 
so that we have proved: 

Theorem 3. The equations of motion (17) are equivalent with the 
variational principle 


Oy 
5, 
—faof ni’or=| [ri ox| ae ee CAN) 
0 
0 Sh Sp 


for every variation 6x! of each Sy in itself, vanishing at the boundary 
of each Sj, and dragged along by d_. 


23) According to 19) this condition means simply that the deformation $1, and the trans- 


formation dr, are interchangeable. 
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A variation of this type leaves the congruence of macroscopical world- 
lines invariant. If we assume each displaced tube of worldlines to be 
filled up with the same amount of matter of each kind as the original 


one, we have 5N?’ =0. Then, putting 
Ne oh et ea nem 12) 


(41) becomes equivalent with 


(43) 


This variational principle shows a great analogy with that of ordinary 
point mechanics, viz 


ty 


§ | Lae=[pidx'l, ee a ee 


to 


From (33) and (28) we obtain by Ph. Th. (58), (60) for an arbitrary 


variation 0x": 


Oy 
vf pau= | (au {ao { aBro" a 
U ‘y 0 


Sp 
9; 
=| do | dB B” (dn Pp—Pr Or A" — 9; On ) = 
0) 


Sp 


Oy 
=| ao { (0. 2°” —N3" dx A’ — PP" 6, H— Bt br, dB), 
(0) 


Sy 


hence 


6, 
ef pau={ao{ (oy, SAY + Ar og NEY +H by Pi” —pd* 6d By). (45) 
U 


0 Sq 


According to (30) this must vanish for every variation dx" vanishing 
at the boundary of U. The integrand of (45) is the quantity which 
would have to vanish if the operator 5, were replaced by a variation 6 
of the variables 0" and 4 alone, without change of the point xp 


(Cf. Ph. Th. (57)). 
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Hence we have proved: 
Theorem 4. The equations of motion (17) are satisfied if and only if 


| 9, 
| dé { (o S2v + ar by, NEY 0" or, P= 6 Coral by e=(0) 


0) Sq 


(46) 


for every variation }x" vanishing at the boundary of U. 

Evidently the last term in the integrand of (46) only defines 6, dj. 
Deforming dB; simply by “dragging along’, we get 6,.d%;—0, and 
this term vanishes. 

The right member of (30) becomes 


Oy 
= ~~ 9, pe ‘ 
[rioxs| [eitox |) + fae [acu oatax, an 
B Sy : e % 


where dOn; are the components of an element of the boundary Oj of S; 
with appropriate exterior orientation. Hence, if dx” vanishes at the 
boundary of each Sy (though not necessarily over Sg and Sy, themselves) 
and if 6,d¥;—0, we obtain, equating the right members of (45) and (47): 
Os 
d 


J aa | (res) Sey aI ele da PHY 8 x!)\=0. 
0 Sq 


Differentiating with respect to 0, and replacing 6, by @ we get 


[os £4 one +0 oPiY—axt & pi) =0, seen (48) 
Sy 


where the invariant LIE-variation 6, has been replaced by the uninvariant 


variation 0 again. 


§ 4. Variational principles of Fermat's type. 
Let us suppose now that the fluid moves in such a way, that a function 
y = 7 (x") exists, such that 


D0; =— 4 07 ©, Can nus ro eee em 40) 
Then the equations of motion (25) become: 


dy 
Sa ; — iW. . + . + . + . . 50 


Transvection with # shows that the density 


ie em eI Tire ne er ye (51) 
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is a constant of the motion: 


da 
—————— oN oe ae Saiyan oe 
9 = 9 (52) 


Moreover variation of the integral 


6, 6 
re dé =fo' p; Cie! a ee ee (53) 


A oh 


with 56—0 at the boundaries, gives by (50), (51): 


6, oh 
= [erdpda= [ dx! dope. as) 


6 50 


Hence, if the variation is performed by “dragging along” 0” (which is 
equivalent with the condition that dx" is dragged along by the operator 
d,), ie. 60% = — dy dx"/dO=0, then the integration in the last member 
of (54) can be carried out by using (52) and we obtain: 
Theorem 573), If (49) is satisfied throughout the fluid, then 
C=ae7pi ev 


is a constant of the motion and 


64 6, 


fe d= 0 [ o- Dida [on pox 


uh Ao 


for any variation 6x" with 6,0" = 074), 
Corollary: If the fluid is homogeneous and if 


Op == Osteen 2) oP eer ee ee 5) 
throughout the motion, then 
5, 
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Proof: By hypothesis n=1. We drop the suffix r and write p’ for 
Op/d4. Then (56) is equivalent with 


A= — 0, 2; HW = — 2; W 0; log (— 2; 94), 


232) The theorem remains valid if to the right member of (49) a term of the form 
—7i is added and if the first statement is replaced by “a a is a constant of the 
motion’, a being defined by (51). If to the right members of (50) and (52) a term y Pi 
and yQ respectively is added, the demonstration remains valid. 


24) In that case the values of @ can be kept constant during the deformation, so that 
6@ vanishes identically, and that (54) is valid. 
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where 2; =a. pi/p’, or 


pr 0; Ar =i 0; A=— q 0; log (—a; dO) 
so that (49) is satisfied with e? = — 2,0. Substitution into (52) leads 
immediately to the proof of the corollary. Evidently (52) holds with 
a——yp’, ie. the condition (56) implies that the equation of continuity 


(cf. § 5 (68)) is satisfied. The equation of continuity, however, is not 
sufficient in order that (56) be satisfied. 

As a particular case let us introduce metric again and assume that an 
equation of state of the form (14) exists with n—=1 (homogeneous 
relativistic fluid) and f; 0. Then (49) becomes 


On A180 On Me een) ote ae oa es (98) 
with 
Bo = Of in = V gi; 0§ NH = c(k Ty) 
and 
s=|7 gu ai ay, Nj—— €1j. 
Hence we have also 


Vip=p Vid—p' ¢ Vid, =D" € 9 0; (p—log %). 


A relation u0;v—0;w, however, implies that w and u are functions 
of v alone and that u=dw/dv. 
As p’, ¢ and % are always 0 (even > 0), » — log % and then also 


p’ 9) =—9;0' =—q=e-+p (cf. () and (10)), hence also @ is a 
function of » alone?) and we have 
dp 
p= log I+ | ob = logy fo sO ee (59) 
Or also, introducing the “index-function” *°) F =a: e?/J), whence 
as ee et (OU 
fog 18 “Al at (60) 


the integral (53) becomes by (49): 


6, 6, 9, 
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25) This means that the quantities q etc. depend upon the coordinates x? through p 


ae dp for all variations 6 of the form 56=6 xJ 0j, but of course not 
g 


for arbitrary variations. Relations of this type with respect to the spatial coordinates alone 
are called by BJERKNES barotropic relations. Cf. V. BJERKNES Eel) jos Cex 
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alone, ie. that dq = 
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Moreover a! becomes 


dF 


— 62 
Pra ( ) 


et Gy i ae 


whence 


a! pj = (a9) py OF 7 = Fi. 


Hence we have proved, the argument being reversible: 

Theorem 6. If a homogeneous fluid has an equation of state of the 
form (14) and moves in such a way, that @ is a function of p alone, 
kT) (dF/dp) is a constant of the motion and 


9 [erdo=9 | Fds=[Fi 0x1 ee 


(where F is given by (60), whereas F;=a Fi;) for any variation 
satisfying the conditions of Theorem 5. 

The variational principle (63), was found in 1924 already by L. P. 
EISENHART?’) by a method based on the geometry of paths. It was 
proved in a different way and discussed further by J. L. SYNGE [2] p. 
393. We obtained it here as a metrical specialisation of (54). 

Finally we apply theorem 5 to a volume, filled with black radiation, 
which can be considered as a homogeneous fluid with 4—0O and 
q=p! =— 4p (cf. RG. § 10). Then (49) is automatically satisfied 
with »y=0, so that J simply becomes 


1={ do=| KT dt = | kT y ds 


It is of importance to note that the proper mass of a volume filled 
with radiation (being greater than the sum of the proper masses of the 
photons, which is zero) is always positive, viz 3p d¥/c?. Therefore the 
average velocity of black radiation is always <c (e.g.=O0 if the 
radiation is enclosed in a box at rest). Hence ds +0, so that i*, To, 
etc. remain finite. Of course this does not apply to a single beam of 
light, which is a limiting case with p—p;0!'=0. Moreover (14) is 
satisfied with 2=0. Hence F=%'=kT>. Finally we remark that 
by theorem 5 a=q, hence yp, hence Ty is a constant of the motion, so 
that we have proved 

Theorem 7: Black radiation moves always according to the variational 
principle 


9, 


9 | da=—4 [pp dx Sean aaverE 


Do 


77) Cf, L, P. EISENHART [1] p. 214—216. 
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which becomes by mettical specialisation, T) being constant along the 
rays, 


| ds= [ir dx4) ee ee 


Hence the “stream-lines’ of black radiation are geodesics. 
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Anatomy. — The cytology of the cortex in the opossum (Didelphys 
virginiana) considered in relation to some general problems of 
cortical evolution. By W. RIESE and G. E. SMytH 1), (Communicated 
by Prof. C. U. ARIENS KAPPERs.) 


(Communicated at the meeting of February 24, 1940.) 


In the extensive literature devoted to the comparative anatomy of the 
cerebral cortex it would appear that in the cytotectonic descriptions too 
little attention is payed to the differentiation of the constituent cells them- 
selves. So the cortical lamination of the opossum has been carefully studied 
by Gray (1924); the excitable areas of the cortex of this animal by GRAY 
and TURNER (1924) and the origin and course of the pyramidal tracts by 
TURNER (1924). Nevertheless little or no information is available on the 
structural differentiation of the cells in the several cortical laminae and in 
representative cortical areas. 

In the present paper we shall deal with some cellular differentiations in 
the cortex of the opossum as an addition to our recent observations on the 
cellular structure of the thalamus and corpus striatum and other subcortical 
centers of this primitive mammal 2). 

In consequence of the lack of systematic phylogenetic studies of cellular 
differentiation in mammals no absolute standard exists whereby the degree 
of differentiation of nerve cells may be tested. In these observations the 
criteria adopted are based on ontogenetic studies, in particular on the 
conclusions arrived at by RIESE (1939) from investigations on the 
development of the brain of the bear (Ursus arctos). They are in close 
agreement with some earlier observations on the development of cortical 
cells in pig and rabbit embryos reported by PATON (1900). 

Realizing the fundamental difference between embryonic and adult cell 
structures, also in lower mammals, the comparisons established in this 
paper are only meant to indicate certain analogies. With this reservation 
the criteria of structural differentiation may be summarized as follows: 

1. The amount of cytoplasm, particularly in relation to the bulk of the 
nucleus, and the sharpness of demarcation between cytoplasm and the 
surrounding substances. The more primitive the cell, the more poorly 
demarcated seems the cell outline. The bulk of the cytoplasm relative to 
that of the nucleus tends to increase with the progressive evolution of 
the cell. 


1) From the Laboratoire de Physiologie générale de la Sorbonne et Laboratoire 
d'Ethologie des Animaux Sauvages au Muséum National d'Histoire Naturelle. 
2) These observations will be communicated soon. 
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2. The intensity with which the cytoplasm is stained by basic solutions 
both absolutely and relatively to the nucleus. This is a measure of the 
amount of chromophil substance present in the cell since the diminution 
of the chromophil substance in the nucleus and its increase in the cytoplasm 
run parallel with the maturation of the cell. 

3. The presence of distinct chromophil granules in the cytoplasm. It 
seems that their presence indicates a high degree of differentiation. 

4. The form of the nucleus and, more particularly, the number of 
nucleoli and their position within the nucleus. The presence of numerous 
nucleoli, more especially when disposed peripherally is caracteristic of a 
primitive type of cell. 

5. Finally the position of the nucleus in the cell, central or excentric, 
and the degree of development of cell processes appear to be of some 
significance in conjunction with the criteria already enumerated. 


Material and methods. The opossum (Didelphys virginiana) on which 
the investigations were made was a young adult male of good nutrition 
which had died suddenly from an acute cardiac affection. 

The brain, after fixation in 10% formol, was embedded in celloidin and 
sectioned serially in the vertical plane. Alternate sections were stained by 
the methods of Nisst (Thionin blue) and SPIELMEYER. 

There is little difficulty in identifying the principal cortical areas with 
those described by GRAY and TuRNER (l.c.). The degree of structural 
differentiation of the nerve cells of the various laminae was examined in 
the following areas: 


1. Area praepiriformis — Palaeo-cortex (KAPPERS) 
2. Cornu ammonis — Archi-cortex (KAPPERS) 
3. Area praeorbitalis .— Neo-cortex 
4. Area parietalis — Neo-cortex 
5. Area striata — Neo-cortex 


No attempt has been made to study the cellular structure in each cortical 
area. The study of some clearly defined and representative areas would 
best serve the purpose of this investigation, which is intended only to find 
some principles indicative of the general line along which the cortex develops. 


1. Area praepiriformis. This cortex is composed of three laminae of 
which the 1st and 3rd are deep and the 2nd is narrow. The 1st lamina is 
very poor in cells and such as occur resemble in every way those of the 
second lamina. 

In the 2nd lamina the cells are extremely numerous, medium sized, 
polygonal and seem to be scantly provided with processes. The cytoplasm 
stains deeply and diffusely. The nucleus stains less deeply than the 
cytoplasm and contains, as a rule, one large and prominent nucleolus. 

In the 3rd lamina the cells are less dense but still very numerous. They 
are somewhat larger than those in L.2 and their size increases slightly 
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in those lying nearer the white matter. Many pyramidal forms are seen, 
but the majority is polygonal. 

Each cell has several prominent processes. The cytoplasm is distinctly 
paler than in L2, and contains numerous small chromophil granules. 
The nucleus is large and round, usually excentric in position and 
provided with a distinct nuclear membrane. It stains lightly and possesses 
a single darkly staining nucleolus and a few threads of chromophil 
substance. The largest cells located near the white matter, often show a 
vacuolated cytoplasm. A considerable glial satellitosis is frequently seen. 

Apparently in this most primitive cortical area the individual cells are 
very mature. It is significant that the cells in L3 are somewhat more 
highly differentiated than in L2 and that, in fact, the degree of 
differentiation shows a progressive increase towards the white matter. 


2. Cornu ammonis. (Archicortex (KAPPERS), Cortex totoparietinus 
bistratificatus holoprotoptychos (ROSE), cortex heterogeneticus rudimen- 
tarius (BRODMANN).) 

The cells of the ammonal area vary in size from that of medium to rather 
large cells. The majority are pyramidal but many ovoid and many pear 
shaped forms occur. The cytoplasm is scanty in relation to the size of the 
nucleus and contains numerous coarse fragments of chromophil substance. 
The nucleus is round, pale and possesses, as a rule, one deeply staining 
nucleolus. The cell processes are well developed: many are directed into 
the lamina zonalis in a vertical and parallel sense. The description agrees 
with that given by ROSE (1927) for Didelphys azarae. The cells of the 
fascia dentata, generally described as “granular” cells (ROSE, l.c.), only 
show a narrow but clearly visible ring of cytoplasm showing a light pink 
stain. The nucleus is round and contains numerous granules of dark 
chromophil substance located peripherally. These cells resemble those of 
the habenular ganglion of this animal. Apparently the cells of the cornu 
ammonis are highly differentiated, showing many of the characters of 
mature ganglion cells, while those of the fascia dentata are of very 
primitive type. 


Area praeorbitalis. The lamination in this area is very distinct. The 
4th lamina (1. granularis interna) is lacking i.e. the cortex is agranular. 
The Ist lamina (1. marginalis) is broad but only contains a few scattered 
cells towards its deeper part. In the 2nd lamina (1. granularis externa) 
the cells are very densely packed. They are mostly of triangular shape 
with the apex directed towards the surface. These cells stain so deeply 
that their internal structure cannot be defined. Each cell possesses several 
long and slender processes. The 3rd lamina (1. pyramidalis) is also broad 
and contains many cells. These cells are larger than in L 2, polygonal and 
possess rather short processes. The cytoplasm is abundant and stains 
lightly. No Niss_ granules are seen. The nuclei are large, round and stain 
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more deeply than the cytoplasm. Each of them contains a net-work of 
chromophil substance and one or two dark nucleoli. In the 5th lamina 
(1. ganglionaris) the cells are generally pyramid shaped, their orientation 
being such that the apex is directed toward the white matter. It is to be 
noted that by their medium size these cells are very different from the 
giant cells of the motor cortex of higher mammals. The outline of the cell 
is well delimited and the processes are long. The cytoplasm is rich in 
tigroid granules only their axon hillock is devoid of these granules. The 
nucleus, located centrally, is pale and contains one deeply stained and 
centrally placed nucleolus. These cells are relatively scanty, in number. 
The 6th lamina (1. multiformis) is poorely delimited towards the white 
matter. The cells are rather scanty. They stain deeply but nevertheless 
their boundaries are rather ill defined. They are polygonal, somewhat 
smaller than those of the 3rd and 5th laminae and their processes are only 
seen with difficulty. Their nuclei, as far as may be observed, only contain 
one dark nucleolus. 

The cytotectonic arrangement in this area as well as the differentiation 
of the cells themselves evidently are highly evolved. As in the areas 
previously described their evolution progresses from without inwards. 


Area parietalis. As in the other areas only an occasional cell is found 
in the deeper part of the ist lamina. In the 2nd lamina the cells are 
closely packed to form rather narrow zones. They are polymorphic, poorly 
defined and each possesses one or more short processes. The cytoplasm 
stains so deeply and uniformly that only occasionally is it possible to see 
the still more deeply stained nuclei with their numerous nucleoli. In the 
3rd lamina the cells are pyramidal in form and sharply outlined. The 
cytoplasm stains with moderate intensity and Nissi granules are absent. 
The nucleus is relatively large, excentric in position. It usually contains 
many nucleoli, but occasionally one nucleolus stands out more prominently 
than the others. From the top of the cells a big process arises directed 
towards the lamina zonalis. The 4th lamina is difficult to define, its cells 
resembling closely those in the 3rd lamina. Generally, they are round or 
ovoid in outline. The cytoplasm is much less deeply stained than in the 
cells in the external granular layer but there is no trace of granules. The 
nucleus accounts for the chief bulk of the cell; it stains more deeply than 
the cytoplasm and contains numerous chromophil granules which show a 
tendency to be grouped at the periphery. Short processes are seen, but 
only occasionally. The 5th lamina is a distinct layer of medium sized nerve 
cells showing a linear arrangement of several layers. The cells are sharply 
defined and have one thick apical process and finer lateral and basal 
processes. The cytoplasm is deeply but irregularly stained as if from the 
presence of coarse masses of chromophil substance. The nucleus stains 
less deeply than the cytoplasm and contains a solitary dark nucleolus. The 
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6th lamina merges with the 5th but its cells are more scanty and their 
orientation is more irregular. They vary in size and, in general, are very 
deeply staining. In those cells in which the nucleus is seen, it is found to 
contain one nucleolus. A moderate glial satellitosis is of common occurrence 
in this layer. 

In this area cytoarchitectonic differentiation is in advance of the 
differentiation of cellular structure. In the 5th lamina the nerve cells are 
less highly evolved than in the corresponding lamina of the area prae- 
orbitalis. It is noteworthy that a considerable difference exists between 
the structure of the cells in the two granular layers. 


Area striata. The Ist lamina, as in the other areas examined, only 
contains scattered nerve cells. Such as occur are small and so intensely 
stained that their internal differentiation cannot be distinguished. The cells 
in the 2nd lamina form a narrow but closely packed zone and their staining 
property is identical with those in the Ist. They are provided with short 
processes, which gives them a spiderlike appearance. The 3rd lamina is 
composed of medium sized pyramidal cells arranged in parallel rows. The 
cytoplasm is deeply and uniformly stained but the nucleus is nearly always 
visible. Their nuclei are to intensely stained to ascertain the internal 
structure. Each cell has two or more rather slender processes. The 4th 
lamina is deep and very distinct from the adjacent laminae. It is composed 
of medium sized fusiform cells, a trifle larger than those of the 2nd lamina. 
In some instances slender spider-like processes are seen. The cytoplasm 
stains far less deeply than that in the 2nd layer and with a slight pink 
colour only. Distinct NissL granules cannot be observed in it. The nucleus 
is relatively larger, generally excentrically placed, and more deeply stained 
than the cytoplasm. Each contains numerous chromophil granules of which 
two or more stand out prominently. In the 4th layer of the area striata of 
the opossum there is no trace of sublamination such as is found in higher 
mammals. Although in low magnification a distinction can be made between 
the 5th and 6th laminae, it is impossible to do so on the basis of a 
difference of cellular structure. The cells in these layers are more widely 
dispersed than in L4 and their size tends to be somewhat larger. The 
cytoplasm stains less uniformly but in all other respects these cells are 
hardly to be differentiated from those of the 4th lamina. A certain richness 
in glial elements exists in the 6th lamina and a marked satellitosis is not 
uncommon, This phenomenon is most evident and most frequent around 
displaced ganglion cells which are to be seen almost as deep as the 
ependyma of the lateral ventricles. 

The most striking feature of this area is the primitive structure of its 
cells, strongly contrasting with its well differentiated cytoarchitecture 
(orientation and degree of density of the cells). Evidently the cells in the 
deeper laminae only are more highly organized. 
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Discussion. 


From our observations it appears that, although an equal degree of 
laminar differentiation is not present throughout the whole cortex, this 
laminar differentiation is in advance of the degree of cellular differentiation 
in the areas examined. In some areas the disparity is very great. In all the 
areas examined the most highly differentiated cells are constantly found in 
the 5th layer. 

These observations agree with the fact observed in the ontogenetic 
development of the cortex in which the formation of the cytoarchitectural 
pattern always preceeds the differentiation of the constitutive cells, and 
with the fact that at a period when the cells composing the lamina 
granularis primaria still are in the neuroblastic stage, certain deeper cells, 
separated from it by a zone poor in cells, already possess the characteristic 
form of ganglion cells. In this layer, the primitive 5th cortical lamina, 
cellular differentiation begins. Nothwithstanding these analogies between 
the mature state of the cortex of a primitive mammal and certain embryonic 
stages in higher mammals these analogies do not, of course, mean that in 
both cases the same causal factors are in operation. 

A comparison between the degree of cellular differentiation in the 
cortical areas examined shows an additional interest viz. that a particular 
degree of cellular differentiation is characteristic of the palaeocortex, 
archicortex and neocortex. Thus, the cells of the area praepiriformis 
(palaeocortex) and those of Ammon’s horn (archicortex) are very highly 
differentiated, whereas the cells of the neocortex, with the exception of the 
area praeorbitalis, are in general poorly differentiated. 

It follows that the subdivision of the cortex into palaeo- and archicortex 
on one hand and neocortex on the other corresponds with the degree of 
cellular evolution: indeed the constitutive cells of the recent areas are less 
highly organized than those of the ancient and more primitive areas. 

Before concluding this discussion two observations of minor interest 
seem worthy of brief mention. 

If the internal structure of the cells composing the second and fourth 
laminae are compared, it will be found that they show considerable points 
of difference. In all the areas examined it is true to say that the cells of 
the internal granular lamina have a slightly more advanced structure than 
those of the external granular lamina. But what chiefly distinguishes them 
from the ‘granular’ cells of higher mammals is the possession of a 
particularly developed cytoplasmic cell body. 

Finally it seems worth while to draw attention to the phenomenon of 
satellitosis as observed in the brain of the opossum. When present this 
phenomenon was confined to the deeper laminae and it reached its 
maximum development around the aberrant nerve cells scattered in the 
white matter towards the ependyma of the ventricle. Although the present 
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observations afford no clue to its significance this peculiar limitation should 
be mentioned. 


Summary and conclusions. 


1. The cytology of the cerebral cortex of the opossum (Didelphys 
virginiana) has been examined in NISSL preparations and the findings are 
discussed in relation to some of the general problems of the evolution 
of the cortex. 

2. In all the areas examined cytoarchitectural differentiation is in 
advance of cellular differentiation. The most highly evolved cells are 
found in the area praepiriformis and cornu ammonis; the least differentiated 
are those of the neo-encephalon. The area praeorbitalis is however 
exceptional in that the cells of the 5th lamina are highly differentiated. 
Irrespective of the area the cells of the 5th lamina are invariably the most 
highly differentiated. Both these observations are closely in accord with 
the facts of ontogenetic development. 

3. In general, the cortical areas considered as primitive possess the 
most highly evolved cellular elements, whilst the more recent areas are 
composed of cells showing a primitive organization. 

4. The cells in the 2nd and 4th laminae are not identical in structure 
but both possess an easily visible cytoplasmic body. 


This work was carried out whilst Dr. SmytH held a Rockefeller 
Travelling fellowship. 
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